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Nonequilibrium Phenomena in Quantum Field Theory: 
From Linear Response to Dynamical Renormalization Group 

Shang-Yung Wang, PhD 
University of Pittsburgh, 2001 



This thesis is devoted to studying aspects of real-time nonequiHbrium dynamics in quantum field 
theory by implementing an initial value formulation of quantum field theory. The main focus is 
on the linear relaxation of mean fields and quantum kinetics in nonequilibrium multiparticle quan- 
tum systems with potential applications to ultrarelativistic heavy ion collisions, cosmological phase 
transitions and condensed matter systems. We first study the damping of fermion mean fields in 
a fermion-scalar plasma with a view towards understanding baryon transport phenomena during 
electroweak baryogenesis. We obtain a fully renormalized, retarded and causal equation of motion 
that describes the relaxation of the mean field towards equilibrium and allows an unambiguous 
identification of a novel damping mechanism and the corresponding damping rate. Secondly, we 
apply and extend the renormalization group method to study nonequilibrium dynamics in a self- 
interacting scalar theory, the 0(4) linear sigma model and a hot QED plasma, with the goals of 
constructing a quantum kinetic description that goes beyond usual Boltzmann kinetics and under- 
standing anomalous (nonexponential) relaxation associated with infrared phenomena. Remarkably, 
within this framework the quantum kinetic equations and equations of motion for mean fields have 
the interpretation as the dynamical renormalization group equation which describes the dynamical 
evolution of a multiparticle system that is insensitive to microscopic details. By solving these equa- 
tions, we effectively integrate out fast motion dynamics, and are left with an effective theory for slow 
motion dynamics. As a by-product, the issue of pinch singularities in nonequilibrium quantum field 
theory is resolved naturally in real time from a quantum kinetic point of view. The final part of this 
thesis presents a real-time kinetic analysis of direct photon production from a quark-gluon plasma 
created in ultrarelativistic heavy ion collisions. We show that the direct photon yield is significantly 
enhanced by the lowest order energy-nonconserving processes originated in the transient lifetime of 
the quark-gluon plasma. In particular, transverse momentum distribution of direct photons, which 
features a power law spectrum in the experimentally relevant regime, is proposed as a nonequilibrium 
signature of the quark-gluon plasma. 
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Chapter 1 
Introduction 



1.1 Motivation 

The theoretical study of noncquiUbrium dynamics in quantum multiparticlc systems dates back to 
1961 when Schwinger published a pioneering paper, Brownian Motion of a Quantum Oscillator 
Conventional applications of quantum theory are restricted mainly to calculating transition matrices 
for scattering processes. In his 1961 paper Schwinger showed, for the first time, how quantum 
theory can be formulated to study initial value problems associated with the dynamical evolution of 
nonequilibrium quantum systems. 

Over the past two decades, nonequilibrium dynamics in quantum field theory has attracted a 
great deal of interest as new experimental techniques in particle and condensed matter physics con- 
tinue to probe novel nonequilibrium quantum phenomena that require field-theoretical descriptions. 
Specific examples of current theoretical interest involving nonequilibrium dynamics of quantum fields 
include, just to name a few, inflationary dynamics in the early Universe, electroweak baryogenesis, 
the chiral phase transition and quark-gluon plasma in ultrarelativistic heavy ion collisions, dynamics 
of phase transition in Bose-Einstein condensation and ultrafast spectroscopy of semiconductors. Such 
a diversity of applications reveals the truly interdisciplinary character of nonequilibrium dynamics 
in quantum field theory. 

In this thesis we develop new field-theoretical techniques to study aspects of real-time nonequilib- 
rium dynamics in quantum field theory with a view towards potential applications to ultrarelativistic 
heavy ion collisions, cosmological phase transitions and condensed matter systems. 

Ultrarelativistic heavy ion collisions 

The goal of ultrarelativistic heavy ion collisions is to create and study a new state of hot and dense 
matter in the laboratory. This program began almost two decades ago with the fixed target heavy ion 
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experiments at BNL Alternating Gradient Synchrotron (AGS) and CERN Super Proton Synchrotron 
(SPS). This new phase of matter, the quark-gluon plasma (QGP) [^j-Q, whose fundamental degrees 
of freedom are the deconfined quarks and gluons is a prediction of quantum chromodynamics (QCD) 
and believed to have existed in the early Universe during the first few microseconds after the big 
bang when the temperature was greater than about 200 MeV 10^^ K). 

Recent assessments of the collected data from CERN SPS seem to provide positive evidence for 
the creation of a new state of matter in Pb+Pb collisions based on a multitude of different observa- 
tions, ranging from anomalous dilepton production, strangeness enhancement, J /tp suppression to 
pion interferometry It is fair to say, however, that while the evidence is very suggestive it is far 
from conclusive. Experiments at higher collision energies provided by the newly operating Relativis- 
tic Heavy Ion CoUider (RHIC) at BNL and the planned Large Hadron CoUider (LHC) at CERN are 
expected to allow for a quantitative characterization of the quark-gluon plasma and detailed studies 
of its early thermalization processes and dynamical evolution. 

The Relativistic Heavy Ion Collider at BNL is currently studying Au+Au collisions by colliding 
two ultrarelativistic, highly Lorentz contracted gold nuclei with center-of-mass energy per nucleon 
pair y/sMN ^ 200 GeV. At RHIC energies the two colliding nuclei penetrate one another and most 
of the baryons are expected to be carried away by the receding nuclei (the fragmentation region). 
The quarks and gluons composing the nuclei collide and transfer a large amount of energy from the 
colliding nuclei to the vacuum, creating a nearly baryon-free region of hot and dense matter in the 
form of energetic quarks and gluons that strongly interact with each other. The hard parton-parton 
scatterings thermalize the quarks and gluons on a time scale of about 1 fm/c (^ 3 x 10^^"* s) and 
produce a deconfined and chirally restored quark-gluon plasma in local thermal equilibrium that lives 
long enough to generate detectable signals. This hot and dense quark-gluon plasma (energy density 
e ~ 4 — 6 GeV/fm'^, corresponding to an initial temperature T ^ 200 — 300 MeV) then expands 
rapidly due to internal pressure and cools down to the deconfinement temperature Tqcd ~ 160 MeV, 
below which the quark-gluon plasma undergoes the QCD phase transition (confinement and chiral 
symmetry breaking) and condenses into a gas of hadrons. If the transition is of first order, the quarks, 
gluons and hadrons coexist in a mixed phase before the phase transition is completed. The hadrons 
continue to scatter from one another, maintaining the pressure and causing further expansion and 
cooling. Eventually, the hadron gas becomes sufficiently dilute that scatterings among the hadrons 
cease and the hadron gas freezes out. Estimates based on energy deposited in the central collision 
region at RHIC suggest that the quark-gluon plasma has a transient lifetime about 10 fm/c and the 
overall freeze-out time of order about 100 fm/c. 

Although there is some evidence that local thermal equilibrium may be achieved at the late stages 
of the heavy ion collisions, colliding nuclei are inherently nonequilibrium quantum multiparticle 
system as manifested in copious production of particle at extremely high energy densities and on 
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unprecedentedly short time scales as well as in rapid expansion and cooling. It is therefore a 
theoretical challenge to understand the full noncquilibrium dynamics of ultrarelativistic heavy ion 
collisions, from the initial state of the heavy ion reaction (i.e., the two colliding nuclei) up to the 
freeze-out of all initial and produced particles after the reaction. 

Over the past two decades, our theoretical picture of the quark-gluon plasma is mainly based on 
thermal equilibrium models which neglect most of the noncquilibrium dynamical effects, but make 
simplified assumptions on the very initial stages of the collisions, e.g., thermalization or plasma 
creation. In the last few years, microscopic (kinetic) [^|-^ and macroscopic (hydrodynamical) 
transport models have been developed to describe various stages of heavy ion collisions. On the 
one hand, while these noncquilibrium approaches have provided predictions that are in reasonable 
agreement with the available experimental results at lower collision energies (mostly from AGS and 
SPS), they are nevertheless phenomenological models which are crucially based on ad hoc semiclassi- 
cal assumptions and approximations such that their validity in extreme situations expected to arise 
during the early states of ultrarelativistic heave ion collisions at RHIC and LHC is questionable at 
best. On the other hand, whereas a coherent picture of the collision dynamics is emerging, a mi- 
croscopic field-theoretical description that allows extracting unambiguous signatures of quark-gluon 
plasma formation remains an open question. 

With the first RHIC results of Au+Au collisions at ^s^^ = 130 GeV coming out recently 
[T^ ], in order to correctly interpret of the experimental data there is a pressing need for a field- 
theoretical description of the quark-gluon plasma created in ultrarelativistic heave ion collisions. 
In particular, there may exist interesting genuine noncquilibrium phenomena which could lead to 
unambiguous signatures of quark-gluon plasma formation. An example is provided by the possible 
formation of disoriented chiral condensates (DCC's), which are regions of misaligned vacuum in 
the internal isospin space A microscopic, first-principle description of novel noncquilibrium 

phenomena in ultrarelativistic heave ion collisions will definitely require a thorough understanding 
of noncquilibrium dynamics in quantum field theory, especially QCD, at finite temperature and 
density. 

Cosmological phase transitions 

The QCD phase transition is not the only phase transition expected to occur in the early Universe. 
On the other hand, based on the evidence from observational dada over the past decades it is now 
widely accepted that the Universe that we observe today is a remnant of several cosmological phase 
transitions that took place at different temperature (energy) scales and with remarkable consequences 
at low temperatures [|l^. The idea of phase transitions in cosmology is closely related to that of 
symmetry braking in particle physics, hence the study of cosmological phase transitions play a 
fundamental role in our understanding of the interplay between cosmology and particle physics in 
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extreme environments. 

Current theoretical ideas beyond the standard model suggest that when the Universe was about 
10"'^^ second old there could have been a symmetry-breaking phase transition at the grand unified 
theory (GUT) scale Tqut ~ 10^^ GeV, above which the GUT symmetry is restored and the strong 
interaction is unified with the electroweak interaction. The GUT phase transition is usually associ- 
ated with an important cosmological stage: infiation, i.e., an epoch of accelerated expansion of the 
Universe Originally introduced by Guth in order to explain the initial conditions for the 

hot big bang model, the idea of infiation has subsequently become the cornerstone of modern cos- 
mology. Current observations of the temperature anisotropics in the cosmic microwave background 
(CMB) seem to confirm predictions of the inflationary cosmology: a fiat Universe, an almost scale 
invariant spectrum of density perturbations that are ultimately responsible for large scale structure 
formation |p^ . 

The next symmetry-breaking phase transition in the early Universe is the electroweak phase 
transition, which occurred at the electroweak scale Tew ^ 100 GeV when the Universe was about 
10~^^ second old. The most remarkable cosmological consequence of the electroweak phase transi- 
tion is the possibility that the observed baryon asymmetry of the Universe may be generated during 
this putative first-order phase transition. As pointed out by Sakharov ]20| long ago, a small baryon 
asymmetry may have been produced in the early Universe if three necessary conditions are satis- 
fied: (i) baryon number violation, (ii) violation of C (charge conjugation symmetry) and CP (the 
composition of parity and charge conjugation), and (iii) departure from thermal equilibrium. The 
standard model and its minimal supersymmetric extensions satisfy all three Sakharov criteria for 
producing a baryon excess As a result, over the last decade electroweak baryogenesis has been 
a popular scenario for the generation of the baryon asymmetry of the Universe and there has been 
considerable interest in nonequilibrium dynamics of the electroweak phase transition. 

The last cosmological phase transition, during which the Universe transformed from a quark- 
gluon plasma to a hadron gas, took place at the QCD scale Tqcd ^160 MeV when the Universe 
was about few microseconds old . Because of the relatively low QCD scale, the QCD phase transition 
perhaps is the only cosmological phase transition that can be (and is currently) studied directly with 
terrestrial accelerators in ultrarelativistic heavy ion experiments ]2^ . Based on a first-order phase 
transition scenario, several potential cosmological implications of the QCD phase transition have 
been proposed, e.g., (i) baryon inhomogeneities that might have affected nucleosynthesis, (ii) solar- 
mass scale primordial black holes that could be part of the cold dark matter (CDM), and (iii) strange 
quark nuggets that could also be a component of the cold dark matter [ p3[ |. 

Most of the theoretical investigations on cosmological phase transitions are based on the pictures 
of phase transition in equilibrium that utilize the finite-temperature effective potential, a field- 
theoretical analogy of the equilibrium free energy in statistical mechanics. Although the effective 
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potential is useful in determining equilibrium properties of the phase transition, e.g., the order 
and critical temperature of the phase transition, it is incapable of describing the full dynamics of 
the phase transition, which is crucial to our understanding of phase transition in a cosmological 
setting in which the system was evolving and therefore not in thermal equilibrium. Whereas the 
importance of nonequilibrium aspects of cosmological phase transition was recognized long ago , 
the self-consistent nonequilibrium description of the inflationary dynamics has only been developed 
recently [ p5[ and field-theoretical descriptions of nonequilibrium dynamics of the electroweak and 
QCD phase transitions are still in their infancy. 

Bose-Einstein condensation and ultrafast spectroscopy 

Recent advances in modern technology have made it possible to explore properties of condensed 
matter systems under unusual laboratory conditions such as at extremely low temperatures or on 
unprecedentedly short time scales. Among many novel developments, two of them have stimulated 
immense theoretical interest: (i) the realization of the Bosc-Einstcin condensation (BEC) |2^, pTf 
in dilute atomic gases in which the atoms were confined in magnetic traps and cooled down to 
temperatures of order fractions of microkelvins, and (ii) the observation of the ultrafast spectroscopy 
in semiconductors in which the hot carriers (electrons and holes) are excited by a femtosecond 
laser pulse [p8|,p9|. Remarkably, the common aspects of these two phenomena that receive intense 
theoretical work but still far from clear are those of the nonequilibrium dynamics. 

The theoretical description of the Bose-Einstein condensation in weakly interacting dilute gases 
has a long history and has accounted for many experimental results pTj . However, a full 
microscopic description of the nonequilibrium dynamics in Bose-Einstein condensates such as the 
condensate formation process, the damping of collective excitations and the relevant time scales, has 
not yet been fully developed. As pointed out by Stoof ||3^], the semiclassical Boltzmann equation 
is unable to treat the buildup of coherence, which is crucial for the phase transition to occur. It is 
therefore evident that such a full microscopic description goes beyond Boltzmann kinetics and calls 
for a deeper understanding of nonequilibrium aspects of phase transitions in quantum multiparticle 
systems. 

The theoretical analysis of the relaxation dynamics in optically-excited semiconductors is usually 
based on the semiclassical kinetic equation of the Boltzmann type. A wealth of recent experimental 
results on the ultrafast femtosecond spectroscopy in bulk gallium arsenide (GaAs) demonstrate 
that for time intervals which are short compared to the optical lattice oscillation period (~ 115 fs 
in GaAs), the relaxation dynamics cannot be described by the semiclassical Boltzmann equation 
in terms of completed, energy-conserving collisions. On such extremely short time scales, time- 
energy uncertainty relation comes into play and virtual (off-shell) processes that do not conserve 
energy will result in important quantum coherent effects. Instead, quantum kinetics has to be 
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used in order to account for the quantum coherent nature of electronic states in the band (e.g., a 
well-defined phase relation between electrons and holes), which is completely neglected within the 
semiclassical Boltzmann description. Furthermore, a comprehensive theoretical framework has to be 
able to describe the buildup of coherence, the dynamics of quantum decoherence and the relevant 
time scales. It is indisputable that the ultimate answer can be provided only by a self-consistent 
nonequilibrium quantum dynamical approach. 

This brief description of timely interdisciplinary physics highlights the necessity for a thorough 
understanding of nonequilibrium dynamics in quantum field theory. 

1.2 New Developments in this Thesis 

The study of nonequilibrium dynamics in quantum field theory plays an important role in our under- 
standing of a variety of nonequilibrium phenomena in ultrarelativistic heavy ion collisions, cosmo- 
logical phase transitions and condensed matter systems. In the literature, however, the conventional 
theoretical framework is largely based on finite-temperature thermal field theory and the semiclassi- 
cal Boltzmann description of kinetics. On the one hand, finite-temperature thermal field theory can 
only study static quantities (of a system in thermal equilibrium) such as the time- independent damp- 
ing rates and mean free paths, which at best can provide only estimates of nonequilibrium properties 
near equilibrium but certainly not the full real-time dynamics. On the other hand, whereas Boltz- 
mann kinetics is capable of describing real-time nonequilibrium dynamics, it relies crucially upon 
the validity of the quasiparticle approximation (well-defined "quasiparticles" with a long lifetime) 
and Fermi's golden rule (energy-conserving processes with time- independent transition rates). 

The main theme of this thesis is to provide, from first principles, a theoretical study of nonequi- 
librium dynamics in quantum field theory directly in real time with the emphasis of extracting 
potential experimental signatures. In particular, the following new developments distinguish the 
work presented in this thesis from the existing work in the literature. 

1. Initial value formulation and linear response. We have developed an initial value for- 
mulation in quantum field theory that allows us to obtain fully renormalized, retarded and 
causal equations of motion for nonequilibrium expectation values of quantum fields and to 
study nonequilibrium quantum dynamics in linear response theory directly in real time. 

2. Dynamical renormalization group. The rcnormalization group (RG) is a powerful tool 
to extract the physics of a multiparticlc system that is insensitive to system details. In the 
literature the renormalization group method is usually confined to the study of static and 
equilibrium properties. We have developed a renormalization group method to study real- 
time nonequilibrium quantum dynamics. We have demonstrated explicitly that the dynamical 
evolution of a quantum multiparticlc system corresponds to a dynamical renormalization group 
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flow in real time with equilibrium being a fixed point of the flow, and that the equation of 
motion for the mean field and the quantum kinetic equation are interpreted as the dynamical 
renormalization group equation which describes the slow motion behavior of a nonequilibrium 
system. 

3. Quantum kinetics directly in real time. We have derived from first principles the quan- 
tum kinetic equation for the (quasi) particle distribution function in scalar and Abelian gauge 
field theories that transcends the semiclassical Boltzmann equation — in the sense that medium 
effects, off-shell (energy-nonconserving) processes and infrared threshold divergences are in- 
cluded consistently. 

4. Resolution of pinch singularities. A technical but important issue of pinch singularities 
in nonequilibrium quantum field theory is clarified and resolved directly in real time from 
a quantum kinetic point of view. A real-time kinetic analysis combined with the dynamical 
renormalization group reveals that pinch singularities signal the breaking down of perturbation 
theory in the long-time limit and provides a consistent resummation scheme to extract the 
nonperturbative long-time dynamics. 

5. Phenomenological application and experimental prediction. Most importantly, based 
on a real-time kinetic approach to direct photon production from the QGP, we have shown that 
emission of hard direct photons from a QGP created at RHIC and LHC energies is significantly 
enhanced by energy-nonconserving effects associated with the transient QGP lifetime. In 
contrast to the exponential falloff predicted by the usual equilibrium calculations the transverse 
momentum distribution of direct photons at midrapidity falls off with a power law, providing a 
distinct experimental nonequilibrium signature of the QGP formation in ultrarelativistic heavy 
ion collisions. 

1.3 Outline of the Thesis 

The thesis is organized as follows. Chapter || is a brief introduction to the theoretical framework 
and most of the techniques that we shall be using in this thesis. This includes short reviews of 
the Schwinger-Keldysh closed-time-path formalism in nonequilibrium quantum field theory and the 
initial value formulation in quantum field theory for studying, in particular, the real-time relaxation 
of mean fields and quantum kinetics. 

In Chapter ^ we study the real-time relaxation of the fermion mean field induced by an adiabati- 
cally switched-on external source in a fcrmion-scalar plasma in terms of the initial value formulation. 
The emphasis is on obtaining a fully renormalized, retarded and causal initial value problem for the 
fermion mean field that allows an unambiguous identification of a novel damping mechanism and 
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the corresponding damping rate. 

In Chapter ^ we presents a novel quantum kinetic approach that goes beyond the usual semiclassi- 
cal Boltzmann kinetics by incorporating directly in real time perturbation theory and the dynamical 
renormalization group resummation. Quantum kinetic equation describing the dynamical evolution 
of (quasi) particle distribution functions is derived in a self-interaction scalar field theory and the 
0(4) linear sigma model and further solved in the linearized approximation near thermal equilib- 
rium. We compare the dynamical renormalization group with the familiar renormalization group in 
Euclidean field theory to highlight the equivalence between the two methods. The issue of pinch sin- 
gularities in nonequilibrium quantum field theory is discussed and a real-time resolution is provided 
from the viewpoint of quantum kinetics. 

In Chapter |^ we study the real-time nonequilibrium dynamics in a quantum electrodynamics 
(QED) plasma at high temperature in the relaxation time and hard thermal loop approximations. 
The dynamical renormalization group approach to real-time relaxation and quantum kinetics allows 
us to extract anomalous nonequilibrium dynamics of photon and fermion mean fields associated with 
infrared divergences in gauge field theories at finite temperature. 

This study is phenomenologically relevant to the production of direct photons from the QGP, 
which is the main subject of Chapter |^. There we present a real-time kinetic description of di- 
rect photon production which reveals that production of hard photons is significantly enhanced by 
the transient lifetime of the QGP, hence providing a distinct experimental signature of the QGP 
formation in ultrarelativistic heavy ion collisions. 

Finally, the main results of this thesis are summarized in Chapter |^. A pedagogical introduction 
to the renormalization group (RG) method in studying asymptotic analysis of ordinary differential 
equations is presented in Appendix 

We would like to provide a road map for the reader on how to read this thesis. The reader 
who prefers skimming this thesis without getting into the technical details is suggested to read the 
abstract, the first chapter, the thesis summary chapter and the introduction and conclusions of 
Chapters ^, U, H and |^. Interested reader can continue to read Chapter || and the main text of 
Chapters ^— ^ for all the technical details. 

To conclude this introductory chapter, we emphasize that whereas this thesis focuses mainly on 
nonequilibrium phenomena related to ultrarelativistic heavy ion collisions, the fundamental work on 
real-time nonequilibrium dynamics in quantum field theory presented in this thesis is truly inter- 
disciplinary and can be adapted to study a variety of nonequilibrium quantum phenomena in the 
realm of condensed matter physics, quantum optics, astrophysics and cosmology. 



Chapter 2 



Nonequilibrium Quantum Field 
Theory 

In this chapter we introduce the Green's function formahsm of nonequihbrium quantum field the- 
ory and the initial value formulation in quantum field theory. These two techniques are closely 
related to each other and constitute the basic theoretical framework which is used to study timely 
nonequilibrium problems in this thesis. 

The Green's function formalism to study nonequilibrium phenomena in quantum field theory 
was originally initiated by Schwinger |jl| , Bakshi and Mahanthappa , and later developed further 
by Keldysh |34| and many others in the 1960s. In the literature, this formalism is usually 

referred to as the Schwinger-Keldysh or closed-time-path (CTP) formalism |3^^-p9|]. 

The essential difference between the usual vacuum quantum field theory and nonequilibrium 
quantum field theory lies in the physical quantity that one is interested in. In the former one is 
interested in obtaining the cross section in a scattering experiment between two beams of particles 
with well-defined momenta. The initial state |in) is prepared at remote past and the final state 
|out) is observed at distant future. The cross section for a given process is related to the transition 
amplitude (5-matrix element) (outjin) between the above asymptotically defined in and out states. 
Using the Lehmann-Symanzik-Zimmermann (LSZ) reduction formula, one can express the 5-matrix 
element in terms of the Green's function (vacuum expectation value of time-ordered product of field 
operators), which in turn has a familiar path integral representation and diagrammatic perturbative 
expansion. 

In nonequilibrium quantum field theory, however, we are interested in the expectation values of 
physical observables. Here, the expectation value is used in the statistical sense and is obtained 
by a weighted sum with the density matrix p. Most importantly, we are want to know the time 
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evolution of these expectation values, which will determine the real-time dynamics of the quantum 
field theoretical system under consideration. The usual S'-matrix formalism designed to calculate 
transition amplitudes between the in and out states simply cannot fulfill this purpose, as there is no 
a priori knowledge about the asymptotic out state at distance future. Hence, a physically intuitive 
description of nonequilibrium phenomena is that of the initial value formulation. The closed-time- 
path formalism of quantum field theory is powerful theoretical framework to describe the expectation 
values of physical observables directly in real time, thus providing a general tool for treating initial 
value problems of nonequilibrium multiparticle dynamics in quantum field theory. 



2.1 Closed-Time-Path Formalism 

The most important quantity in statistical mechanics is the equilibrium density matrix, which con- 
tains all the information of the physical system, likewise the most important ingredient in nonequi- 
librium quantum field theory is the density matrix p specified at an initial time io- In the Heisenberg 
picture the full time dependence is contained in the fields and the density matrix p is not time de- 
pendent as it was in the Schrodinger picture. The expectation value of an operator with one time 
argument (Oh(0) defined by 

{Ouit))^TT[pOu{t)]/TTp. (2.1) 

The questions we would like to ask are the following: how to calculate {Oii{t)) and what is its time 
evolution? 

To be specific, let us consider a system described by the time independent Hamiltonian H and 
the initial density matrix p. The time evolution of Oii(t) is determined by the Heisenberg equation 
of motion 

ijOnit)^[Ouit),Hl (2.2) 

whose formal solution is 

0^{t) = U{to,t)OU{t,to), (2.3) 

where O is the operator at time to, which by definition is the same as the time independent 
Schrodinger operator and U{t,to) = g-«^f(t-to) jg fj^^ time evolution operator in the Schrodinger 



picture. Substituting Eq. ( |2.3| ) into Eq. {2A) and using the property U{to,t)U{t,to) = 1, we can 
rewrite the expectation value {Onit)) as 

/.^ TT[pUito,t)OUit,to)] 
^'^«(*)^= Tr[pUito,t)Uit,to)] • 

The numerator and denominator of the above expression have a simple interpretation as the time 



evolution along a closed-time contour. As illustrated in Fig. 2.1, this closed contour starts from to 



to t' along the forward branch C+ and back to to along the backward branch C_, where any point on 
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Figure 2.1: The contour C along the time axis used to evaluate the generating functional for nonequi- 
librium Green's functions. It consists of a forward branch C+ running from to t' and a backward 
branch C_ from t' back to to- The cross denotes insertion of an operator. 

the forward branch is understood at an earlier instant than any point on the backward branch. The 
only difference between the numerator and the denominator is that in the numerator the operator 
O is inserted at time t. The above observation can be easily generalized to the expectation value of 
product of operators with one time argument (Oaiti) 011(^2) • ■ • OuiUi)), thus leading to multiple 
insertion of operator O at times tn, . . ., ^2 and ti. As usual, the insertion of the operator O can 
be achieved by introducing external source coupled to the operator in the time evolution operator, 
constructing the generating functional and taking variational derivatives of the generating functional 
with respective to the source. Note that, however, unlike the usual situation where there is only 
time forward evolution, we now have both forward and backward time evolution. Hence we need 
to introduce two different sources and J^, respectively, for the forward and backward time 
evolution. 

For later convenience in the construction of the nonequilibrium Green's functions, we will intro- 
duce sources coupled to the field operator $ in the time evolution operators.^ The nonequilibrium 
generating functional is defined as 



where U{t,t'; J) denotes the time evolution operator in the presence of the external source J. Note 
that in general J+ 7^ J~, so that Z[J+, J^] depends on two different sources. If these are set equal, 
one has Z[J,J] — Trp, which is equal to unity after proper normalization, being a statement of 
unitarity. 

The functional Z[J^ , J~] can be represented as a path integral by imposing boundary conditions 
in terms of complete sets of eigenstates of the Schrodinger field operator $. Making use of the 
completeness of eigenstates and the path integral representation for the time evolution operators, 
one obtains the following path integral representation for Z[J^ , J^]: 



Here we discuss the case for (real) Bose fields, generalization to the case for Fermi fields is straightforward. 



Z[J+, J-] = Tt[p U{to, t'; J-) U{t\ to; J+)], 



(2.5) 
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X exp 



Jto ^ 



(2.6) 



where (^^ ) refers to field defined on the forward C+ (backward C_) branch and £[<!>] is the 
Lagrangian density. We note that in above functional expression <i>+ and $^ are not independent 
variables, but are linked through the boundary conditions at time t' in the future. To circumvent 
this difficulty, we will take t' ^ oo for all practical purpose and treat $+ and <i>^ as independent 
variables. 

The generating functional Z[J^ , J^] can be written in a compact path-ordered form 

(■02 



where 



and 



Z]Jr] = 



D<j>i D<j>2 {(t>i\p\h) 

$+ for t e C+ 
for t e C- 



I?$c exp 



Jc 



; / dtd^x(^C[^c] + Jc^, 



for t e C+ 
for t e 




dt= dt- 

C J to C+ 



dt 



(2.7) 



(2.8) 



(2.9) 



to C_ 



accounting for the opposite direction of integration along the backward time branch. 

The above results are very general and are valid for any initial density matrix p. In fact, as can 



be seen from Eqs. (2^) and ( p.7[ ), the only effect of the initial density matrix p is to specify the 
boundary conditions for the path integral at to through the matrix element {4ii\p\(j)2) , thus in turn 
to impose the boundary conditions on the noncquilibrium Green's functions. This means that the 
equations of motion for the Green's functions t > are not influenced by the presence of the 
initial density matrix. 

We note that it is always possible to express the matrix element of the initial density matrix as 
an exponential of a polynomial in the fields M,EO|: 



l|p|02) = e 



-Kl 



where 



K 



d^xKa{yi)(l)a{^) + / d^xd^y0a(x)is:ab(x,y)0b(y) 



(2.10) 



(2.11) 



with a,h = 1,2. In the above expression the constant K can be absorbed into the normalization and 
the various coefficient functions Ka{^^)^ Kab{'^,y), etc., have the interpretations of initial conditions 
on the one-point (mean field), two-point correlation functions, etc., which contain all the information 
of the initial density matrix p. For an initial density matrix which is diagonal in the basis of the 
(quasi)particle number operators, it can be shown that all the initial correlation functions higher 
than two point vanish. 
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Perturbation theory and Feynman rules 

The most convenient feature of the CTP formalism of nonequilibrium quantum field theory is that it 
is formally analogous to standard quantum field theory, except for the fact that the fields have contri- 
butions from both time branches. In particular, as in usual field theory, one obtains from variations 
of the generating functional Z[J~^ , J~] (or, equivalently, Z[Jc]) the nonequilibrium Green functions, 
which, however, now include all correlations between points on either forward and backward time 
branches. The path-ordered n-particle (2n-point) Green's function is defined by 

Gcixi ...Xn,x[...x'J = i" {%^Xi) . . . ^XnMx'J . . . <^{x[)}, (2.12) 

with 



(2.13) 



where %■ is the contour ordering operator along the closed contour C. Although for physical observ- 
ables the time arguments are on the forward branch, both forward and backward branches will come 
into play at intermediate steps in a self-consistent calculation. 

For instance, there are four types of single-particle (two-point) Green's functions (i.e., propaga- 
tors) to be given by (here we have suppressed the space arguments for notational simplicity) 

G++{t,t') = ±j(T$(t)$(t')) 

= G>{t,t')e{t - t') + G<{t,t')e{t' - t), 
G—{t,t') = ±z(T$(t)$(«')) 

= G> {t, t')e{t' -t) + G< {t, t')e{t - 1'), 

G-+{t,t') = ±t{<fm{t')) = G>it,t'), 

G+-(t,t') = im')m) = G<{t,t'), (2.14) 

where T (T) is the (anti)time-ordering operator, the upper (lower) sign refers to Bose (Fermi) fields, 
6{t — t') is the Heaviside step function and $ = for bosons but ^ = ^t-yO £qj. Djj-ac fermions. 
These nonequilibrium Green's functions are not completely independent of each other, instead they 
are related by the relation 

G++it,t') + G—{t,t') - G-+{t,t') - G+-{t,t') = 0, (2.15) 

where use has been made of the identity 6{t — t') + 6{t' — t) = 1. Furthermore, we can define the 
retarded and advanced Green's functions, respectively, by 

GR{t,i!) = G++{t,i!)-G+-{t,t') 

= [G>{t,t')-G<{t,t')]e{t-t'), 

GA{t,t') = G++{t,t')-G-+it,t') 

= [G<{t,t')-G>{t,t')]e{t' -t). (2.16) 
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We note that the Wightman functions G^(t, t') determine uniquely the nonequihbrium single-particle 
Green's functions, thus playing important roles in nonequihbrium quantum field theory. In equi- 
librium G{x,x') = G{x — x') and the space-time Fourier transform of the Green's functions have 
specific analytic properties that allow to extract dispersion relations that are important to our later 
discussion. 

One can then construct a diagrammatic perturbative expansion of the Green's functions much 
in the same manner as in usual field theory by decomposing the full Lagrangian into the free field 
(noninteracting) part £o and the interaction part Cint and treating the latter as perturbation. The 
doubling of the sources, fields and integration contours in the nonequilibrium generating functional 



results in an effective nonequilibrium CTP Lagrangian density [see Eq. (2.6)] 

/:cTp[$+, = , (2.17) 

which in turn leads to the following nonequilibrium Feynman rules in perturbation theory: 

1. There are two types of interaction vertices, corresponding to fields defined on the forward 
and backward time branch. The vertices associated with fields on the forward branch are the 
usual interaction vertices, while those associated with fields on the backward branch have the 
opposite sign. 

2. There are four kinds of free field single-particle Green's functions (propagators), corresponding 



to the possible contractions of fields among the two branches [see Eq. ( 2.14 )]. Besides the usual 
time-ordered propagators which are associated with correlations of fields on the forward branch, 
there are antitime-ordered propagators associated with correlation of fields on the backward 
branch as well as the unordered Wightman functions associated with correlation of fields on 
different branches. 

3. The combinatoric factors, rules for loop (time and momentum) integrals, etc., remain the same 
as in usual field theory. 

In this thesis we study two different types of nonequilibrium phenomena in quantum field theoret- 
ical systems: the relaxation of mean fields in the linearized approximation (i.e., linear response the- 
ory) and the quantum kinetic equation that describes the time evolution of the single- (quasi) particle 
distribution function, hence we will specify below in Sec. |2.2| the initial density matrix p used in 
each case. 



2.2 Initial Value Problems in Quantum Field Theory 

The initial value formulation provides a convenient description for studying the time evolution 
of dynamical systems, ranging from Hamilton's equations, kinetic equations and hydrodynamical 
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equations in classical physics to the time-dependent Schrodinger equation in nonrelativistic quantum 
mechanics. The fundamental ingredients of initial value problems are the equations of motion, which 
govern the dynamical evolution, and the initial conditions, which contain all the dynamics prior to 
some initial time to from which on the evolution of the state is followed. The advantage of the initial 
value formulation is that the dynamics of the system for times t > to will be determined once the 
equations of motions and initial condition are prescribed. 

When one studies relativistic quantum field theory, however, one generally does not use the initial 
value formulation, instead a quite different approach to relativistic quantum field theory based on 
5-matrix elements was developed. The main reason for this is that the majority of problems to 
which the theory has been applied do not require knowledge of the detailed time evolution of the 
system. Indeed, as mentioned in the beginning of this chapter, the typical quantity of experimental 
interest is the cross section, which is most easily calculated in terms of the S'-matrix elements. Thus, 
an initial value formulation in relativistic quantum field seems superfluous. 

Nevertheless, as emphasized in the Introduction, recent development in cosmology, ultrarela- 
tivistic heavy ion collisions and condensed matter physics reveals that there are a wide variety of 
physically interesting questions which require not only a quantum field-theoretical description but 
also a thorough understanding of the time evolution of these systems. Clearly, one must depart from 
the usual description in terms of the time-independent S'-matrix element and treat the dynamics 
with the full time evolution. 

The rest of this chapter is devoted to an introduction to two important initial value problems 
for studying nonequilibrium dynamics in quantum field theory: linear relaxation of the mean fields 
and quantum kinetic theory. 

2.2.1 Linear relaxation of mean fields 

The mean fields under consideration are the expectation value of field operators in the nonequilibrium 
state induced by the external time-dependent perturbations. Our strategy to study the irreversible 
relaxation (damping) of these mean fields as initial value problems is to prepare these mean fields 
for a system via adiabatic switching-on of an external perturbation from the remote past. Once the 
external perturbations are switched off at time t = 0, the induced mean fields must relax towards 
equilibrium and we aim to study this nonequilibrium dynamics in linear response theory [^-^3| 
directly in real time. 

To illustrate all the relevant physics and avoid the complexity associated with the issue of ther- 
malization, we assume the system at to = ^oo is described by the free (quasi)particle Hamiltonian 
Hq and in thermodynamical equilibrium at a temperature T. We further assume that the inter- 
action Hamiltonian _ffi,it is adiabatically turned on while we prepare the mean fields via adiabatic 
switching-on the external perturbation. This allows the induced mean fields to be "dressed" adia- 
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batically by the interaction and external perturbation during the preparation time. Hence in this 
case the initial density matrix p at time = — oo has the form p = exp(— /3i7o)/Tr[exp(— /?i7o)], 
where /3 = 1 /T is the inverse temperature. 

The initial value formulation for studying the linear relaxation (damping) of the mean fields 
begins by introducing c-number external sources J{x) coupled to generic (Bose or Fermi) quantum 
fields $(x).0 In the presence of external perturbation the effective CTP Lagrangian density becomes 

/:cTp[$+(a;),$~(a;)] ^£[<S>+{x)] + [j(a;) $+(x) + H.c] - ($+ ^ $"). (2.18) 

The presence of external sources will induce responses of the system. The expectation value of $(x) 
induced by J{x) in a linear response analysis is given by^ 



0(x,<) = (<&±(x,t))j 

= ±f dt' [ d^x'GR{K,t;x',t')J{x',t'), (2.19) 

Jto J 

where the upper (lower) sign in the second equality refers to Bose (Fermi) fields and G/j(x, t; x', i') 



is the retarded Green's function defined in Eq. (2.16) 



Gfl(x,t;x',t') = ±z[($±(x,t)$+(x',t'))-(<I>^(x,i)^"(x',t'))] 

= [G>{x,t;x',t')-G<{x,t;x',t')]9{t~t'). (2.20) 
In the above, [ • , • ]zp denotes commutator (— ) for bosons and anticommutator (+) for fcrmions and, 



as before, $ = $'1^ for bosons but 'I' = ^I't'yO for Dirac fermions. We emphasize that in Eq. (2.20) 
the expectation values are calculated in the CTP formalism with full Lagrangian density C, but in 
the absence of external source as befits the linear response approach. 

As explained above, a practically useful initial value formulation for the real-time relaxation of 
the mean field is obtained by considering that the external source is adiabatically switched on in 
time from to = — oo and suddenly switched off at t = 0, i.e., 

J(x,t) = J(x)e''*6l(-t), £^0+. (2.21) 

For t > 0, after the external perturbation has been switched off, the mean field will relax towards 
its equilibrium value and our aim is to study this relaxation directly in real time. 



^The reader should not confuse the external source ,/ with the source introduced previously in the construction 
of the nonequilibrium generating functional Z[J+, J~]. The external sources J introduced here, just like those 
introduced in usual linear response theory, are external perturbations to the system, hence it is set to be equal on 
both forward and backward time branches. 

^We assume here that in the absence of the external source J{x), the mean field (pix) vanishes identically in thermal 
equilibrium. This is true for fermion and gauge mean fields and for scalar mean filed in the unbroken symmetry phase. 
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For the study of linear relaxation of the mean fields, we assume that the system is in thermal 
equilibrium at a temperature T at time t =^ —oo with 0(x, t = —oo) = 0. The retarded and 
equilibrium (i.e., translational invariant in time) nature of Gj^(x, t; x', t') and the adiabatic switching 
on of J(x, t) entail that 

0(x, t = 0) = <^o(x), </>(x, t < 0) = 0, (2.22) 

where 4>o{'x.) is determined by J(x) [and vice versa, V'o(x) can be used to fix J(x)] and, here and 
henceforth, an overdot denotes time derivative. We note that 0(x, f = 0) is unspecified even though 
(/)(x, t < 0) = 0. This is because the external source is instantaneous switched off at i = 0. In fact 
there could be initial time singularities associated with our choice of the external source given in 



Eq. (2.21). ft has been shown in Ref. |4j] that the long-time behavior is insensitive to the initial 
time singularities, hence we will not address this issue here as it is beyond the scope of this thesis. 

In order to study the relaxation of the mean field ^(x, t) for T > 0, we need to relate the linear 
response problem to an initial value problem for the equation of motion of the mean field. This can 
be achieved by considering the following formal equation of motion for the mean field 

I?o 0(x,t) = ± J(x,i), (2.23) 

where the upper (lower) sign refers to Bose (Fermi) fields and the (intcgro-differential) operator D 



is the inverse of the retarded Green's function Gr. The source J{x) is given by Eq. (2.21) and (/)(x) 



satisfies the condition Eq. ( ^.22 ) at t = 0. For t < 0, on the one hand, Eq. (2.23) is exactly the linear 



response problem Eq. ( ^.1S| ) satisfying the condition Eq. ( ^.22| ) at < = 0; for t > 0, on the other 
hand, it describes the time evolution of the mean field 4'{x, t) with the initial condition specified by 



Eq. (2.22) 



It is at this stage where the nonequilibrium CTP formalism provides the most powerful frame- 
work. The real-time equation of motion for the mean field (pix) can be obtained via the tadpole 
method , which automatically leads to a retarded and causal initial value problem for the mean 
field (j){x). 

The central idea of the tadpole method is to write the field into the c-number expectation value 
plus quantum fluctuations around it, i.e., writing 

<i>±(x,t) = 0(x,t)-fx^(x,t), with (x±(x,t))=0. (2.24) 

The equation of motion for the mean field 0(x, t) can be obtained to any order in perturbation 
theory by requiring the tadpole condition (x^(x, t)) = to all orders in perturbation theory [psf . 

2.2.2 Quantum kinetic theory 

The foundation of modern kinetic theory dates back to 1872 when Boltzmann published his famous 
kinetic equation for dilute gases pq]. The Boltzmann equation is a nonlinear integro-differential 
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equation for the single-particle distribution function, a microscopic quantity introduced by Boltz- 
mann that describes the density of gas molecules in phase space. The success of Boltzmann's kinetic 
description for dilute gasses not only provides the first theoretical understanding of macroscopic 
physics in terms of the underlying microscopic dynamics, but also makes kinetic theory one of the 
most powerful approaches to transport phenomena in macroscopic systems. 

A necessary criterion for the validity of Boltzmann's kinetic description is that there must be a 
wide separation of time scales in the system. A classic example is furnished by the derivation of 
the classical Boltzmann equation: Boltzmann introduced, with great intuition, "molecular chaos" 
(assumption that particles in a dilute gas are not correlated), which crucially depends upon the fact 
the duration of each individual collision process is much shorter than the average time that elapse 
between two successive collisions. Based on Boltzmann's classical picture of completed collisions, one 
can easily write down the semiclassical Boltzmann equation, also known as the Uehling-Uhlenbeck 
equation |^7[| , by modifying the collision term using Fermi's golden rule (which gives the quantum 
mechanical transition probability per unit time) and taking into account the quantum statistics of 
the particles (bosons or fermions). Whereas the semiclassical Boltzmann equation is appealing in 
describing a variety of transport problems, Boltzmann kinetics is an approximation only and needs 
to be justified from first principles, especially for quantum systems under extreme nonequilibrium 
conditions which have been a subject of great theoretical and experimental interest over the last 
decade. 

Although a wide separation of time scales is a necessary condition for the validity of the (semi)- 
classical Boltzmann kinetics, it is far from sufficient. In particular when there is competition of rele- 
vant time scales or quantum transient kinetics at early stages is physically important, a full quantum 
kinetic equation, of which usual Boltzmann kinetics is a semiclassical limit, must be obtained. The 
most widely used method in the derivation of the quantum kinetic equation is that of nonequilibrium 
Green's functions originally developed by Kadanoff and Baym IsH]. This method involves Wigner 



transforms of two-point Green's functions and eventually a gradient expansion |48, 49|. Such gradient 
expansion assumes that the center-of-mass Wigner variables are "slowly varying" , but it is seldom 
clear at this level which are the fast and which are the slow scales involved. A "coarse-graining" 
procedure is typically invoked that averages out microscopic scales in the kinetic description, leading 
to irreversible evolution in the resulting kinetic equation. Such an averaging procedure is usually 
poorly understood and justified a posteriori. 

In this thesis, we propose to study quantum kinetics by focusing on a systematic, first-principle 
derivation of quantum kinetic equations from the underlying quantum field theory. Specifically, we 
aim to provide a kinetic description of an initially prepared nonequilibrium quantum system by 
studying the time evolution of the single-particle distribution functions for the relevant degrees of 
freedom directly in real time. 
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This novel approach begins by defining a suitable (quasi) particle number operator whose expec- 
tation value nk{t) is interpreted as the single particle distribution function for the relevant degrees 
of freedom of momentum k in the presence of the medium, and by separating the total Hamiltonian 
into the free part that commutates with the (quasi)particle number operator and a interaction part 
that describes interaction between these (quasi)particles. As a consequence of the medium effects, 
the above defined (quasi) particles are in general different from the corresponding free particles in 
the vacuum. In the Heisenberg picture, the time evolution of the distribution function dn\^{t)/dt is 
obtained by using the Heisenberg equation of motion. Choosing the initial density matrix p at a 
finite initial time t = toto be diagonal in the basis of the (quasi) particles number operators but with 
nonequilibrium distributions n\^{tQ), one can write down a perturbative expansion for dn\^{t)/dt as 
an initial value problem in terms of the initial distribution function. The solution of this initial 
value problem contains secular terms that grow with time and invalidate the perturbative expansion 
at large times as a consequence of the fact that perturbation theory neglects the change in the 
distribution function as time goes on. 

We introduce a novel dynamical renormalization group that allows a consistent resummation of 
secular terms in real time and leads to an improvement of the perturbative expansion that is valid 
on kinetic time scales. This real-time initial value formulation has the advantage that it displays 
the relevant microscopic and kinetic time scales, leads to a well-defined coarse-graining procedure, 
allows a systematic improvement on the kinetic description and, if necessary, consistently includes 
nonexponential relaxation in a physically intuitive manner. 

To conclude this chapter, we emphasize that the main theme of this thesis is describing the 
dynamical development of a nonequilibrium multiparticle system, which evolves from an initially 
prepared quantum state. Since the nonequilibrium system is characterized by its dynamical evolu- 
tion in time, a physically intuitive description of nonequilibrium dynamics is that in terms of the 
initial value formulation. In this thesis we propose to study real-time relaxation of mean fields and 
quantum kinetics in quantum multiparticle systems as initial value problems. The CTP formalism 
of nonequilibrium quantum field theory is a powerful theoretical framework designed to describe 
directly in real time the expectation values of physical observables in the presence of medium, thus 
providing a general tool for treating initial value problems of nonequilibrium multiparticle dynamics 
in quantum field theory. 



Chapter 3 

Fermion Damping in a 
Fermion- Scalar Plasma 



3.1 Introduction 

The propagation of quarks and leptons in a medium of high temperature and/or density is of fun- 
damental importance in a wide variety of physicaUy relevant situations. In stellar astrophysics, 
electrons and neutrinos play a major role in the evolution of dense stars such as white dwarfs, 
neutron stars and supernovae pO[ |. The propagation of quarks during the nonequilibrium stages 
of the electroweak phase transition is conjectured to be an essential ingredient for baryogenesis 
at the electroweak scale both in non-supersymmetric and supersymmetric extensions of the stan- 
dard model Furthermore, medium effects can enhance neutrino oscillations as envisaged in the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect and dramatically modify the neutrino electromag- 
netic couplings [S^ . 

In-medium propagation of particles is dramatically different from that in the vacuum. The 
medium modifies the dispersion relation of the excitations and introduces a width to the propagating 
excitation |^,54,4^,^ that results in damping of the amplitude of the propagating mode. Whereas 



the propagation of quarks and leptons in a QED or QCD plasma has been studied thoroughly [ p5| , p3[ , 
a similar study for a scalar plasma has not been carried out to the same level of detail. Recently 
some attention has been given to understanding the thermalization time scales of boson and fermion 
excitations in a plasma of gauge |Q and scalar bosons It has also been shown that fermion 
thermalization is an important ingredient in models of baryogenesis mediated by scalars [ p8[ . Most 
of the studies of fermion thermalization focus on the mechanism of fermion scattering off the gauge 
quanta in the heat bath (Landau damping). Although the scalar contribution to the fermion self- 
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energy to one-loop order has been obtained a long time ago p^ , scant attention has been paid to a 
more detailed understanding of the contribution from the scalar degrees of freedom to the fermion 
relaxation and thermalization. As mentioned above, this issue becomes of pressing importance in 
models of baryogenesis and more so in models in which the scalars carry baryon number | |5^ ] . 

In this chapter we focus on several aspects of propagation of fermion excitations in a fermion- 
scalar plasma. Specifically, we offer a detailed and general study of fermion relaxation and thermal- 
ization through the interactions with the scalars in the plasma directly in real time for arbitrary 
scalar and fermion masses, temperature and fermion momentum. More importantly, we focus on a 
novel mechanism of damping for fermion excitations that occurs whenever the effective mass of the 
scalar particle allows its kinematic decay into fermion pairs. This phenomenon only occurs in the 
medium and is interpreted as an induced damping due to the presence of scalars in the medium. It 
is a process different from coUisional damping and Landau damping which are the most common 
processes that lead to relaxation and thermalization. This process results in new thermal cuts in the 
fermion self-energy and, for heavy scalars, a quasiparticle pole structure in the fermion propagator 
that provides a finite width to the fermion excitation. The remarkable and perhaps nonintuitive 
aspect of our analysis is that the decay of the scalar in the medium leads to damping of the fermion 
excitations and their propagation as quasiparticle resonances. 

The effective real-time Dirac equation in the medium allows a direct interpretation of the damping 
of the fermion excitation and leads to a clear definition of the damping rate. By analyzing the 
quasiparticle wave functions we obtain an all- order expression for the damping rate that confirms 
and generalizes recent results for the massless chiral case ]59|,|60|]. 

In order to provide a complementary understanding of the process of induced decay of the heavy 
scalars in the medium and the resulting fermion damping, we study relaxation of the fermion distri- 
bution function by using a semiclassical Boltzmann equation. Linearizing the Boltzmann equation 
near the equilibrium distribution, we obtain the relation between the thermalization rate for the 
distribution function in the relaxation time approximation and the damping rate for the ampli- 
tude of the fermion mean fields to lowest order in the Yukawa coupling. This analysis provides a 
real-time confirmation of the oft quoted relation between the interaction rate (obtained from the 
Boltzmann kinetic equation in the relaxation time approximation) and the damping rate for the 
mean field |5j,^. More importantly, this analysis reveals directly, via a kinetic approach in real 
time how the process of induced decay of a heavy scalar in the medium results in damping and 
thermalization of the fermion excitations. A study of the relation between the interaction rate and 
the damping rate has been presented recently for gauge theories within the context of the imaginary 
time formulation Our results provide a real-time confirmation for the scalar case. 



This chapter is organized as follows. In Sec. 3.2 we obtain the renormalized effective in-medium 
Dirac equation for the fermion mean field in real time starting from the linear response to an external 
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(Grassmann) source that induces the mean field. The renormahzation aspects are addressed in 
detail. We then study in detail the structure of the renormalized fermion self-energy and establish 
the presence of new cuts of thermal origin. In Sec. |3.3| we present a real-time analysis of the evolution 
of the fermion mean field and clarify the difference between complex poles and resonances (often 
misunderstood). An analysis of the structure of the self-energy and an interpretation of the exact 
quasiparticle spinor wave functions allows us to provide an all-order expression for the damping 



rate of the fermion mean fields. In Sec. B.4 we present an analysis of the evolution of the fermion 
distribution function in real time by obtaining a (Boltzmann) kinetic equation in the relaxation time 
approximation. We clarify to lowest order the relation between the damping rate of the fermion 



excitations and the interaction rate of the distribution function. In Sec. 3^ we summarize our 
results. 



3.2 Effective Dirac Equation in the Medium 

As mentioned in the Introduction, whereas the damping of collective and quasiparticle excitations 
via the interactions with gauge bosons in the medium has been the focus of most attention, imder- 
standing of the influence of scalars has not been pursued so vigorously. 

Although we are ultimately interested in studying the damping of fermion excitations in a plasma 
with scalars and gauge fields within the realm of electroweak baryogenesis in either the Standard 
Model or generalizations thereof, we will begin by considering only the coupling of a massive Dirac 
fermion to a scalar via a simple Yukawa interaction. The model dependent generalizations of the 
Yukawa couplings to particular cases will differ quantitatively in the details of the group structure 
but the qualitative features of the effective Dirac equation in the medium as well as the kinematics 
of the thermal cuts that lead to damping of the fermion excitations will be rather general. 

We consider a Dirac fermion 5* with the bare mass Mq coupled to a scalar (f> with the bare 
mass mo via a bare Yukawa coupling j/o- The bare fermion mass could be the result of spontaneous 
symmetry breaking in the scalar sector, but for the purposes of our studies we need not specify its 
origin. The Lagrangian density is given by 

4>] = ^{i^- Mo)"^ + ^ap09^</> - ^ml<P^ - Ci[(P] - yo^(t)'f + + ^V + J<P- (3.1) 

The self-interaction of the scalar field accounted for by the term Cj [(j>] need not be specified to lowest 
order. In the above expression, rj and j are the respective external fermionic and scalar sources that 
are introduced in order to provide an initial value problem for fermion relaxation. We now write 
the bare fields and sources in terms of the renormalized quantities (referred to with a subscript r) 
by introducing the renormahzation constants and counterterms: 

^ = Zl/^^r, 0=^y'</'r, V^Z-'/^rjr, J^Z-'/^Jr, 
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y = yaZV'^Z^/Zy, ml = + m^) /Z^, Mq = {Sm + M) /Z.^. 



(3.2) 



With the above definitions, the Lagrangian can be rewritten as (here and henceforth, we have 
suppressed the subscript r for notational simphcity) 

1 „ . , 1 



1 



(3.3) 



where m and M are the renormahzed masses, and y is the renormahzed Yukawa couphng. The 
terms with the coefficients 



1, Scj, — Zcji - 



1, Sm Z.^Mq - AI, 5rn ^ Z^ml~ m^, 6y 



1, 



and SCj are the counterterms to be determined consistently in the perturbative expansion by choos- 
ing a renormahzation prescription. As it wiU become clear below this is the most natural manner 
for obtaining a fully renormahzed Dirac equation in a perturbative expansion. 

Our goal is to understand the relaxation of inhomogeneous fcrmion mean field V'(X; i) = (^(x, t)) 



induced by external source that is adiabatically switched-on dX t — —00. As explained in Sec. 2.2 



the equation of motion for the mean field can be obtained to any order in perturbative theory via 
the tadpole method Ba] by writing 



*±(x,i) =V'(x,t)+x^(x,t), with (x±(x,t))-0. 



(3.4) 



The essential ingredients for perturbative calculations are the following free real-time Green's 
functions (in momentum space): 
(i) Scalar Propagators 



Go++(k,t,i') 



G,-ik,t,t') 
G>{k,t,t') 



G^ (k, t, t') 6{t - t') + G^ (k, t, t') e{t' - t), 

G^ (k, t, t') e{t' -t) + G^ (k, t, t') e{t - 1'), 

G>{k,t,t'), G+-{k,t,t')^G<{k,t,t'), 



I 

2wk 



G^{k,t,t') = i d^xe-'^-''{(f>{0,t')(f>{x,t)) 



2wk 



riB 



where Wk = vk^ + m^ and nB((jj) — l/(e'^" — 1) is the Bose-Einstein distribution, 
(ii) Fermion Propagators (zero fermion chemical potential) 



(3.5) 



S++{k,t,t') 



5'o>(k, t, t') e{t -t') + 5o< (k, t, t') 0{t' - t). 
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S-+{k,t,t') 



S<{k,t,t') 



S> (k, t, t') 9{t' -t) + S< (k, t, t') 9{t - t'), 
S>{k,t,t'), S+-{-k,t,t')^S<{k,t,t'), 
(vl.(x,t)^(0,t')) 

(7°<:^k - 7 • k + M)[l - npiOji,)] e-'^''(*-*') 
(7"wk + 7 • k - M) nF(cDk) e-'^-^*-*') 



2wk 



2<Dk 



■(7V + 7-k-M) [l-np{Ldi,)]e 



(3.6) 



where o^k = Vk^ + and npito) ~ l/{e^^ + 1) is the Fcrmi-Dirac distribution. We note that 



these free propagators given in Eqs. (3.5) and (|3.6|) are thermal because, as discussed in Sec. 2.2 



the initial state is chosen to be in thermal equilibrium and the interaction is assumed to be turned 
on adiabatically. 

In momentum space, we find the effective real-time Dirac equation for the mean field of momen- 
tum k reads 

[{ifdt - 7 • k - + (S^ {ifdt - 7 ■ k) - 5m] V^(k, t) 

dt' S(k, t - t') 7/;(k, t') = -77(k, i), (3.7) 

where dt = d/dt, E(k,t — t') is the retarded fermion self-energy and 

iP{k,t) = I d^xe''^-''ip{x,t). 



Using the real-time free scalar and fermion propagators given by Eqs. (3.5) and ( ^.6| ), respectively, 
we find to one-loop order, that S(k, t — t') is given by 



S(k, t-t')^ 17° (k, i - i') + 7 • k (k, t-t') + (k, t-t'), 



(3.8) 



with 



S(°)(k,t-t') 



E(i)(k, t-t') 



S(2)(k, t-t') 



-y 



d^q 1 



(27r)3 2wp 



cos[(wp -\-uJc^){t - t')] [1 -I- n_B(wp) - ?^F(^^q)] 

-I- cos[(cjp - cjq)(i - t')] [riBii^p) + nF((I'q)] , 
2 / d^q k • q 



sin[(wp -I- Wq)(t - t')] [1 + fiB(wp) 



(27r)3 2fc2cjp(Ilq 

n_F(<I'q)] - sin[(tjp - tl'q)(t - t')] [nB(cjp) -f nF(wq)] 
d^q M 



sin[(wp + Wq)(i - t')] [1 + riBiujp) - npiusq)] 



(27r)3 2a;pa;q 
sin[(a;p - ujq)(t - t')] [^^(wp) -t- nF(wq)] , 



(3.9) 
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where p = q + k and k = |k|. 

As mentioned before, the source is taken to be switched on adiabatically from t = —oo and 



switched off at i = to provide the initial conditions [see Eq. (2.22)] 

V'(k,< = 0) = ?A(k,0), ip{k,t<0)^0. (3.10) 

Introducing an auxiUary quantity a(k, t — t') defined as 

T,{k,t-t') = dt'a{k,t-t') (3.11) 

and imposing rjlk, < > 0) = 0, we obtain the equation of motion for f > to be given by 

[(i7°at -j-k-M)+S^ {i-f"dt - 7 ■ k) - a(k, 0) - Sm] Mt) 

+ 1 dt' <7^{t-t')^]^{t') ^Q. (3.12) 

"'0 

This equation of motion can be solved by Laplace transform as befits an initial value problem. The 
Laplace transformed equation of motion is given by 

[i-f°s - 7 • k - A-f + (5^ (i7"s - 7 • k) - (5ji/ - cr(k, 0) + s a(s, k)] -0(5, k) 

= [if + z5^7" + 5(s, k)] V(k, 0), (3.13) 

where 0(s,k) and CT(s,k) are, respectively, the Laplace transforms of ?A(k, i) and (T(k, t) 

1*00 1*00 

i;{s,k)= (ite-"*V(k,i), a(s,k)=/ dt e'"* a{k,t), 
Jo Jo 

with Re s > 0. 

3.2.1 Renormalization 

Before proceeding further, we address the issue of renormalization by analyzing the ultraviolet 
divergences of the self-energy. As usual the ultraviolet divergences are those of zero temperature 
field theory, since the finite temperature distribution functions are exponentially suppressed at large 
momenta . Therefore the ultraviolet divergences are obtained by setting to zero the distribution 
functions for the scalar and fermion. 



With Eq. (3.11), one can write cr(k, t — t') as 



a(k, t-t') = ij^a'-"^ (k, t - t') + 7 • ka(i) (k, t-t') + a^^) (k, t - t'). (3.14) 
A straightforward calculation leads to the following ultraviolet divergences 
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where, CT*^*^(s,k) {i = 0, 1,2) arc the Laplace transform of cr'^*) (k, t), A is an uhraviolct momentum 
cutoff, K is an arbitrary renormahzation scale, and "~" denotes only the divergent contribution is 
included. Therefore, the counterterms 6^ and 6m are chosen to be given by 



0^ = — m — + fanite, Oa/ = „ in — 



167r2 K ' 87r2 

and the components of the self-energy are rendered finite 



K 



finite, 



a{s, k) + ij'^Sii, = finite, (T(k, 0) + 7 • kJ^ + (5^/ = finite. 



(3.16) 



(3.17) 



The finite parts of the counterterms in Eq. (3.16) are fixed by prescribing a renormahzation scheme. 
There are two important choices of counterterms: (i) determining the counterterms from an on- 
shell condition, including finite temperature effects and (ii) determining the counterterms from a 
zero temperature on-shell condition. Obviously these choices only differ by finite quantities, but the 
second choice allows us to separate the dressing effects of the medium from those in the vacuum. For 
example, by choosing to renormalize the self-energy with the zero-tempcrature counterterms on the 
mass shell, the poles in the particle propagator will have unit residue at zero temperature. Then, in 
the medium the residues at the finite temperature poles (or the position of the resonances) are finite, 
smaller than one and determined solely by the properties of the medium. Thus the formulation of 
the initial value problem as presented here yields an unambiguous separation of the vacuum and 
in-medium renormahzation effects. 

Hence we obtain the renormalized effective Dirac equation in the medium and the corresponding 
initial value problem for the fermion mean field 



i7°s - 7 • k - M - E(s, k) V(s, k) = [i7° + ct(s, k)] ^(k, 0), 



(3.18) 



where I](s,k) — (T(k, 0) — sa{s,'\<i) is the Laplace transform of the renormalized retarded fermion 
self-energy [see Eq. ( [3.17| )], which can be written in its most general form 



I](s,k) = i7°s£:(°)(s,k) + 7-k5(i)(s,k) +Mf(2)(s^k). 



The solution to Eq. (3.18) is given by 

^{s, k) = i {1 + S{s, k) [7 • k -f Af + E(0, k)] } V(k, 0), 
where S'(s,k) is the fermion propagator in terms of the Laplace variable s 



(3.19) 



(3.20) 



5(s,k) 



i7°s - 7 • k - M - E(s, k) 
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°s[l - f (0) (s, k)] - 7 • k[l + f (1) (s, k)] + M[l + £(2) j^)] 



S2[l - £:(0)(s,k)]2 + fc2[l +f(l)(s^k)]2 -f M2[l +f (2)(s,k)]2 



(3.21) 



We note that the square of the denominator in Eq. ( 3.21 ) is recognized as 



det 



i'j^s - 



7 • k - M - I](s, k) 



(3.22) 
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The real-time evolution of t) is obtained by performing the inverse Laplace transform along 
the Bromwich contour in the complex s-plane parallel to the imaginary axis and to the right of all 
singularities of ip{s,'k). Therefore to obtain the real time evolution we must first understand the 
singularities of the Laplace transform in the complex s-plane. 



3.2.2 Structure of the self-energy and damping processes 

To one-loop order, the Laplace transform of the components £''-'(s,k) (i — 0, 1,2) of the fermion 
self-energy S(s, k) [see Eq. ( [3.19| )] can be written as dispersion integrals in terms of spectral functions 
p«(fco,k) 



£(i)(s,k) > = 



+ 00 



dkf] 



p(")(fco,k) 








w 




A:op(2)(fco,k) 




I M J 



(3.23) 



with the one-loop spectral functions given by the expressions 



p(°)(fco,k) = -y" 



1 



(27r)3 2tjp 



(5(fco - t^p - t^q) [1 + nB(tJp) - riF(Wq)] 



p(i)(fco,k) 



p(2)(fco,k) 



-|-J(fc0 - ^^p +^q) [^J-B(Wp) + n-_F(Wq)] 

k • q 



-y 



(27r)3 2fc2tJp(I)q 



J(fco - Wp - ojq) [1 -f ns(t^p) - ?^F(Wq)] 



- (5(fco - t^p + t^q) [n-B(t^p) + J^Fl^q)] 
d^g 1 



6{ko - LUp - ujq) [1 + riBiujp) - n_F(wq)] 



(3.24) 



- S{ko -ujp+ ujq) [nsi^^p) + n-_F(wq)] . 

The analytic continuation of the retarded self-energy E(s, k) and its components f*^'^ (s, k) in the 
complex s-plane (physical sheet) are defined by 



E(tj,k) = i;(s = -icj + 0+,k) = Rei;(tJ,k) +ilml](a',k), 
£«(tj,k) = £«(s = -iw + 0+,k) = Re£(*)(w,k) +ilm£«(w,k). 



Using the relation 



1 



P 



where p denotes the principal value, one finds 



(3.25) 



(3.26) 



Re£(o)(s,k) 
Re5(i)(s,k) 
Re£(2)(s^k) 



dko 



p(o)(fco,k) 






fcopW(fco,k) 


H 




fcop(2)(fco,k) 

J 




I A/ j 



(3.27) 



and 



Im£(")(a;,k) = sgn(u;) U") (|^|, k) + (-|c.|, k)' 

^ a; L 
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ImfW(a>,k) = |sgnH[p«(|c.|,k)-pW(-|c.|,k) 



(3.28) 



We note that the real (imaginary) parts of the retarded self-energy are even (odd) functions of cu. 

The denominator of the analytically continued fermion propagator in Eq. ( p. 21 ) can be written 
in a compact form tj^ — cj^ — n(cLi, k) where 



Il(uj,k) 



'^2£(0) (cu, k) + (^, k) + Af2£(2) (^, k) 

-Lj^ [£^°Hu;,k)f + k^ [£(i)(^,k)]2 + Af2 [£^^\cj,k)f 



(3.29) 



We recognize that the lowest order term (but certainly not the higher order terms) of this effective 
self-energy can be written in the familiar form |5l,B2lE3| 



(3.30) 



U{uj,k) = -2L^£^°\uj,k) + e£'^^\uj,k) + M^£'^^\cu,k) 
= -iTr[(7%-7.k + M)S(c^,k)] . 



The imaginary part of Il{Lu,k) evaluated on the fermion mass shell will be identified with the 



fermion damping rate (see below). The expression given by Eq. ( 3.30 ) leads to the familiar form of 
the damping rate, but we point out that Eq. ( |3.30 ) is a lowest order result. The full imaginary part 



must be obtained from the full function n(a>,k) in Eq. (3.29) and the generalization to all orders 
will be given in a later section below. 

For fixed M and m, the delta function constraints in the spectral functions p'^*-' (w, k) can only be 
satisfied for certain ranges of u. Since the imaginary parts are odd functions of lo we only consider the 



case of positive uj. The delta function 6{\llj\ — ujp — Cjq) has support only for > y/k'^ + {m + Af)^ 
and corresponds to the usual two-particle cuts that are present at zero temperature corresponding 
to the process f ^ s + f. Hence, the two-particle cuts (both for positive and negative uj) do not give 
a contribution to Imn(ti;,k) on the fermion mass shell and the only contributions are from terms 



proportional to nB{wp) + nF{ujq) in Eq. (3.28). We find 

,2 I '^^1 "-b(^p) + '^F(Wq) 



Imn(a;,k) = -vry sgn(cj) 



(27r) 



2WpWq 



(|cj|wq - k • q - 



X 5{\uj\ - UOp+ tJq) + (|ci;|cjq + k • q -|- A/^) 5{\uj\ -I- Wp - CJq) 



(3.31) 



In the above equation the first delta function 5{\lo\ — tjp + Wq) determines a cut in the region 



\u)\ < y/k'^ + {m — A/)2 and originates in the process s f + f , whereas the second delta function 
(5(|a;| + Wp — Wq) determines a cut in the region < \uj\ < k and originates in the process f + s f . 
We note that the first cut originates in the process of decay of the scalar into fermion pairs and 
the second cut (w^ < k'^) is associated with Landau damping. As to be shown below, both delta 
functions will result in restriction in the range of the momentmn integration in Imn(a;,k). 
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For m > 2M the scalar can kincmatically decay into a fermion pair, and in this case the fermion 



pole is embedded in the cut < + [m — M)^ becoming a quasiparticle pole and only the first 
cut contributes to the quasiparticle width. This is a remarkable result, the fermions acquire a width 
through the induced decay of the scalar in the medium. This process only occurs in the medium 
(obviously vanishing at T = 0) and its origin is very different from either collisional broadening or 
Landau damping. A complementary interpretation of the origin of this process as a medium induced 
decay of the scalars into fermions and the resulting quasiparticle width for the fermion excitation 



will be highlighted in Sec. iA using the kinetic approach to relaxation. 

The width is obtained to lowest order from Imn((I'k,k), and the on-shell delta function (5(tl'k — 
+ ^q) is recognized as the energy conservation condition for the process s f + f. To lowest 

order we find the following expression for Imn((Dk,k) for arbitrary scalar and fermion masses with 

m > 2M and arbitrary fermion momentum and temperature: 



Imn(a;k, k) 



'Try 



d^q Wk^^q — k • q — 



(27r)3 



2cjqCJp 



In 



1 + e^l^'^'i 



q=q- 



where 



2M2 



k 1 



2M" 



± 



4M2\ 



(3.32) 



(3.33) 



with q (z {q ,q'^) being the support of (5(ciJk 



for fixed k. 



3.3 Real-Time Evolution of Mean Fields 

The real time evolution is obtained by performing the inverse Laplace transform as explained above. 
This requires analyzing the singularities of '!/'(s,k) given in Eq. ( |3.21 ) in the complex s-plane. It is 



straightforward to see that the putative pole at s = has vanishing residue, therefore the singularities 
are those arising from the fermion propagator S{s, k). For m < 2Af , the fermion poles at lu — ujp{k) 
are real and isolated away from the multiparticle cuts along the imaginary axis s — —iuj with 



< uj-ik) and > a;+(fc), where uj±{k) — ^k"^ + (m ± Af)^. In this case, the inverse Laplace 
transform can be performed by deforming the Bromwich contour, circling the isolated poles, and 



wrapping around the cuts as depicted in Fig. 3.1 (a). One finds the inverse Laplace transformation 



is dominated by the pole contribution and, to lowest order, the fermion mean field ■(/'(k, i) oscillates 
at late times with a constant amplitude and frequency Lo^iji). 

On the other hand, for m > 2M (i.e., when the scalar particle can decay into fermion pairs) the 
fermion poles become complex and are embedded in the lower cut, and one must find out if they 
are complex poles in the physical sheet (the domain of integration) or moves off to the unphysical 
(second) sheet. 
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Figure 3.1: The analytic structure of the fermion propagator S{s, k) and the complex contour used 
for evaluating the mean field tp(\i,t) for (a) m < 2M and (b) m > 2M. 



3.3.1 Complex poles or resonances 

For TO > 2M the fermion poles become complex and are embedded in the cut s = —icu with 
|w| < w_(fc). The position of the complex poles are determined from the zeros of — (I>^ — Tl{u>, k) 
in the denominator of the analytically continued fermion propagator for lu ~ i-Up{k) — ij{k), where 
LiJp{k) ['y{k)] is the real (imaginary) part of the complex pole. In the following discussion, we consider 
the narrow width approximation, \"f(k)/LUp(k)\ <C 1, where the physical concept of a propagating 
mode is still meaningful. 



With the expressions for the discontinuities in the physical sheet given by Eq. (3.25), the position 
of the complex pole is determined by 



[cjp(fc) -i7(/c)]2 -(I>2 _Ren(wp(fc),k) +isgn[7(A:)]Imn(wp(/c),k) = 0. 



(3.34) 



Retaining terms at most linear in ^{k), one finds the real and imaginary parts of Eq. (3.34), respec- 
tively, become 

ojlik) =ujI + ReU{ujp{k),k), (3.35) 

and 



7(fc) = sgn[7(fc)] 



Imn(wp(/c),k) 

2LUp{k) 



(3.36) 



Since Ren(ti;,k) is an even function of uj, here and henceforth, we choose ujp{k) to be the positive 



solution of Eq. (3.35). To lowest order, one finds 

Ren(wk,k) 



Wp(fc) 



2ujk 



cDk + — Tr[(70wk - 7 • k + M) ReS(cDk, k)]. 
4wk 



(3.37) 



To this order, the solution for j{k) is obtained by approximating LUp{k) « Wk hi Eq. (3.3i 
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A close inspection, however, shows that Eq. ( ^.36| ) cannot have a solution because Inin(a-', k) 
is an odd function of lu and Imn((I'i(;, k) < 0. Therefore, there is no complex pole in the physical 
sheet. Indeed, this is a fairly well-known (but seldom noticed) result: if the imaginary part of the 
self-energy on the mass shell is negative, then there is no complex pole in the physical sheet and the 
pole has moved off into the unphysical (second) sheet. 

In the case that Imn(tJk,k) > 0, complex poles appear in the physical sheet, but in such case 
there are two poles with both signs for 7 (A:) (one corresponds to a decaying exponential in time 
and the other a growing exponential in time), which is the signal of an instability, not of damping. 
However since we confirm that Imn([Iik, k) is negative in the case under consideration, the complex 
poles are in the unphysical sheet and correspond to resonances. 



3.3.2 Scalar decay implies fermion damping 

We have shown that there are no complex poles in the physical sheet in the complex s-plane (the 
integration region), hence the only singularities are the cuts along the imaginary axis s = —iuj with 
\uj\ < w_(fc) and \lu\ > cj+(fc). The contour of integration can be deformed to wrap around these 
cuts as depicted in Fig. (b). Since the resonance is below (and away from) the multiparticle 
cut > Lj+(fc), in the narrow width approximation and consistent with perturbation theory the 
contribution from the cut < Lu^{k) becomes the dominant one while that from the multiparticle 
cut is always perturbatively small. 



It proves convenient to write the product in Eq. (3.20) in the compact form 

g(g,k) [7-k-KAf + S(0,k)j ^ 9 "^l^^'^L (3-38) 



n(s,fc)' 



which defines A/'(s,A:), and to change variables to s = —ioj ± 0+ on the right-hand (+) and left- 
hand (— ) side of the cut. After some algebra, one finds the following contribution to the real time 
evolution of the mean field 



TT 



Col - Ren(w, k)]2 + Imn(u;, k)2 



ImAf (lu, k) [uj"^ -Ldl- Ren(tj, k)] 
;2 - (1-2 _ Ren(tj, k)]2 + Imn(w, k) 



V'(k,0). (3.39) 



The term proportional to ReA/'(a;, k) features a typical Breit-Wigner resonance near the real part 
of the complex pole, where uj^ — uj^ — nij(wp,k) = 0, since, for m > 2M, the imaginary part of 
the self-energy at this value of to (perturbatively close to ±cjk) is nonvanishing. On the other hand, 
the term proportional to ImA/'(w,k) is a representation of the principal part in the limit of small 
Ren(ti;,k) and is therefore subleading. 

The sharply peaked resonances at uj = ±ujp(k) dominate the integral and give the largest con- 
tribution to the real-time evolution of ip(k,t). In the limit of a narrow resonance, the uj integral 
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is performed by taking the integration limits to infinity and approximating near the resonances at 

LO — ±UJp{k) 

1 Imn(tj,k) Z{k) -f{k) 



where 



Z{k) = 



1 - 



9Ren(w, k) 



1 -1 



7(fc) 



Zik) 



Imn(tjp(fc),k). 



(3.40) 



(3.41) 



The wave function renormahzation constant Z{k) is finite, since the effective self-energy n(aj, k) has 
been rendered finite by an appropriate choice of counterterms. Only when the counterterms are 
chosen to provide a subtraction of the self-energy at the position of the resonance uj w ±ujp{k) will 
result in Z{k) = 1. On the other hand, as noted above, if the counterterms are chosen to renormalize 
the theory on the fermion mass shell at zero temperature then Z{k) describes the dressing of the 
medium and is smaller than one. As will be shown shortly, the width of the Breit-Wigner resonance 
7(/c) is the damping rate of the fermion mean field. 

The integration in the variable u can be performed under these approximations (justified for 
narrow width), leading to the real-time evolution 
Z{k) 



^{k,t) 



2ul{k) 



RcAA(cjp(fc), k) e-^^p^*^)* + ReAA(-cjp(fc), k) e*"^^'^)* e-'^('=)* V'(k, 0). (3.42) 



Gathering the results for Imn(Ljp(fc), k) given in Eq. ( 3.32| ), we find the damping rate for the fermion 
mean fields at one-loop order to be given by (for m > 2M) 



327rfctDk 



AM' 



In 



1 + e-'^"^" 



(3.43) 



q=q 



where q are given in Eq. (3.33) 



Figures 3.2 and 3.3 display the behavior of ^(k) for several ranges of the parameters. We have 
chosen a wide range of parameters for the ratios of the scalar to fermion masses (to/A/) and the 
temperature to fermion mass (T/M) to illustrate in detail the important differences. The damping 
rate features a strong peak as a function of the ratio k/M. This peak is at very small momentum 
when the ratio of scalar to fermion mass is not much larger than 2, but moves to larger values of 



the fermion momentmn when this ratio is very large. Figure 3.3 displays this feature in an extreme 



case (to/A/ = 800) to highlight this behavior. The height of the peak is a monotonically increasing 
function of temperature as expected. This is one of the important results of this work: the decay 
of the heavy scalar into fermion pairs results in a induced damping of the amplitude of fermionic 
excitations. 



3.3.3 All-order expression for the damping rate 

Since there is no complex pole solution in the physical sheet for to > 2A/, there are no solutions 
of the effective in-medium Dirac equation. However, we can define the spinor wave function of the 
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Figure 3.2: The fermion damping rate plotted as a function of the momentum for m/M — 4 and 
T/M = 150 (soHd), 100 (dotted), and 50 (dashed). 



resonance by considering the solutions of the in-medium Dirac equation with only the real part of 
the self-energy at a; = ±tjp(fc) |5^,^. The form given by Eq. ( 3.1£ ) for the self-energy is general and 



not restricted to perturbation theory, hence our analysis below is valid to all orders in perturbation 
theory. 



Using the expression for the self-energy given by Eq. ( 3.1£ ) for s = ^ia'p(/c)-|-0+, and introducing 
the following variables 

W = c^p(fc)[l-Ref(°)(cc>p(fc),k)], 
K = k[l + Re£(i)(wp(fc),k)], 

M = M[l + Ref(2)(cjp(fc),k)], (3.44) 

we find the resonance wave functions ^'s(±Wp(fc), k) obey the following effective in- medium Dirac 
equation 

(±7V - 7 • /C - X) *,(±u;p(fc), k) = 0, (3.45) 

where s = 1,2 labels the spin. Furthermore, using the properties that Re£^*''(w,k) are even 
fmictions of uj and depend on fc, we can define two in-medium particle (i.e., positive energy) 
solutions C/s(k) = ^'s(wp(fc),k) and two in-medium antiparticle (i.e., negative energy) solutions 
Ki(k) — ^'s(— ij-'p(fc), — k). These in-medium Dirac spinors satisfy 

{j°W-j-K-M) Us{k) = 0, 

(7°W - 7 • a: + X) (k) = 0, (3.46) 
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Figure 3.3: The fermion damping rate plotted as a function of the momentum for m/M — 800 and 
T/M = 1000 (sohd), 800 (dotted), and 600 (dashed). 



and can be constructed from the usual free Dirac spinors through the replacement Wk W, k ^ /C, 
and M ^ M. 

Combining Eqs. ( 3.19 ), (|3.25 ) and ( 3.2£ ), we obtain, from the definition of the damping rate given 
by Eq. ( |3.4C| ), the following expression for the width of the resonance to all orders in perturbation 
theory: 

Z{k) 



lik) 



4wp(fc) 
MZjk) 

'2LUp{k) 



Tr [(7">V -j-K + M) ImE(wp(fc), k)] 
[t/s(k)ImE(a;p(fc),k)t/4k)] , 



(3.47) 



where we have alternatively written the expression for the width in terms of the exact resonance 
spinors in the medium. This result confirms those found in Refs. |5^, and leads to the often 
quoted expression for the width in lowest order ||54,43|. 



3.4 Kinetics of Fermion Relaxation 

To clarify and confirm independently the result of the previous section that decay of the scalar leads 
to a quasiparticle width of the fermionic excitations, we now provide an analysis of the relaxation of 
the fermion distribution function via a kinetic Boltzmann equation. This analysis will also provide, 
directly in real time, a firm relationship between the damping rate and the interaction rate in the 
relaxation time approximation (i.e., linearization near equilibrium). 

Because here we are only interested in the kinetics of fermions, we assume that the scalars 
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are in thermal equilibrium. Let us denote the distribution function for fermions of momentum k 
and spin s by fis,k(0- the kinetic approach, the time derivative of this distribution function is 
obtained from the Boltzmann equation. Since for a fixed spin component the matrix elements for 
the transition probabilities are rather cumbersome, we study the spin-averaged fermion distribution 
function defined as n-k,{t) — ^ nk,s(i)- 

For m > 2M, two processes are responsible for the change in the fermion populations in a 
fermion-scalar plasma: s f + /, which provides the "gain" term in the balance equation, and 
f + f ^ which provides the "loss" term. Using the standard approach to obtain the semiclassical 
Boltzmann equation, we find 

-77 "kW = Try / ——3 — ^ d{uJk + LUq~UJp) 

at J (277)'' ZWqWkWp 

X KK)[1 - rlk(i)][l - nq(t)] - [1 + nBK)]nkW^qW] , (3.48) 

where p = k + q. It is easy to check that the above equation has an equilibrium solution given by 
^k(^) = "F('^k) for all momentum k. 

Let us now consider small departure for fermions from this equilibrium solution and study the 
relaxation to equilibrium in the relaxation time approximation. This is achieved by introducing a 
disturbance in the distribution function for fermion of (fixed) momentum k so that 

nk{t) =nFiiUk)+Snk{t), (3.49) 



where i5nk(i)/"-F(^k) <C 1. Retaining terms at most linear in 5nk(i), one obtains from Eq. ( |3.48 ) 
an equation for 5nk(t) 

^Snk{t)^^r{k)dMt), (3.50) 
at 

where r(fc) is the interaction rate, whose inverse characterizes the time scale for the fermion distri- 
bution to approach equilibrium 

r(k) = ny jTT-^ [nB(Wp) + nF(wq)J5(wk-Wp + Wq) 



y^TO^T / 4M2 



1677fca;k 



1 TT- In 



1 + 



q=q+ 



(3.51) 



with given in Eq. ( 3.33| ) 



Upon comparing the interaction rate with the damping rate found in the previous section [see 
Eq. ( 3.43| )], we provide a real-time confirmation of the result 



r(fc)=27(fc). (3.52) 



The kinetic analysis confirms that the damping of the ferniionic quasiparticle excitations in the 
medium is a consequence of the decay of the heavy scalar. Furthermore, this analysis is carried out 
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directly in real time and clearly establishes the relation between the interaction rate in the relaxation 
time approximation and the damping rate of the mean field at least to lowest order. 

Recently a detailed investigation between the damping and interaction rates for chiral fermions 
in gauge theories at finite temperature has been reported in Ref. | |6l[ | within the framework of the 
imaginary-time (Matsubara) formalism of finite-temperature field theory. Our analysis provides a 
complementary confirmation of this result in real time both for the relaxation of the mean field and 
that of the spin-averaged fermion distribution function. 

3.5 Conclusions 

In this chapter we have focused on studying the propagation of fermion excitations in a fermion- 
scalar plasma, as a complement to the more studied issue of fermion propagation in a hot (Abelian 
or non- Abelian) plasma. Our motivation was to provide a real-time analysis of the propagation of 
fermions that could eventually be used in other problems such as neutrino oscillations in the medium 
and noncquilibrium processes in electroweak baryogenesis. 

The first step of the program is to obtain the effective in-medium Dirac equation for the fermion 
mean field directly in real time. This is achieved by relating the problem of linear response to an 
initial value problem for the fermion mean field induced by an external Grassmann source that is 
adiabatically switched on from remote past and switched off at time t = 0. The resulting effective 
Dirac equation for the mean field is fully renormalized, retarded and causal, thus allowing a direct 
study of real-time relaxation of the mean field. We study in detail the structure of the renormalized 
fermion self-energy and establish the presence of new cuts of thermal origin in the fermion propagator. 
We found that when the scalar mass is large enough that decay of the scalar into fermion pairs is 
kinematically allowed, this decay process results in an induced damping of the fermion excitations 
and their propagation in the medium as quasiparticle resonances. Solving the effective Dirac equation 
by Laplace transform, we obtained the time evolution of the fermion mean field, which leads to a clear 
identification of the damping rate. We calculate the damping rate in the narrow width approximation 
to one loop order for arbitrary values of the fermion and scalar masses (provided the scalar is heavy 
enough to decay), temperature and fermion momentum. 

An all-order expression for the fermion damping rate (in the narrow width approximation) is 
obtained from an analysis of the structure of the self-energy and the exact quasiparticle spinor wave 
functions. A kinetic approach based on a semiclassical Boltzmann kinetic equation for the spin- 
averaged fermion distribution function reveals that the interaction rate (obtained in the relaxation 
time approximation) is twice the damping rate of the mean fields at least to lowest order in the 
Yukawa coupling. We emphasize that this relation is established from a real-time analysis both for 
the mean field and the distribution function. 
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Although the expression for the damping rate was obtained for arbitrary values of the scalar and 
fermion masses, temperature and fermion momentum, a deeper analysis is required if the theory 
undergoes a second-order (or weakly first-order) phase transition. The reason being that if the 
fermion masses are a result of spontaneous symmetry breaking in the scalar sector, near a second- 
order (or weak first-order) phase transition both the scalar mass and the chiral breaking fermion 
mass vanish. In this case the kinematic region in momentum space for which the energy conservation 
delta functions are fulfilled shrinks and one must understand if, for soft fermion momentum k ~ yT, 
a resummation akin to the hard thermal loops (HTL's) in hot gauge field theory is required |]63| . 



Chapter 4 



Dynamical Renormalization Group 
Approach to Quantum Kinetics 



4.1 Introduction 

The ultrarelativistic heavy ion experiments at the BNL Relativistic Heavy Ion CoUider (RHIC) 
and the forthcoming CERN Large Hadron CoUider (LHC) have the potential of providing clear 
evidence for the formation of the quark-gluon plasma (QGP) j^,^ and hopefully to study the chiral 
phase transition p^ . Since this is perhaps the only opportunity to study phase transitions that 
are conjectured to occur in particle physics with earth-bound accelerators, a great theoretical effort 
parallel to the experimental work has been devoted to understanding the signatures of the QGP and 
the chiral phase transition. 

A essential part of this theoretical program is to establish, from first principles, a consistent quan- 
tum kinetic description of transport phenomena in a hot nonequilibrium multiparton system. Such 
a kinetic description has the potential of providing a detailed understanding of experimental signa- 
tures such as photon and dilepton production, charmonium suppression, strangeness enhancement, 
freeze-out of hadrons, and collective flow |^,|6^ that will lead to an unambiguous determination of 
the formation of the QGP and the observables of chiral phase transitions in ultrarelativistic heavy 
ion collisions. The kinetic description is also of fundamental importance in the understanding of 
the emergence of hydrodynamics in the long- wavelength limit of quantum fleld theories | |6^ . During 
the last few years, there have been significant advances ranging from first-principle derivations of 
kinetic and transport equations in QCD |67,^, and scalar field theories [^-73|, numerical 



codes that describe the space-time evolution of heavy ion collisions such as the screened perturbative 
QCD (pQCD) parton model [^|7^, the parton cascade model Q and the hydrodynamical transport 
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model [g to lattice simulations of nonequilibrium dynamics in quantum field theories |7^,[76[ and clas- 
sical gauge theories ||77[-[79|. Furthermore, an effective quantum kinetic equation has been advocated 
as a description of the dynamics of soft collective excitations in hot non-Abelian plasmas |8^-^,Q. 

The typical approach to derive quantum kinetic equation involves a Wigner transform of the 
two-point Green's functions and a gradient expansion [^^,67, 68|,|8^ (a gauge covariant Wigner 
transform in the case of gauge theories) and often requires a quasiparticle approximation p8|,|68|,[73|, 
[85| ]. The gradient expansion assumes that the center-of-mass Wigner variables are slowly varying, 
but it is seldom clear at this level which are the fast and which are the slow scales involved. A 
"coarse-graining" procedure is typically invoked that averages out microscopic scales in the kinetic 
description, leading to irreversible evolution in the resulting kinetic equation. Such an averaging 
procedure is usually poorly understood and justified a posteriori. The quasiparticle approximation 
consists in neglecting the broadening of the single-particle states when computing the collision terms. 
A derivation of quantum kinetic equation for a hot QCD plasma along these lines has recently been 
reported recently in Ref . [Q , however the collision terms obtained in the quasiparticle approximation 
turn out to be infrared divergent. 

The goal of this chapter is to provide a novel and alternative derivation of quantum kinetic equa- 
tions directly in real time from the underlying quantum field theory by implementing a dynamical 
renormalization group (DRG) resummation. The strategy to be followed is a generalization of the 
dynamical renormalization group method introduced in Ref. pit] to study anomalous relaxation in 



Abelian gauge theories, but adapted to the description of quantum kinetics |71, 8^. The starting 
point is the identification of the distribution function of the quasiparticles which could require a 
resummation of medium effects (the equivalent of hard thermal loops p^]). The equation of motion 
for this distribution function is solved in a perturbative expansion in terms of nonequilibrium Feyn- 
man diagrams. The perturbative solution in real time gives rise secular terms that grow rapidly 
with time and invalidate the perturbative expansion beyond a particular time scale (recognized a 
posteriori as the kinetic time scale). The dynamical renormalization group implements a systematic 
resummation of these secular terms and the resulting renormalization group equation is the quantum 
kinetic equation. 

The validity of this approach hinges upon the basic assumption of a wide separation between the 
microscopic and the kinetic time scales. Such an assumption underlies every approach to a kinetic 
description and is generally justified in weakly coupled theories. Unlike other approaches in terms 
of a truncation of the equations of motion for the Wigner distribution function, the main ingredient 
in the approach presented here is a perturbative diagrammatic evaluation of the time evolution of 
the distribution function in real time JzH improved via a renormalization group resummation of the 
secular divergences. 

An important bonus of this approach is that it illuminates the origin and provides a natural 
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resolution of pinch singularities [^7 , 88 1 found in perturbation theory out of equilibrium. The pertur- 
bative real-time approach combined with the renormalization group resummation reveal clearly that 
these are indicative of the nonequilibrium evolution of the distribution functions. In this framework, 
pinch singularities are the manifestation of secular terms. 



This chapter is organized as follows. In Sec. 4.2 we study in detail the dynamical renormalization 
group approach to quantum kinetics in the familiar case of a self-interacting scalar theory, including 
in addition a nonequilibrium resummation akin to the hard thermal loops to account for the effective 

we 



scalar masses in the medium and therefore the relevant microscopic time scales. In Sec. 4.3 



discuss the main features of the dynamical renormalization group method and compare it to the 
more familiar renormalization group in Euclidean field theory. 



In Sec. 4.4 we apply these techniques to obtain the kinetic equations for pions and sigma mesons 
in the 0(4) linear sigma model in the chiral limit. In relaxation time approximation we obtain the 
relaxation rates for pions and sigma mesons. This case allows us to highlight the power of this the 
dynamical renormalization group approach to study threshold effects on the relaxation of resonances, 
in particular the crossover between two different relaxation regimes as a function of the momentum 
of the resonance. This aspect becomes phenomenological important in view of the recent studies by 
Hatsuda et al. ]89| ] that reveal a dropping of the sigma mass near the chiral phase transition and an 
enhancement of threshold effects with potential observational consequences in ultrarelativistic heavy 



ion collisions. In Sec. 4.5 we discuss the issue of pinch singularities found in calculations in the so- 
called real-time formalism of nonequilibrium field theory and provide a resolution via the dynamical 



renormalization group. In Sec. 4.6 we summarize our results and discuss further implications. 



4.2 Self-Interacting Scalar Theory 

The starting point is the following Lagrangian for a hot scalar theory: 

m^\{d,q^f~lm^cb'~^_cb\ (4.1) 

where m and A are the physical parameters defined in the vacuum (or, at zero temperature). Since 
we focus on medium (finite-temperature) effects, in the following discussion we ignore the vacuum 
(zero-temperature) ultraviolet divergences that can be absorbed into the renormalization of m and 
A. 

The first step towards understanding the kinetic regime is the identification of the microscopic 
time scales in the problem. In the medium, the bare particles are dressed by the interactions 
becoming "quasiparticles" . One is interested in describing the relaxation of these quasiparticles. 
Thus the important microscopic time scales are those associated with the quasiparticles and not 
the bare particles. If a kinetic equation is obtained in some perturbative scheme, such a scheme 
should be in terms of the quasiparticles, which already implies a resummation of the perturbative 
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expansion accounting for the medium effects. This resummed perturbative expansion is precisely 
the rationale behind the hard thermal loop (HTL) program developed by Braaten and Pisarski [ ]62[ 
in hot gauge field theory and also behind the self-consistent treatment in nonequilibrium scalar field 
theory proposed by Lawrie | |7C[ |. 

Since we are interested in obtaining a kinetic equation for the quasiparticle distribution function, 
the most natural initial state corresponds to an initial density matrix p prepared at time t = 
to that is diagonal in the basis of free quasiparticles number operator (defined below) but with 



nonequilibrium quasiparticle distribution functions nk(io) (see Sec. 2J1). In this picture the width of 
the quasiparticles arises from their interaction and is related to the relaxation rate of the distribution 
function in the relaxation time approximation. This point will become more clear below. 

With such a choice of initial density matrix, perturbative expansions are carried out with the 
following nonequilibrium free quasiparticle Wightman functions (in momentum space): 

G> (k, t, t') = + «k(io)] e-^--^*-*') + nk(to) e^^-^*-*')" , 

Go<(k,t,i') = ^ Lk(io)e-*"''^*-*'^ + [l + nk(to)]e^"-(*-*'^l , (4.2) 

ZWk L J 



where Wk = -^/k^ + m^jj is the energy of the free quasiparticle and m^s is the quasiparticle effective 
mass to be determined self-consistently below. We note that it is easy to check the relation 

Go>(k,t,t') = G<(k,t',t), (4.3) 

which will be useful in following calculations. 

In the self-interacting scalar field theory the HTL resummation can be implemented by writing 
m in the Lagrangian as 

m^^m^g-^W, (4.4) 

where meg is the effective quasiparticle mass which enters in the effective propagators for free 
quasiparticles [see Eq. (4.2)], and 5m? is a counter term which will cancel a subset of Feynman 



diagrams in the perturbative expansion and is considered part of the interaction Lagrangian. The 
effective mass rricff will in general depend on the state of the system (e.g., distribution functions of 
the scalars) through the medium effects. In particular, in thermal equilibrium TOcff depends on the 
temperature and is called the thermal mass | |4^ . 

We note that one could also introduce counter term for the coupling constant and proceed with 
a resummed perturbative expansion in terms of the in-medium effective propagators and vertices, 
akin to those in the HTL program for hot gauge theory However, to lowest order that we 

are interested in only the mass counter term is relevant. 

The implementation of this resummation scheme to lowest in perturbation theory is achieved by 
requiring that the mass counter term 5m? cancel the one-loop tadpole contribution to the scalar 
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self-energy |90 . This leads to the following self-consistent gap equation for m^fj I? 



(4.5) 



where nk(io) is the initial (nonequilibrium) distribution of the quasiparticles at time t — to. We 
emphasize that rrics depends on the initial quasiparticle distribution functions. This observation 
will become important later when we discuss the time evolution of the quasiparticle distribution 
functions and therefore of the effective mass. 

Assuming that the system is in thermal equilibrium, one can make an assessment of the (thermal) 
effective mass TOoff(T). Upon substituting the replacement of equilibrium distribution functions 
nu{to) nB{y/k'^ + TO^jj(r)) in Eq. (4^), one finds that for T ^ mcff{T) the solution of the gap 
equation is given by [po|,|9l|] 



^-Lm,s(^T)+0{mlsiT) In 



"T-cff(T) 

T 



(4.6) 



In particular, for T ^ a/AT 3> m (justified in the weak coupling limit), one can neglect the mass 
m and obtains 



Ar2 



l + OiX 



l/2^ 



(4.7) 



We note that this effective mass determines the important time scales in the medium but is not the 
position of the quasiparticle pole (or, strictly speaking, resonance). 

In the massless case, m.cs{T) serves as an infrared cutoff for the loop integrals hence 



the leading term of Eq. (4.7) provides the correct microscopic time scale at high temperature. 
Furthermore, when the temperature is much larger than the vacuum (zero temperature) mass to, 
the self-consistent resummation is needed to incorporate the physically relevant microscopic time 
and length scales in the perturbative expansion, which, in the weakly coupling limit, is given by 

tmic ~ l/TOoft(r) ~ 1/VXT. 



4.2.1 Quantum kinetic equation for quasiparticles 

Having explained the notion of quasiparticles in the medium and the microscopic time scales asso- 
ciated with them, we now proceed to obtain the kinetic equation for the quasiparticle distribution 
functions. Here we present an alternative to the usual derivation in which correlation functions are 
written in terms of relative and center-of-mass space-time coordinates and a Wigner transform in 
the relative coordinates is performed. In this approach the kinetic equation is obtained in a gradient 
expansion assuming that the dependence on center-of-mass coordinates is weak. Here, however, we 
will not assume such a situation, but instead analyze which are the relevant time scales over which 
a coarse graining procedure must be implemented. 

In the presence of the medium the quasiparticles will have an effective mass rrics resulting from 
medium effects, much in the same manner as the temperature dependent thermal mass in the 
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equilibrium situation discussed above. This effective mass rrics will be very different from the vacuum 
mass m and must be taken into account for the correct assessment of the microscopic time scales. 
In the kinetic approach, this is achieved by writing the Hamiltonian of the theory as H = Hq + Hi^t, 
where 

with n(x, t) — d(f)(ii, t) / dt being the conjugate momentum of 0(x, i). In the above definition, the 
free Hamiltonian Hq describes free quasiparticles of effective mass mcff and is diagonal in terms 
of free quasiparticle creation and annihilation of operators a^(k) and a(k), respectively, and the 
mass counter term has been absorbed in the interaction Hamiltonian -ffint- The pole mass of the 
quasiparticles will acquire corrections in perturbation theory, but these will remain perturbatively. 
With this definition, the lifetime of the quasiparticles will be a consequence of interactions. In this 
manner, the nonequilibrium equivalent of the hard thermal loops (in the sense that the distribution 
functions are non-thermal), which in this theory amount to local (momentum independent) terms, 
have been absorbed in the definition of the effective mass. This guarantees that the microscopic 
time scales are explicit in the quasiparticle Hamiltonian. 

As discussed above, we consider that the initial density matrix p at time t = ta is diagonal in the 
basis of free quasiparticles number operator, but with nonequilibrium initial distribution functions 
nk(io)- Hence, the Heisenberg quasiparticle field operators at time t (in momentum space) can be 
written as 

0(k,t) = -^[a(k,O+at(-k,t)] , U{k,t)^i^[aH-k,t)-a{k,t)], (4.9) 
where a'' {k, t) and a(k, t), respectively, are the creation and annihilation operators in the Heisenberg 



picture and = y^k^ + m^g as defined in Eq. (12) above. The expectation value of Heisenberg 



quasiparticle number operators nk(i) is interpreted as the time dependent quasiparticle distribution 
function and can be expressed in terms of the field 0(k, t) and the conjugate momentum H(k, t) as 
follows 

nk(i) = (at(k,t)a(k,t)) 

= ^[(H(k,t)n(~k,i))+cc-^(</,(k,i)0(-k,t))]-i, (4.10) 

where the bracket ( • ) denotes the nonequilibrium expectation value with the initial (Gaussian) 
density matrix specified by the initial distribution functions rtk(to)- 
The interaction Hamiltonian Hi„± in momentum space is given by 

A f-X dh 



(q,t)0(-q,i). (4.11) 



1^1 (2^)3 

dm? f (fq 
(2^ 
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Taking the derivative of rtk(i) with respect to time and using the Heisenberg field equations, one 
finds 



A 



12wk 

5m? 

2tJk 



([(/)3(k, t)] n(-k, t)) + (n(k, t) [03(k, i)]) 
(^(k,<)n(-k,t)) + (n(k,i)0(-k,t)) 



where we use the compact notation 
[</)3(k,t)]^ 



n '^(qi' ^) '^^(k - qi - q2 - qs). 



(4.12) 



(4.13) 



In a perturbative expansion care is needed to handle the conjugate momentum and the scalar field 
at the same time because of the Schwinger terms. This ambiguity is avoided by noticing that 



(n(k,t)[03(_k,t)]) = lmi^Tr[[(/)3(_k,<)]+pr(k,t')] 
= lm|;([0^(-k,i)]+r(k,t')), 



(4.14) 



where we used the cyclic property of the trace and, as usual, the ± superscripts for the fields refer 
to field defined in the forward (+) and backward (— ) time branch in the CTP formalism. 

We now use the canonical commutation relation between 11 and <j) and define the mass countert- 
erm 5m? — A A/6 to write the above expression as 

A 



Ihn ^ (^2 <)]+ 0- (_k, t')) + A [(0+(k, t)r (-k, i!)) 



12wk 

-(0+(k,i')0"(-k,i)) 



3i 



(0+(q,t)r(-q,i)) 



(4.15) 



The expectation values on the right-hand side of Eq. (4.15) can be obtained perturbatively in weak 
coupling expansion in A. Such a perturbative expansion is carried out in terms of the nonequilibrium 



Feynman rules with the free quasiparticle Wightman functions given by Eqs. (4.2). 

Before proceeding further, an important point to notice is that the free propagators entering in 
the perturbation expansion are the resummed propagators, which include the proper microscopic 
scales as the contribution of the hard thermal loops have been incorporated consistently by summing 
the tadpole diagrams in the scalar self-energy. As a result, the terms with A are required to 
cancel the tadpole contribution in Eq. ( 4.15| ) to all orders in perturbation theory. Thus after A is 
properly chosen in order to cancel the tadpole diagrams in the scalar self-energy, from the formidable 
expression given in Eq. ( [4.15 ) only the first term remains (with the understanding that no tadpole 
diagrams being included) . Physically, this requirement guarantees that the mass in the propagators 
is the effective mass that includes the microscopic time scales. 

To lowest order the condition that the tadpole contributions are canceled leads to the following 
condition on A 

(fik . . 
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Figure 4.1: The Feynman diagrams that contribute to the quantum kinetic equation for a self- 
interacting scalar theory up to two-loop order. 



which in turn entails that the effective mass nics is determined self-consistently by the gap equation 



given in Eq. (4^). We highlight that the requirement that the term proportional to A in Eq. ( 4.15 ) 
cancels the tadpole contributions is equivalent to the HTL resummation in the equilibrium case [po| , 
[9l| ] and makes explicit that iries is a functional of the initial nonequilibrium quasiparticle distribution 
functions nk(^o)- 



At order A, the right-hand side of Eq. ( l.K ) vanishes identically. This is a consequence of the 
fact that both the free quasiparticle Hamiltonian Hq and the initial density matrix p are diagonal in 
the basis of free quasiparticles number operator. Hence, the relaxation of quasiparticle distribution 
functions is solely due to the interaction between quasiparticles. 

To order , we find that the final form of the kinetic equation is given by 



—nu(t)^ — lim (U^(k,t)]+ (j)-(-k,t')) , 

dt ' t'^t dt' XL*^ ^ ' ^ ' ^/ ' 



(4.17) 



dt ' Bwk dt' 

with the understanding that no tadpole diagrams being included as they are automatically canceled 



by the terms containing A and the last term in Eq. (4.15). To two-loop order, the Feynman diagrams 



that contribute to Eq. (4.17) is displayed in Figure 4.1 and the resultant time evolution of the 
quasiparticle distribution function is given by 



dt 



6wk 



7Vi(io) cos[Awi(i ~t")]+ 37V2(io) cos[Aw2(i - t")] 
-|-3AA3(io)cos[Acj3(i - t")] +7V4(io)cos[Acj4(i - t")] 



(4.18) 



where 



Awi = t^k + t^qi + t^qa + t^qa : = ^k + CJqj + CJq^ - W, 



Awa = o^k — t^qi — 



q2 ^ "^qa ' 



l^UJi = Wk ^ iJ^qi — 



q2 "^qa ' 



(4.19) 



and 



M2{t) 



[1 + n^{t)] [1 + nq, {t)] [1 + 71q, (<)] [1 + 71q3 (t)] - n^{t) U^, (t) Hq, (<) 71q3 (t) , 
[1 + nk(i)] [1 + «qi it)] [1 + nq^ (<)] f^qa (t) - «k(<) «qi (t) ^q, (t) [1 -|- U^^ (i)] , 
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U3{t) = [1 + nk(i)] nqi(t) nq,{t) [1 + nq3(0] - ^kW [1 + ?^ql (*)][! + "q^ W] '^q,, W> 

AA4(t) - [1 + nk(i)] nqi(t) nq. W nq3(i) - nk(i) [1 + nq,(t)][l + nq,(i)][l + nq,{t)]. 

(4.20) 



The kinetic equation (4.18) is retarded and causal. The different contributions have a physical 
interpretation in terms of the 'gain minus loss' processes in the plasma. The first term describes the 
creation of four particles minus the destruction of four particles in the plasma, the second and fourth 
terms describe the creation of three particles and destruction of one minus destruction of three and 
creation of one, the third term is the scattering of two particles off two particles and is the usual 
Boltzmann term. 

As will be discussed in detail below, such a perturbative expansion will be meaningful for times 
t <C ikin, where ^kin is the kinetic time scale for the nonequilibrium distribution functions. For 
small enough coupling we expect that ikin will be large enough such that there is a wide separation 
between the microscopic and the kinetic time scales that will warrant such an approximation (see 
below) . 

Since the propagators entering in the perturbative expansion of the evolution equation are 
in terms of the quasiparticle distribution functions at the initial time to, the time integration in 



Eq. ( 4.2C ) can be done straightforwardly leading to the following equation: 

m ^ ,, sin[(u;-^k)(^-to)] 

-T:nk(t) = — dujn[uj,k;Afj{to)] -, r , (4.21) 

at 6 7r(a; - Wk) 

where 

7^[c.,k;AG■(^o)] = ^ / n (2nr2c.^. {2^)'S'{k - qi - q2 - qa) [^(to) S{Au;,) 

+ 3Af2{to) 6{Auj2) + iAfsito) SiAcu^) + Af^ito) <5(Acj4)] , (4.22) 

with j = 1,...,4. 

We are now ready to solve the kinetic equation derived above. Since 7?.[aj, k; A/'j(to)] is fixed at 



initial time to, Eq. (4.21) can be solved by direct integration over t, thus leading to 

n^it) ^ n^{to) + — dujTZ[uj,k;J^j{to)] ^. 4.23 

It is noted that the above expression gives the time evolution of the quasiparticle distribution 

function to lowest order in perturbation theory, but only for early times. To make this statement 

more precise, let us consider the following approximation 

1 — cos[(w — ijJk)(i — io)l / s X 
lim ^ = {t-to)6{oj^LUk). (4.24) 

t-to— »oo Tr(uj — Wkj 

Thus, provided that 7?.[aj, k; A/'j(to)] is finite at = ojk, one finds that nk(t) at large times is given 

by 

nk(i) = nk(to) + '^['^k, k; A/'j(to)] (^ ~ ^o) + nonsecular terms, (4-25) 
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where the nonsecular terms denotes terms that are bound at aU thnes. In the above expression, 
the term that grows (hnearly) with time is referred to as a secular term. The approximation above, 
replacing the oscillatory terms with resonant denominators by tdiuj — ujjs) is the same as that invoked 
in ordinary time-dependent perturbation theory leading to Fermi's golden rule in elementary time- 
dependent perturbation theory . 

Clearly, the presence of secular terms in time restricts the validity of the perturbative expansion 
to a time interval t — to <^ tkin with 

iM„(fc) « ,,J''fir(f (4.26) 



Since the time scales in the integral in Eq. ( 4.23 ) are of the order of or shorter than t^ic ^ 1/m-cff 



the asymptotic form given by Eq. ( 4.25 ) is valid for t — to ^ ^mic- Therefore for weak coupling there 



is a regime of intermediate asymptotics in time defined by 

imic<i-io<ikin(fc), (4.27) 

such that perturbation theory is valid and the corrections to the quasiparticle distribution function 
is dominated by the secular term. 

We note two important features of this analysis: 



(i) In the intermediate asymptotic regime (4.27) the only explicit dependence on the initial 



time ^0 is in the secular term, since TZlcdu, k; A/'j(to)] depends on only implicitly through the initial 
distribution functions. These observations will become important for the analysis that follows below. 



(ii) The fmiction(al) TZ[uJk^ k; A/'j(to)] given by Eq. (4.22) vanishes if the initial distribution func- 
tions are the equilibrium ones as a result of the energy conservation delta functions and the equi- 
librium relation 1 + nsii^) — e^'^nB{i^j)- In this case there are no secular terms in the perturbative 
expansion. 

To highlight the significance of the second point above in a manner that will allow us to establish 
contact with the issue of pinch singularities in a later section, we note that the secular term in 



Eq. (4.25) corresponds to the net change of quasiparticles distribution function in the time interval 



t ~ to- To see this more explicitly, we note that 7?.[ciJk, k; A/'j(io)] can be rewritten as 



A i 
—n[uJi,,k;J^j{to)] = — [l + nu(io)]S<(c^k,k) -nk(io)S>(c^k,k) , (4.28) 



where S>(cjk,k) — I]<(wk,k) — 2 i Iml]fl(wk, k) with Ei?(a;k,k) the on-shell retarded self-energy 
calculated to two- loop order in terms of the initial distribution functions rik(to) l l7l[ . Thus, the 
absence of secular term for a system in thermal equilibrium [for which nk(io) = nsiuJk)] is a 
consequence of the Kubo-Martin-Schwinger (KMS) condition for the equilibrium self-energy 

k) = e^'^^<{Lj, k). (4.29) 
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Furthermore, it can be recognized that the first (second) term in Eq. (4.28), corresponds to the 
"gain" ("loss") part in the usual Boltzmann type collision term. Hence, one can readily interpret 
A^7?.[wk, k; A/j(to)]/3 as the net production rate of quasiparticles per unit time. 

4.2.2 Dynamical renormalization group equation 

Secular divergences are an ubiquitous feature in the perturbative solution of differential equations 
with oscillatory behavior. A resumniation method based on the idea of the renormalization group 
(RG) that improves the perturbative solutions of differential equations was introduced by Goldenfeld 
and collaborators to study pattern formation in condensed matter systems.]^ As explained in 
details in Refs. ||9^,Q, this powerful method not only allows a consistent resummation of the secular 
terms in the perturbative expansion but also provides a consistent reduction of the dynamics to the 
slow degrees of freedom. In a recent development, this method has been generalized to the realm of 
nonequilibrium quantum field theory as a dynamical renormalization group to resum the perturbative 
series for the real time evolution of nonequilibrium expectation values |^,|9^. This generalization 
is a major step since in nonequilibrium systems the equations of motion for expectation values are 
nonlocal and, as will be shown below, require a resummation of Feynman diagrams. 



Having recognized that a time scale emerges in the perturbative expansion (4.25) due to the 
presence of the secular term, we now implement the dynamical renormalization group resummation 
of secular divergences to improve the perturbative expansion following the formulation presented in 
Ref. @. 

To lowest order, this resummation is achieved by introducing the "renormalized" initial distri- 
bution fmictions ?ik(''") which are related to the "bare" initial distribution function 7ik(io) via a 
renormalization constant Zk(T, tg) given by 

nk(to) = 2k(r,io)nk(r), Zk(r, io) = 1 + y 2k(T, io) + O(A^), (4.30) 

where r is an arbitrary "renormalization scale" and Z]s^{T,to) will be chosen to cancel the secular 



term on a time scale r. Substitute Eq. (4.30) into Eq. (4.25), to order we obtain 

nk(i) = "k(T) + y [zi,{T,ta)nu{T) + {t - to)n[uJk,'k;J\fj{T)]] + OiX'^). (4.31) 



Thus, to this order, the choice 



.Uf\ U f ^ ^K,k;AA,■(r)] 

z\^{T:to) = -{t - to) 1—^ (4.32) 

nk(r) 



leads to 

A2 



nk(i) = nk(r) +'—{t-T) 7^K, k; AG-(t)] + ©(A^). (4.33) 



^For a pedagogical introduction to the renormalization group method, see Appendix ^ 
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Whereas the original perturbative solution was only valid for times such that the contribution from 
the secular term remains very small compared to the initial distribution function at time tg, the 



renormalized solution Eq. (4.33) is valid for time intervals t — t such that the secular term remains 
small, thus by choosing r arbitrarily close to t we have improved the perturbative expansion ]9^,p^. 

To find the dependence of nk(T) on r, we make use of the fact that the "physical" distribution 
function rtk(t) does not depend on the arbitrary scale r as a change in the renormalization point 
T is compensated by a change in the "renormalized" initial distribution functions nk(T). This r- 
independence leads to the dynamical renormalization group equation, which to lowest order is given 

by 

(4.34) 



d 

^nk(r) - y 7^[L^k,k;AG■(r)] = 0. 



This renormalization of the distribution function also affects the effective mass of the quasipar- 
ticles since rn^g is determined from the self-consistent equation (4.5), which is a consequence of 



the tadpole cancelation consistently in perturbation theory. Since the effective mass is a functional 
of the distribution function it will be renormalized consistently. This is physically correct since 
the in-medium effective masses will change under the time evolution of the distribution functions. 



Choosing the arbitrary scale r to coincide with the time t in Eq. (4.34), we obtain the resummed 
quantum kinetic equation: 

d^qi d^q2 d^qs 



—n (t) = — 

dt^^^ ' 6tJk J (27r)32wq, (27r)32wq, {2ttY2uj^^ 



(27r)3^3(j^_q^_q^_q^) 



5 (Wk + t^qi + ^Cl2 + t^qa ) A/"! (t ) + 3 (5(tJk + t^qi + t^qa - ^qa ) -^^2 {t) 



-|-3(5(wk - 



^q2 



+ tJqg ) 7V3 (t ) + (5 (tJk - t^qi - t^qs - ^qa ) J^iit) 



(4.35) 



where the Mj(t) are given in Eq. (4.20). 

To avoid cluttering of notation in the above expression, we have not made explicit the fact that 
the frequencies Wq — \/ + m^jj depend on time through the time dependence of mcs, which is 
in turn determined by the time dependence of the quasiparticle distribution function. Indeed, the 
dynamical renormalization group resummation leads at once to the conclusion that the cancelation 
of tadpole terms by a proper choice of A requires that at every time t the effective mass is the 
solution of the time- dependent gap equation 



A f d^q n^{t) 



(27r)3 ^q^+mi^(t)-' 



(4.36) 



where nq(t) is the solution of the kinetic equation (4.35). Thus, the full quantum kinetic equation 
that includes a nonequilibrium generalization of the hard thermal loop resummation in this scalar 



theory is given by Eq. ( [4.35 ) with the replacement a 
self-consistent time-dependent gap equation ( 4.36| ). 



iOqit), where (jJq{t) are the solutions of the 
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The quantum kinetic equation (4.35) is therefore more general than the famihar semiclassical 
Boltzmann equation for a scalar field theory in that it includes the proper in-medium modifications 
of the quasiparticle masses. This approach provides an alternative derivation of the self-consistent 
method proposed in Ref. [7C| ]. 

It is now evident that the dynamical renormalization group systematically resums the secular 
terms and the corresponding dynamical renormalization group equation extracts the slow evolution 
of the nonequilibrium system. 

For small departures from equilibrium the time scales for relaxation can be obtained by lineariz- 
ing the kinetic equation ( [4.35 ) around the equilibrium solution at t = <o- This is the relaxation 
time approximation which assumes that the quasiparticle distribution function for a fixed mode of 
momentum k is perturbed slightly off equilibrium such that nk(io) — nsii^k) + (5rik(^o)i while all 
the other modes remain in equilibrium. 

Recognizing that only the delta function S{ujk — '-^m ~ + ^qa) that multiplies the scattering 
term 7V3(t) in Eq. ( p!35|) is fulfilled, one finds that the linearized kinetic equation reads 

j^5nkit) = -T{k)6nkit), (4.37) 

where T{k) is the scalar interaction rate 

"^^'^ ^ ^ J (2.)32c.,, (2.)32a;,. {2.^2.,, ^''^ - qi - - qa) 

X (5(wk - - -f Wqg ) [1 + Hb (tJqi )] [1 + Hfi {(^q^ )] n-B (Wqj ) 

-"-B(WqJnB(LJqJ [1 +nB(LJq3)] , (4.38) 



where Wk = + m^jj(r). Solving Eq. (4.37) with the initial distribution ^nk(io)i one finds that 
the quasiparticle distribution function in the relaxation time approximation evolves in time in the 
following manner 

5nk(i) = 5n^{to) g-^^f*-*"), (4.39) 

which gives the time scales for relaxation close to thermal equilibrium. 

In the case of soft momentum (k <C mofr(T) <C T) and high temperature (AT^ 3> rn^), one 
obtains 

1 32^24^ 

*-(^"«)^f(A:^":A^- (^-^o) 

For very weak coupling (as we have assumed) , the kinetic time scale is much larger that the micro- 
scopic one iniic ~ l/jTioff ^ since ikin/^mic ^ 1/A S> 1. This verifies the assumption of the 
separation of microscopic and kinetic time scales in the weak coupling limit. 
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4.3 Comparison to the Euclidean Renormalization Group 

In the section we establish a very close relationship between the dynamical renormalization group 
discussed above and the usual renormalization group in Euclidean field theory by demonstrating that 
the rationale behind the two methods are identical and that there exists an one-to-one correspondence 
between them. 

First we review the idea of the Euclidean renormalization group in zero-temperature field theory. 
Let us consider a (massless) scalar field theory in four dimensions with an upper momentum cutoff 
A described by the Lagrangian density 

m^l{d,<l>f-^cl^\ (4.41) 

where Aq is the bare coupling. As an example, we calculate the two particles to two particles 
scattering amplitude in perturbation theory. The one-particle- irreducible (IPI) four-point function 
(i.e., two particles to two particles scattering amplitude) at an off-shell symmetric point to one-loop 
order in Euclidean space is given by 

r(4)(p,p,p,p) ^ Ao - hi In - + 0{Xl), (4.42) 
z p 

where p is the Euclidean four-momentum. Clearly perturbation theory breaks down for A/p > e^^^° . 
Let us introduce the renormalized coupling constant at a momentum scale k as usual as 

Xo^Zx{K)XiK), 2xiK)^l + zi{K)XiK)+OiX^), (4.43) 

and choose zi{k) to cancel the logarithmic divergence at an arbitrary renormalization scale k. Then 
in terms of A(k) the scattering amplitude T^'^'>{p,p,p,p) becomes 

r(4)(p,p,p,p) = XiK) + h'i^) + 0{X% (4.44) 

2 K 

with r^*) (k, K, K, k) — X{k). The scattering amplitude does not depend on the arbitrary renormal- 
ization scale K and this independence implies k dV^^^p, p, p, p) / Ok — 0, which in turn to lowest order 
leads to the renormalization group equation 

^^^/j(A), f3{X)^h\K) + 0{X'), (4.45) 
dK 2 

where /3(A) is referred to as the renormalization group beta function. Solving this renormalization 
group equation with an initial condition X{p) = X that determines the scattering amplitude at some 



value of the momentum and choosing k = p in Eq. ( 4.44|) , one obtains the renormalization group 



improved scattering amplitude (at an off-shell symmetric point) 

T^^\p,p,p,p;Xip))^Xip), (4.46) 
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with X{p) the sohition of the renormahzation group equation 

Hp) = ^ • (4.47) 

l-(3A/2)ln(p/p) 

The connection between the renormalization group in momentum space and the dynamical renor- 
mahzation group in real time (resummation of secular terms) used in the previous section is imme- 
diate through the identification 



io^ln(A/p), i^ln(p/p), r <^ ln(K/p), 



which when replaced into Eq. (4.42) illuminates the formal equivalence with secular terms. 

This simple analysis highlights how the dynamical renormalization group does precisely the same 
in the real time as the renormalization group in zero-temperature Euclidean field theory. Much in 



the same manner that the renormalization group improved scattering amplitude given by Eq. (4.46) 
is a resummation of the perturbative expansion, the kinetic equation obtained from the dynamical 
renormalization group improvement represents a resummation of the perturbative expansion in real 



time. The lowest order renormalization group equation ( 4.45 ) resums the leading logarithms, while 



the lowest order dynamical renormalization group equation resums the leading secular terms. 
Dynamical renormalization group: fixed points and coarse-graining 

The similarity between the renormalization of distribution functions and the renormalization of 
couplings suggests that the collision terms of the quantum kinetic equation can be interpreted as beta 
functions of the dynamical renormalization group and that the space of distribution functions can 
be interpreted as the space of coupling constants. The dynamical renormalization group trajectories 
determine the fiow in this space, therefore fixed points of the dynamical renormalization group 
describe stationary solutions with given distribution functions. Thermal equilibrium distributions are 
thus fixed points of the dynamical renormalization group. Furthermore, there can be other stationary 
solutions with non-thermal distribution functions, e.g., those describing turbulent behavior |p6| . 
Linearizing the renormalization group equation around the fixed points corresponds to linearizing 
the kinetic equation near equilibrium and the linear eigenvalues are related to the relaxation rates 
in the relaxation time approximation. 

One can establish a closer relationship to the usual renormalization program of field theory in its 
momentum shell version with the following alternative interpretation of the secular terms and their 
resummation. 

The initial distribution at a time to is evolved in time perturbatively up to a time to + At 
such that the perturbative expansion is still valid, i.e., ^kin S> Ai with tkin the kinetic time scale. 
Secular terms begin to dominate the perturbative expansion on a time scale At ^ tniic with t„iic the 
microscopic time scale. Thus if there is a separation of time scales such that tkin ^At^ imic, then 
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in this intermediate asymptotic regime perturbation theory is rehable but secular terms appear and 
can be isolated. 

A renormalization of the distribution function absorbs the contribution from the secular terms. 
The "renormalized" distribution function is used as an initial condition at to + At to iterate forward 
in time to to + 2At using the perturbative expansion, but with the propagators in terms of the 
distribution function at the time to + At. This procedure can be carried out "infinitesimally" (in the 
sense compared with the kinetic time scale) and the differential equation that describes the changes 
of the distribution function under the intermediate asymptotic time evolution is the dynamical 
renormalization group equation. 

This has an obvious similarity to the renormalization group in terms of integrating in momentum 
shells. The result of integrating out degrees of freedom in a momentum shell are absorbed into a 
renormalization of the couplings and an effective theory at a lower scale but in terms of the effective 
couplings. This procedure is carried out infinitesimally and the differential equation that describes 
the changes of the couplings under integrating out degrees of freedom in these momentum shells is 
the renormalization group equation. 

An important aspect of this procedure of evolving in time and "resetting" the distribution func- 
tions is that in this process it is implicitly assumed that the density matrix at later times remains 
diagonal in the basis of free quasiparticle number operators. Clearly, if at the initial time the den- 
sity matrix was diagonal in this basis, because the interaction Hamiltonian does not commute with 
the density matrix, off-diagonal density matrix elements will be generated upon time evolution. In 
resetting the distribution functions and using the propagators in terms of these updated distribution 
functions we have neglected off-diagonal correlations, e.g., correlations of the form (o(k, t)a(k, f)) 
and its complex conjugate will be generated upon time evolution. In neglecting these off-diagonal 
terms we are introducing coarse-graining. Indeed, two stages of coarse-graining had been introduced: 
(i) integrating in time up to an intermediate asymptotics and resumming the secular terms neglect 
transient phenomena, i.e., averages over the microscopic time scales, and (ii) off-diagonal matrix ele- 
ments (in the basis of free quasiparticle number operators) had been neglected. The coarse graining 
also has an counterpart in the language of the renormalization group in Euclidean field theory, i.e., 
neglecting the irrelevant couplings that are generated upon integrating out shorter scales. Keeping 
all of the correlations in the density matrix would be equivalent to Wilson's approach to renormal- 
ization group, in which all possible couplings are included in the Lagrangian and all of them are 
maintained in the renormalization on the same footing. 
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4.4 Application: 0(4) Linear Sigma Model 

In this section we apply the dynamical renormalization group method to derive kinetic equations for 
pions and sigma mesons in the 0(4) linear sigma model, which is an effective theory of low energy 
QCD describing the hadronic world. We recall that the sigma condensate (cr) represents the quark 
condensate (qq) in the sense that they both have the same transformation properties, corresponding 
to two flavors of massless quarks. At high temperatures, quarks and gluons exist in a deconfined, 
chirally symmetric phase. At some critical temperature of order 160 MeV a transition to a hadronic 
phase occurs. In this confined and chiral symmetry broken phase the quark condensate is nonzero. 

The Lagrangian density of the 0(4) linear sigma model in the limit of exact chiral symmetry, 
i.e., without an explicit symmetry-breaking term, is given by 

£[7^, a] = i {d,nf + 1 {d,af ^^{n' + a'- flf , (4.48) 

where tt — (tt^, tt^, tt'^) and /jr ^ 93 MeV is the pion decay constant. At high temperature T > Tc, 
where Tc /jr is the critical temperature, the 0(4) symmetry is restored by a second order phase 
transition [ p7[ . 

In the symmetric phase pions and sigma mesons are degenerate, and the linear sigma model 



reduces to a self-interacting scalar theory, analogous to that discussed in Sec. 4.2. Thus, we limit our 
discussion here to the low temperature broken symmetry phase in which the temperature T <C /tt- 
Since at low temperature the 0(4) symmetry is spontaneously broken via the sigma condensate 
(cr), we shift the sigma field cr(x, t) — (t(x, t) + v, where v is temperature dependent and yet to be 
determined. 

In thermal equilibrium v is fixed by requiring that ((t(x, t)) = to all orders in perturbation 
theory for temperature T < Tc- In nonequilibrium this split must be performed on both forward 
(-I-) and backward (— ) time branches, writing the sigma field cr*(x, i) as 

cr±(x,t) =a±(x,i)-f (a±(x,i)) 0, (4.49) 

where ( • ) denotes the nonequilibrium expectation value. The expectation value v is determined by 
requiring that the expectation value of (t(x, t) vanishes to all orders in perturbation theory via the 
tadpole method p5| . One finds to one- loop order the equation that determines v is given by 

v[v^~f^ + {7v^)+3{a')]=0. (4.50) 

We emphasize that in perturbation theory v depends on the initial nonequilibrium distribution 
functions of the pions and sigma mesons, which in turn implies that after dynamical renormaliza- 
tion group resummation v acquires implicit time dependence through the time dependence of the 



distribution functions (see Sec. 3.2). Once the solution of this equation for v is used in the per- 
turbative expansion to one loop order, the tadpole diagrams that arise from the shift in the field 
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cancel automatically to this order. This feature of cancelation of tadpole diagrams that would result 
in an expectation value of a by the consistent use of the tadpole equation persists to all orders in 
perturbation theory. 



A solution of Eq. (4.50) with v ^ signals broken symmetry and massless pions (in the strict 
chiral limit). Therefore once the correct expectation value v is used, the one-particle-reducible (IPR) 
tadpole diagrams do not contribute in the pcrturbative expansion of the kinetic equation. To this 
order the inverse pion propagator for zero energy and momentum reads 

G-\uj = 0, k = 0) = -A y ~fl + (tt^) + 3(^2)] , (4.51) 

which vanishes whenever u 7^ by the tadpole condition Eq. ( 4.50| ), hence the Goldstone theorem is 
satisfied and the pions are the Goldstone bosons. 

The study of the relaxation of sigma mesons (resonance) and pions near and below the chiral 
phase transition is an important phenomenological aspect of low energy chiral phenomenology with 
relevance to heavy ion collisions. Furthermore, recent studies have revealed interesting features 
associated with the dropping of the sigma mass near the chiral transition and the enhancement 
of threshold effects with potential experimental consequences The kinetic approach described 
here could prove useful to further assess the contributions to the width of the sigma meson near the 
chiral phase transition, this is an important study on its own right and we expect to report on these 
issues in the near future. 

With the purpose of comparing to recent results, we now focus on the situation at low tempera- 
tures under the assumption that the distribution functions of sigma mesons and pions are not too far 
from equilibrium, i.e., cool pions and sigma mesons. At low temperatures the relaxation of pions and 
sigma mesons will be dominated by the one-loop contributions, and the scattering contributions will 
be subleading. The scattering contributions are of the same form as those discussed for the scalar 
theory and involve at least two distribution functions and are subdominant in the low temperature 
limit as compared to the one-loop contributions described below. 

Since the linear sigma model is renormalizable and we focus on the medium (finite-temperature) 
effects, we ignore the vacuum (zero-temperature) ultraviolet divergences which can be absorbed into 
a renormalization of f^^. For a small departure from thermal equilibrium, one can approximate (tt^) 
and ((7^) by their thermal equilibrium values 

4 ' ^ ' J (2^)3c.q 



= "M-q), (4.52) 



respectively, where tjq ~ ■y/q^ + m^ with the sigma mass = 2\v^ being determined self- 
consistently. In the low temperature limit T <C /tt, one finds v"^ — f^[\ — 0{T'^ / f^)] and TOq- = 
V2A/7r[l — 0(T^//^)]. Thus in our discussion below where T ^ /tt, we set v — f-^ and rria = \/2A/7r. 
The main reason behind this analysis is to display the microscopic time scales for the mesons: 
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^mic — l/"^cr and ij^jjp = 1/fc with k being the momentum of the pion. The validity of a kinetic 
description wiU hinge upon the kinetic time scales being much longer than these microscopic scales. 

Upon the above replacement, one can rewrite the Lagrangian for the cool linear sigma model to 
lowest order as 

^[tt, <y]^l {d^^f + I {d^<jf - - A/, {an' + - ^ + '^')' , (4-53) 

where the bar over the shifted sigma field has been omitted for simplicity of notation. 

Our goal in this section is to derive the kinetic equations describing pion and sigma meson 
relaxation to lowest order. The unbroken 0(3) isospin symmetry ensures that all the pions have 
the same relaxation rate, and sigma meson relaxation rate is proportional to the number of pion 
species. Hence for notational simplicity the pion index will be suppressed. We now begin to study 
the quantum kinetics of cool pions and sigma mesons. 

4.4.1 Quantum kinetics of cool pions and sigma mesons 

Without loss of generality, in what follows we discuss the relaxation for one isospin component, say 
TT^ , but we suppress the indices for simplicity of notation. As before, we consider the case in which 
at an initial time t — to, the density matrix is diagonal in the basis of free quasiparticles number 
operators, but with nonequilibrium initial distribution functions n^{to) and n^{to). 

Pions 

The pion field operator and the corresponding conjugate momentum in the Heisenberg picture can 
be written as (in momentum space) 

7r(k, = [a-(k, t) + at (-k, t)] , n,(k, t) = -^■y/| [a,(k, t) - flt (_k, t)] , (4.54) 



The expectation value of pion number operator can be expressed in terms of 7r(k, t) and n7r(k, t) as 

nlit) = {ai{k,t)a^{k,t)) 

= ^ [(n,(k, t)n^i-k, t)) + fc2 (^(k, i)7r(-k, t))] ~ i. (4.55) 

Use the Heisenberg equations of motion, to leading order in A, we obtain (no tadpole diagrams are 
included since these are canceled by the choice of v) 

|<(.) = -^lhn^|, / (^(.-(k-q,t).+ (q,0.-(-k,O). (4.56) 

The expectation values can be calculated perturbatively in terms of nonequilibrium vertices and 
Green's functions. To 0{X) the right-hand side of Eq. ( 4.56| ) vanishes identically. Fig. 4.2a shows 
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t I \ *' + t ' f 

+ +t" + -f" 



(a) 



(b) 



Figure 4.2: (a) The Feynman diagrams that contribute to the quantum kinetic equation for the pion 
distribution function at one foop order. The sohd Hne is the sigma meson propagator and the dashed 
hue is the pion propagator, (b) The only contribution on-shell is the decay of a sigma meson into 
two pious minus the reverse process. 

the Feynman diagrams that contribute to order A^. After some algebra, one obtains 

-r.<{t) = / 7^ / dt" M(io)cos[Awi(i-t")]+AA2(to)cos[Au;2(i-t")] 

at k J [2'Ky qivp Jt^ L 

+ J\f3{ta)cos[Auj3{t-t")]+J\f4{to)cos[Auj4{t-t")] , (4.57) 

where p = k + q, 

Acji = k + q + ujp, Auj2 — k ~ q~ ujp, 

Auj^ — k — q + Up, Alu4 = k + q — ujp, (4.58) 

m) = [l + nimi + nimi+n;{t)]^nl{t)nl{t)nl{t), 
= [l + nl{t)]nl{e)nl{t)^nl{t)[l + nl{tW+nl{t% 

M^it) = [l + nUtm + nl{t)]n^pit)-nl{t)nl{t)[l+n^p{t)]. (4.59) 

The different contributions have a very natural interpretation in terms of "gain minus loss" processes. 
The first term in brackets corresponds to the process — + ct + tt + tt minus the process ct + tt + tt —^ 0, 
the second and third terms correspond to the scattering tt + cr — > tt minus tt — > tt + cr, and the last 
term corresponds to the decay of the sigma meson cr ^ tt + tt minus the inverse process vr + tt — s- cr. 

Just as in the scalar case, since the propagators entering in the perturbative expansion of the 
kinetic equation are in terms of the distribution functions at the initial time, the time integration 



and 
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can be done straightforwardly leading to the following equation: 

^nlit) = /^°"rfc.7^.[c.,k;AG(to)] ^^(^ ' fe)^ " ^°)] , (4.6O) 
dt tt{uj - k) 

where 7?,7r[w, k; A/j(io)] is given by 

P f (Pa 1 r 

7^,[^,k;AG■(^o)] = -f irrk (5(A^i)M(to) + <5(Ac^2)AA2(io) 

+ 5(Acj3)AA3(io) + <^(Acj4)AA4(to)l. (4.61) 



can be solved by direct integration over t with the given initial condition at t^, leading to 



nl{t)=nl{t,) + \^ ^°°d^7^.[c.,k;AG-(to)] ^ ^^2^ ■ (4.62) 

J -00 nluj-ky 



Potential secular term arises at large times when the resonant denominator in ( [4.62[ ) vanishes, i.e. 



uj ^ k. A detailed analysis reveals that 7?.7r[w, k; A/}(to)] is regular at lu = k, hence upon using 



( 4.24 ) one finds that at intermediate asymptotic time k{t — to) ^ 1, the time evolution of the pion 



distribution function reads 

nJJ(i) = nl.{to) + \^ TZjrlk, k; A/'j(io)] {t — to) + nonsecular terms, (4.63) 

where 7?.7r[A:, k; A/}(to)] does not depend on to explicitly. 

At this point we would be tempted to follow the same steps as in the scalar case and introduce the 
dynamical renormalization of the pion distribution function. However, much in the same manner 
as the renormalization program in a theory with several coupling constants, in the case under 
consideration the tt field and the a field are coupled. Therefore one must renormalize all of the 
distribution functions on the same footing. Hence our next task is to obtain the kinetic equations 
for the sigma meson distribution functions. 

Sigma mesons 

The field operator for the sigma meson and the corresponding conjugate momentum in the Heisen- 
berg picture can be written as (in momentmn space) 

a(k, t) = [a,(k, t) + ali-k, t)] , H,(k, t) = -i^ [a,(k, t) ~ at (-k, t)\ , (4.64) 



with ~ ^yW + rn^. Again, for notational simplicity we suppress the pion isospin index. The 
expectation value of sigma meson number operator can be expressed in terms of cr(k, t) and Ho.(k, t) 
as 

<(t) = {ai{k,t)a„{k,t)) 

= ^[{IlAk,t)n,i-k,t))+ujl{a{k,t)a{-k,t))]-^. (4.65) 
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Using the Heisenberg equations of motion to leading order in A, and requiring again that the tadpole 
diagrams are canceled by the proper choice of we obtain 



d ai.\ 3A/^ d 

—n-i.(t) = lim — — 

dt ' ujk t'~.t dt' 



d\ 



(27r)3/2 

+ 3(a+(k-q,t)a+(q,i)a-(-k,0) 



(7r+(k-q,i)7r+(q,<)a-(-k,t')) 



(4.66) 



where the factor three for the first term accounts for three isospin components of the pion field. 



To 0{\) the right hand side of Eq. ( 4.66 ) vanishes identically. Fig. Oa depicts the one-loop 
Feynman diagrams that enter in the kinetic equation for the sigma meson to order A^. To the 
same order there will be the same type of two loops diagrams as in the self-interacting scalar theory 
studied in the previous section, but in the low temperature limit the two-loop diagrams will be 
suppressed with respect to the one-loop diagrams. Furthermore, in the low temperature limit, the 
focus of our attention here, only the pion loops will be important in the relaxation of the sigma 
mesons. A straightforward calculation leads to the following expression 



dt 



Kit) 



d^q 
(2^ 



dt" 



— (M^ito)cos[ALul{t-t")] 
pq \ 



+ Af^itQ) cos[AwJ(i - t")] + Af^ito) cos[AwJ(i - t")] 
-A(r(io)cos[A<(i-i")]^ ^ 



(4.67) 



where p = k -t- q. 



Aoj'. 



7r((T) 



t^k + g i^q) + P (Wp), = tJk + g (t^q) - P (Wp), 



- <? (t^q) +P (Wp), 



(4.68) 



and 



<^\t) = [l + <(t)][l + <('^)(t)][l + n^(-)(i)] -n£(t)<('^)(t)n^('^)(i), 
[1 + Kit)] [1 + {t)]n<"^ (t) - Kit) n^(-) (t) [1 + n^(-) (t)], 
[1 + Kit)] n<'^\t) [1 + n<'^\t)] <(i) [1 + <-)(t)] n^M(i), 
Mr>{t) = [1 + K{t)] n<'^\t) nl^^\t) - Kit) [1 + <('^)(t)][l + <-)(t)]. 



Af;^^\t) 
U;^^\t) 



(4.69) 



Although the above expression is somewhat unwieldy, the different contributions have a very 
natural interpretation in terms of the "gain minus loss" processes. In the first parentheses (the pion 
contribution) the first term corresponds to the process (T + tt+tt minus the process (t+tt + tt 0, 
the second and third terms correspond to the scattering tt — > tt + cr minus tt + ct — s- tt, and the last 
term corresponds to the decay of the sigma meson ct ^ tt + tt minus the inverse process tt 4- vr — > cr. 
Likewise, in the second parentheses (the sigma meson contribution) the first term corresponds to 
the process — > c + cr + ct minus the process cr + a + o" — > 0, the second and third terms correspond 
to annihilation of two sigma mesons and creation of one sigma meson minus the inverse process, and 
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Figure 4.3: (a) The Feynman diagrams that contribute to the quantum kinetic equation for the 
sigma meson distribution function at one-loop order. The sohd hue is the sigma meson propagator 
and the dashed hne is the pion propagator, (b) The only contribution on-shell is recombination of 
two pions into a sigma meson minus the decay of a sigma meson into two pions. 



the last term corresponds to annihilation of a sigma meson and creation of two sigma mesons minus 
the inverse process. 

Since the propagators entering in the perturbative expansion of the kinetic equation are in terms 
of the distribution functions at the initial time, the time integration can be done straightforwardly 
leading to the following equation: 



sin[{uj - ^k)ft - to)] 
7r(w - Wk) 



(4.70) 



where 



pq 



9 



(4.71) 



Just as before 7?,cr[w, k; A/}(to)] is fixed at initial time ^o, Eq. (4.7C) can be integrated over t with 
the given initial condition at Iq, thus leading to 



Kit) = <(to) + 



dLjn^[LO,k;J\f,{to)]- 



(4.72) 



7r(cj-tJk)^ 

At intermediate asymptotic times m^{t — to) S> 1, potential secular term arises when lu fa uj^ 



in Eq. ( 4.72 ). We notice that, although TZa[co,'k;Afj{tQ)] has (infrared) threshold singularities at 
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Lu — ±fc, it is regular on the sigma meson mass shell. This observation will allow us to explore a 
crossover behavior for very large momentum later. 

Since TZa-i'-^^ k; A/'j(to)] is regular near the resonance region lo = iciJk, the behavior at intermediate 
asymptotic times is given by 

'^k(^) = n-k(^o) + A^7^CT[t^k,k;7VJ(^o)](^ - ^o) + nonsecular terms. (4.73) 

We note that the perturbative expansions for the pion and sigma meson distribution functions 
contain secular terms that grow linearly in time, unless the system is initially prepared in thermal 
equilibrium. To resum the secular terms, we must now renormalize both Eqs. ( 4.63| ) and (4.73) 



simultaneously, because, as discussed in Sec. O, this situation corresponds to that with two relevant 
couplings in Euclidean field theory. The renormalization is achieved by introducing the renormalized 
initial distribution functions 'r^}'^^ (t) for pions (sigma mesons) that are related to the bare initial 
distribution functions '^^^'^^(^o) via the renormalization constants Z^^'^\T,to) 

nl^'^^ {to) = (r, to) nl^"^ (r), Z^^'^' (r, to) = 1 + z^^'^^ (r, to) + ©(A^), (4.74) 

where r is an arbitrary renormalization scale at which the secular terms will be canceled. As in the 
case of the scalar theory, the renormalization constants z^^'^''(r, to) are chosen so as to cancel the 
secular terms at the arbitrary scale r. Upon substitute Eq. ( |4.74| ) into Eq. (4.62) and choosing 



^k' ^(r,to) = -(r-to)^i^^^^(^--^, (4.75) 

"k iv 

consistently up to 0{\^), one obtains 



nJ-^Ht) = <(^)(r) + (t-T)7^^(^)[fc(c.k),k;AG■(T)] +0(A4). (4.76) 

The independence of n^^"^^ (t) on the arbitrary renormalization scale r leads to the simultaneous set 
of dynamical renormalization group equations to lowest order: 

|-<(r) = A2 7^,[fc,k;AG■(T)], 

-^<(t) = A2 7^,K,k;AG■(T)]. (4.77) 

These equations have an obvious resemblance to a set of usual renormalization group equations for 
"couplings" and n^, where the right-hand sides are the corresponding "beta functions". 

As before choosing the arbitrary scale r to coincide with the time t and keeping only the terms 
whose delta functions have support on the mass shells, we obtain the quantum kinetic equations for 
pions and sigma mesons: 

d _ 7rA2/2 f d^q S{k + q-ujp) 



-<(t),^^(t)[l + n;(t)]}, (4.78) 



62 



-<(i)[l + n^(i)][l + n^(t)]}. (4.79) 



The processes that contribute to Eq. (4.78) are depicted in Fig. 2b and those that contribute to 



Eq. ( |4.79| ) are depicted in Fig. 3b. 



4.4.2 Relaxation rate for pions and sigma mesons 

Thermal equilibrium is a fixed point of the dynamical renormalization group equations, i.e., a sta- 



tionary solution of the quantum kinetic equations ( 4.7§| ) and (4.79). Near equilibrium a Linearized 
kinetic equation for pions (sigma mesons) can be obtained in the relaxation time approximation, in 
which only distribution function for pions (sigma mesons) of momentum k are slightly perturbed off 
equilibrium whereas all the other modes are in equilibrium. In the relaxation time approximation 
for pions (sigma mesons), one obtains 

|K^'^^W--r.(.)(fc)K^^Ht), (4.80) 

where (fc) is the pion (sigma meson) relaxation rate, which is identified with twice the damping 
rate of the corresponding mean field. 



Linearizing Eq. (4.7S) in the relaxation time approximation, we obtain the pion relaxation rate 
to be given by 

ttX^JI f d^q neiq) -nsiujp) 



A2/2r I - e-/3(™'/4/c+fe) 
" "4^ 1 - e-/3™J/4fe ■ (^-^^^ 

This is a remarkable expression because it reveals that the physical processes that contribute to pion 
relaxation are the decay of a sigma meson to two pions and its inverse process, i.e., a ^ tt + tt. 
The sigma mesons present in the medium can decay into pions and this increases the number of 
pions, but at the same time pions recombine into sigma mesons, and since there are more pions 
in the medium because they are lighter the loss part of the process prevails producing a non-zero 
relaxation rate for the pion distribution function. This is an induced phenomenon in the medium in 
the very definitive sense that the decay of the heavier sigma meson induces the relaxation of the pion 
distribution function, it is a collisionless process. We note that such induced relaxation of pions is 
analogous to the relaxation of fermions in a fcrmion-scalar plasma discusses in Chap. |[ the latter 
is induced by the decay of a massive scalar into light fcrmion pairs. 
For soft, cool pions (fc ^ T <C /t), the relaxation rate reads 



63 



The exponential suppression in the above expression is a consequence of the heavy sigma mass. Our 
results of the pion relaxation rate are in agreement with the pion damping rate found in Rcf. p8| . 
These results (accounting for the factor two necessary to relate the relaxation rate to the damping 
rate) also agree with those reported recently in Rcf. |Q , wherein a related and clear analysis of pion 
and sigma meson damping rates was presented. 

Likewise, the sigma meson relaxation rate in the relaxation time approximation is found to be 
given by 

d^q [l + riBiq) +nB{p)] 



mi 



(27r)3 
2T 



1 



, In- 

k 1 



pq 

g-/3(i^k+fe)/2 



(4.83) 



g-/3(L^k-fe)/2 

We note that the first temperature- independent term in T^(k) is the usual zero-temperature sigma 
meson decay rate | 100 |, whereas the finite-temperature terms result from the same processes a ^ 
TT + TT that also determine TTr{k). For soft, cool sigma meson (fc <C T <C /tt), one obtains 



^cothf!^) 



(4.84) 



which agrees with the decay rate for a sigma meson at rest 101, 102]. 

On the other hand, let us consider the theoretical high temperature and large momentum limit 
k 3> m„ > T and ujk — k <ti T. In this limit Ta{k) becomes logarithmic (infrared) divergent. The 
reason for this divergence is that, as was mentioned above, TZa[uj^ k; A/}(to)] has an infrared threshold 
singularity aX ld = k arising from the terms proportional to {to), which accounts for the emission 
and absorption of collinear massless pions. In the presence of this threshold singularity, we can 



no longer apply Fermi's golden rule (4.24). Instead, we must analyze the long-time behavior of 



Eq. (4.67) more carefully. 

Understanding the influence of threshold behavior of the sigma meson on its relaxation could 
be important in view of the recent proposal by Hatsuda and collaborators p9[ that near the chiral 
phase transition the mass of the sigma meson drops and threshold effects become enhanced with 
distinct phenomenological consequences. 



4.4.3 Threshold singularities and crossover in relaxation 

In order to understand how the (infrared) threshold divergence modifies the long-time behavior, 
let us focus on the sigma mesons with large momentum k 3> > T. This situation is not 
relevant to the phenomenology of the cool pion-sigma meson system for which relevant temperatures 
are T <^ m„. However studying this theoretical limiting case will reveal important insight on 
how threshold divergences invalidate the usual Fermi's golden rule analysis, which leads to energy 
conservation delta functions in the intermediate asymptotic regimes. This issue will become physical 
relevant in the case of gauge field theories studied below (see Chap. ||). 
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To present this case in the simplest and clearest manner, we will proceed in the relaxation time 
approximation. Keeping only the infrared divergent contribution in TZa-[ijj,'k;Mj{tQ)], i.e., the term 



proportional to Af^{tf)), one finds that Eq. (4.75) simplifies to 

cos[(ti; — LUk)(t 



+00 

dujTZ{uj,k) — 



to)] 



where 



7^(w, k) 



7r(cj - ujk)^ 
37rA2/2 f d^q l + nB{q)+nBip) 



(4.85) 



2wk 

87r(jJk 



(27r)3 pq 

2T 1- e-'3('^+fc)/2 



■S{UJ 



(4.86) 



which has a infrared threshold singularity at = fc. At intermediate asymptotic times 171^(1— to) ^ f , 



the integral over lu in Eq.(4.85) is dominated by the region near lo k. In the limit k 2> > T, 
by approximating 



7^(w,k) 



47rfc2 



In 



2T 



UJ - 



+ 0{u-k), 



the Ld integral (denoted as X) for ma{t — to) ^ I can be evaluated in a closed form 



47rfc2 



(4.87) 



(4.88) 



where T = T/{uj]^ - fc) « 2kT/ml and t = {uJk - k){t - to) 
variables and 

T(x) = X (^In 2T + ci{x) - ^ 
with ci(x) being the cosine integral function 



mj(t — to) /2k are dimensionless 



(4.89) 



ci(a;) = — 



+00 



dt 



cost 



(4.90) 



The function J-{x) has the following asymptotic behaviors 
J^ix) = 



(4.91) 



x [ln2a;T + 7- 1 + (a;2)] for a; < 1, 
X [In 22^ + (.T-2)] for a; > 1, 

where 7 = 0.577 ... is the Euler-Mascheroni constant. Thus, for fixed and large k one can see that 
there is a crossover time scale tc ~ 2fc/m^ on which the time dependence of J-^(t) changes from the 
tint dependence for t ~ tg <ti t^ to the linear dependence for t — t^ 3> tc- Consequently, in the large 
momentum limit as the sigma meson mass shell approaches the threshold, this crossover time scale 
becomes very long such that an "anomalous" (nonlinear) secular term of the form tint dominates 
during most of the time whereas the usual linear secular term ensues at very large times. 

The secular terms can be resumed by implementing the dynamical renormalization group method 
[see Eq. ( 4.73| )] with the following choice of renormalization constant (in terms of dimensionless time 
variable) 



z^irM) 



47rfc2 



T{T-to). 



(4.92) 
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The resultant dynamical renormalization group equation reads 



J o \2 J•2^^ J 



whose solution is given by 



Sn^{t) — 6n^{tQ) exp 



(4.94) 



In the large momentum limit, using Eq. (4.91) we find that the crossover in the form of the secular 
terms results in a crossover in the sigma meson relaxation; an "anomalous" (nonexponcntial) re- 
laxation will dominate during most of the time and the usual exponential relaxation ensues at very 
large times. 

This above discussion has revealed several important features highlighted by a consistent resum- 
mation via the dynamical renormalization group: 

(i) Threshold infrared divergences result in a breakdown of Fermi's golden rule. The secular 
terms of the perturbative expansion are no longer linear in time but include logarithmic contributions 
arising from these infrared divergences. 

(ii) The concept of the damping rate is directly tied to exponential relaxation. The infrared 
divergences of the damping rate reflect the breakdown of Fermi's golden rule and signal a very 
different relaxation from a simple exponential. 

(iii) Whereas the usual calculation of damping rates will lead to a divergent result arising from the 
infrared threshold divergences, the dynamical renormalization group approach recognizes that these 
threshold divergences result in secular terms that are non-linear in time as discussed above. While in 
relaxation time approximation linear secular terms lead to exponential relaxation and therefore to an 
unambiguous definition of the damping rate, non-linear secular terms lead to novel nonexponcntial 
relaxation phenomena for which the concept of a damping rate may not be appropriate. 

This discussion of threshold singularities and anomalous relaxation has paved the way to studying 
the case of gauge field theories in next chapter, wherein the emission and absorption of magnetic 
photons that are only dynamically screened lead to a similar anomalous relaxation |84]. 



4.5 Pinch Singularities and their Resolution 

An important difference between the approach to nonequilibrium evolution described by quantum 
kinetic equations advocated in this work and that often presented in the literature is that we work 
directly in real time not taking Fourier transforms in time. This has to be contrasted with the 
so-called real-time formalism (RTF) of (non)equilibrium quantum field theory, in which there are 
also a closed-time-path contour and four propagators but the propagators and quantities computed 
therefrom are all in terms of temporal Fourier transforms. 
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In thermal equilibrium the Fourier representation of these four propagators for a scalar field are 
given by 

1 



le 



2mnB{\ko\)S{K' -mis), 



G^+{K) = 27n [9{ko) + nadfeol)] SiK' - mi^), 
G+-{K) = 2m [e{-k„) + nB{\ko\)] S{K^ - mlg), 



(4.95) 



where K — (fco,k;) is the four- momentum with = k^ ^ k'^, whereas out of equilibrium the 
distribution functions are simply replaced by non-thermal ones, i.e., nsdAol) '^k;(io)- Using the 
integral representation of the Heaviside step function 



(4.96) 



one can easily show that the propagators given in Eq. ( 4.95| ) and the ones obtained by replacing the 
thermal equilibrium distributions by the nonequilibrium ones are, respectively, the temporal Fourier 



transforms of those given in Eqs. (3.5) and (4.2). The Fourier transforms of the free retarded and 
advanced propagators are obtained similarly and read 

1 



G^'/\k) 



(4.97) 



K'^ -mlff±isgn{ko)e' 

Several authors have pointed out that the calculations using the CTP formalism in terms of 



the standard form of free propagators in Eq. (4.95) or those obtained by the replacement of the 
distribution functions by the nonequilibrium ones, lead to pinch singularities |87,88, 103 107]. 

In a consistent perturbative expansion both the retarded and advanced propagators contribute 
and pinch singularities arise from the product of these, for example for a scalar field this product is 
of the form 



G^iK)G^{K) = 



[K^ -mls + isgn{ko)e][K^ 



off 



isgn(A:o)e] 



(4.98) 



For finite e this expression is regular, whereas when e 0+ it gives rise to singular products such 
as [5{K^ — ™cff)]^ ^ discussed in Refs. |^,^, 103-107|. Singularities of this type are ubiquitous in 
nonequilibrium and are not particular to scalar field theories. 

A detailed analysis of these pinch terms reveals that they do not cancel each other in perturbation 

10^-[l05|]. Indeed, this severe problem 



theory unless the system is in thermal equilibrium |87 
has cast doubt on the validity or usefulness of the CTP formulation to describe nonequilibrium 
phenomena |^8|] . Although this singularities have been found in many circumstances and analyzed 
and discussed in the literature often, a systematic and satisfactory treatment of these singularities is 
still lacking. In Ref. |107| it was suggested that including an in-medium width of the quasiparticles 
to replace Feynman's "ie-prescription" does provide a physically reasonable solution, however this 
clearly casts doubt on the consistency of any perturbative approach to describe even weakly out of 
equilibrium phenomena. 
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Recently some authors have conjectured that pinch singularities in perturbation theory might 



be attributed to a misuse of Fourier transforms (for a detailed discussion see Refs. [103-105|). As 
an illustrative and simple example of these type of pinch singularities, these authors discussed the 
elementary derivation of Fermi's golden rule in time-dependent perturbation theory in quantum 
mechanics. In calculating total transition probabilities there appears the square of energy conser- 
vation delta function, which arises due to taking the infinite time limit of scattering probabilities. 
In this setting, such terms are interpreted as the elapsed scattering time multiplied by the energy 
conservation constraint rather than a pathological singularity. 



A close inspection of Eq. ( 4.98|) reveals that the pinch term is the square of on-shell condition 
for the free quasiparticle, which implies a temporal Fourier transform in the infinite time limit and 
of the same form as the square of the energy conservation constraint for the transition probability 
obtained in time-dependent perturbation theory. 



By assuming that the interaction duration time is large but finite, Niegawa |104| and Greiner 



and Leupold |105] showed that for a self- interacting scalar field the pinch part of the distribution 



function can be regularized by the interaction duration time as |104, 105 1 



pinch 



(t)^(t-to)rr(io), 



(4.99) 



where "~" denotes that only the pinch singularity contribution is included, t — to is the interaction 
duration time, and rj^°*(io) is the net gain rate of the quasiparticle distribution function per unit 
time 



K'Hto) = TT- [l + nk(to)]S<(c^k,k)-nk(io)S>(c^k,k) 



(4.100) 



In the above expression, the on-shell self-energies I]^(wk,k) are calculated in terms of the initial 



distribution functions nk(to) 1 104, 105 1. 

Upon comparing Eqs. ( 4.9E ) and ( 4.100| ) with Eqs. ( 4.25 ) and ( 4.28 ), respectively, we clearly see 
the equivalence between the linear secular terms in the perturbative expansion and the presence of 



pinch singularities in the usual CTP formalism. Furthermore, in the discussion following Eq.(4.25) 
we have recognized that secular terms are not present if the system is in equilibrium, much in 
the same manner as the case of pinch singularities as discussed originally by Altherr Thus 
our conclusion is that pinch singularities are a temporal Fourier transform representation of linear 
secular terms. 

The dynamical renormalization group provides a systematic resummation of these secular terms 
and a consistent formulation to implement the renormalization of the distribution functions suggested 



by Niegawa. [104|. We emphasize that the dynamical renormalization group approach explains the 
physical origin of the pinch singularities in terms of secular terms and Fermi's golden rule, and 
provides a consistent and systematic resummation of these secular terms that lead to the quantum 
kinetic equation as a renormalization group equation that determines the time evolution of the 
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distribution function. This result justifies in a systematic manner the conclusions and interpretation 



obtained in Ref. |107] where a possible regularization of the pinch singularities was achieved by 
including the width of the quasiparticle obtained via the resummation of hard thermal loops. 

Furthermore, we stress that the dynamical renormalization group is far more general in that it 
allows to treat situations where the long time evolution is modified by threshold (infrared) singulari- 
ties in spectral densities, thereby providing a resolution of infrared singularities in damping rates and 
a consistent resummation scheme to extract the asymptotic time evolution of the distribution func- 
tion. The infrared singularities in these damping rates is a reflection of anomalous (nonexponential) 
relaxation as a result of threshold effects. 

The pinch singularities signal the breakdown of perturbation theory, just as the secular terms in 
real time, however, the advantage of working directly in real time is that the time scale on which 
perturbation theory breaks down is recognized clearly from the real-time perturbative expansion 
and is identified directly with the kinetic time scale. The dynamical renormalization group justifies 
this identification by providing a resummation of the perturbative series that improves the solution 
beyond the intermediate asymptotics. 

The resolution of pinch singularities via the dynamical renormalization group is general. As 
originally pointed out in Ref. pq] the pinch singularities typically multiply expressions of the form 



Eq. (4.100) which vanish in equilibrium, just as the linear secular terms multiply similar terms in 



the real-time perturbative expansion, as highlighted by Eq. (4.28). These terms are of the typical 
form gain minus loss, in equilibrium they vanish, but their nonvanishing simply indicates that the 
distribution functions are evolving in time and it is precisely this time evolution that is described 
consistently by the dynamical renormalization group. 



4.6 Conclusions 

The goal of this chapter is to provide a novel method for obtaining quantum kinetic equations 
from a field theoretical and diagrammatic perturbative expansion that is improved via a dynamical 
renormalization group resummation directly in real time. 

The first step of this method is to use the microscopic equations of motion to obtain the evolution 
equation of the quasiparticle distribution function, i.e., the expectation value of the quasiparticle 
number operator in the initial density matrix. This evolution equation can be solved in a consistent 
diagrammatic perturbative expansion and one finds that the solution for the time evolution of the 
distribution function features secular terms, i.e., terms that grow in time. In perturbation theory 
the microscopic and kinetic time scales are widely separated, hence there is a regime of intermediate 
asymptotics in time, within which (i) perturbation theory is valid and (ii) the secular terms dominate 
the time evolution of the distribution function. 
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The second step is to implement the dynamical renormalization group to rcsum the secular 
terms. A renormalization of the distribution function absorbs the contribution from the secular 
terms on a given renormalization time scale, thereby improving the perturbative expansion. The 
arbitrariness of this renormalization scale leads to the dynamical renormalization group equation, 
which is interpreted as the quantum kinetic equation. In relaxation time approximation linear one 
recognizes that secular terms correspond to usual exponential relaxation, whereas nonlinear secular 
terms correspond to anomalous (nonexponential) relaxation. 

We first test this new dynamical renormalization group method within the familiar self-interacting 
scalar field theory, not only reproducing but also generalizing the previous results in the literature. 
We then move on to apply the this method to study quantum kinetics of a gas of cool pions and 
sigma mesons described by the 0(4) linear sigma model in the chiral limit. In the relaxation 
time approximation the exponential relaxation of the pions and sigma mesons is described by the 
corresponding relaxation rates, which are in agreement with the damping rate found recently for 



the same model 99, 102 1. This particular model also reveals a crossover behavior in sigma 
meson relaxation in the large momentum limit as a result of threshold singularities associated with 
the emission and absorption of massless pions. In the relaxation time approximation we find a 
crossover between purely exponential and anomalous nonexponential relaxation with an exponent 
of the form tint. The crossover time scale is found to be dependant on the momentum of the sigma 
resonance. The anomalous relaxation is a novel result and could be of phenomenological relevance 
in view of recent suggestions of novel threshold effects of the sigma resonance near the chiral phase 
transition p^ , this possibility is worthy of a deeper study and generalization of the dynamical 
renormalization group method to reach the critical region is in progress. 

Furthermore, we have established a very close relationship between the usual renormalization 
group and the dynamical renormalization group approach to kinetics. We have shown that the 
dynamical renormalization group equation is the quantum kinetic equation, the collision terms are 
the equivalent of the beta functions in the Euclidean renormalization group. Fixed points of the 
dynamical renormalization group are identified with stationary solutions of the kinetic equation and 
the exponents that determine the stability of the fixed points are identified with the relaxation rates 
in the relaxation time approximation. We have also suggested that in this language coarse-graining 
is the equivalent to neglecting irrelevant couplings in the Euclidean renormalization program. This 
identification brings a new and rather different perspective to kinetics and relaxation that will 
hopefully lead to new insights. 

There are many advantages in the dynamical renormalization group approach to quantum kinet- 
ics: 

(i) This method is based on straightforward quantum field theoretical diagrammatic perturbation, 
hence it allows a systematic calculation to any arbitrary order in perturbation theory. It includes 
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nonequilibrium medium effects through a consistent resummation of the secular terms, thereby 
providing nonequihbrium generalizations of the usual hard thermal loop resummation in quantum 
kinetic theory. This feature is worked out in detail in the scalar field theory. 

(ii) It allows a detailed understanding of crossover behavior between different relaxation phenom- 
ena directly in real time. This is important in the case of wide resonances where threshold effects 
may lead to nonexponential relaxation on some time scales, and also in the case of relaxation near 
phase transitions where soft collective excitations dominate the dynamics. 

(iii) It describes nonexponential relaxation directly in real time whenever threshold effects are 
important, thus providing a real-time interpretation of infrared divergent damping rates in gauge 
field theories. We consider this one of the most valuable features of the dynamical renormalization 
group which makes this approach particularly suited to study relaxation in gauge field theories in a 
medium where the emission and absorption of soft gauge fields typically lead to threshold infrared 
divergences. 

(iv) This method provides a simple and natural resolution of pinch singularities often found in 
nonequilibrium field theory when the distribution functions are out of thermal equilibrium. Pinch 
singularities are the temporal Fourier transform manifestation of the secular terms, and their res- 
olution is achieved via consistent resummation of secular terms implemented by the dynamical 
renormalization group. 

We envisage several important applications of the dynamical renormalization group method 
primarily to study transport phenomena and relaxation of collective modes in gauge theories where 
infrared effects are important, as well as to study relaxation phenomena near critical points where 
soft collective fluctuations dominate the dynamics. An important aspect of this method is that 
it does not rely on the quasiparticle approximation and allows a direct interpretation of infrared 
phenomena directly in real time. We will return to this issue in Chap. I, where we study in detail 
nonequilibrium dynamics of photons and fcrmions in a QED plasma at high temperature. 



Chapter 5 



Nonequilibrium Dynamics in a 
QED Plasma at High Temperature 



5.1 Introduction 

The study of nonequilibrium dynamics in hot Abehan and non-Abehan plasmas under extreme 
conditions is of fundamental importance in the understanding of the formation and evolution of a 
novel phase of matter, the quark-gluon plasma (QGP) 1^, ||, expected to be produced in current 
ultrarelativistic heavy ion experiments at BNL Relativistic Heavy Ion Collider (RHIC) and future 
programs at CERN Large Hadron Collider (LHC). Estimates based on energy deposited in the central 
collision region at RHIC energies for y^SNN ^ 200 GeV suggest that the lifetime of the deconfined 
QGP is of order 10 fm/c At such unprecedentedly short time scales, an important aspect is an 

assessment of thcrmalization time scales and the potential for nonequilibrium effects associated with 
the rapid expansion and finite lifetime of the plasma and their impact on experimental observables. 
Lattice QCD is simply unable to deal with these questions because simulations are restricted to 
thermodynamic equilibrium quantities, hance a nonequilibrium field-theoretical approach is needed 
for an consistent description of the formation and evolution of the QGP. 

An important and pioneering step in this direction was undertaken by Geiger who applied 
transport methods combined with perturbative QCD (pQCD) to obtain a quantitative picture of 
the evolution of partons in the early stages of formation and evolution of the plasma. The consistent 
study of the evolution of partons in terms of pQCD cross sections that include screening corrections 
to avoid the infrared divergences associated with small angle scattering lead to the conclusion that 
quarks and gluons thermalize on time scales about 1 fm/c [108, 7^. The necessity of a deeper 



understanding of equilibrium and nonequilibrium aspects of the quark-gluon plasma motivated an 
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intense study of the Abelian and non-Abclian plasmas in extreme environments. A major step 
towards a consistent description of nonperturbative aspects was taken by Braaten and Pisarski , 
who introduced a novel resummation method that reorganizes the perturbative expansion in terms 
of the effective degrees of freedom associated with collective modes rather than bare particles. This 
hard thermal loop (HTL) resummation program is now at the heart of most treatments of equilibrium 
Abelian and non- Abelian plasmas at high temperature psj . 

Thermal field theory provides the tools to study the properties of plasmas in equilibrium p2| , 
|4^ ], but the consistent study of nonequilibrium phenomena in real time requires the methods of 
nonequilibrium field theory p9[| . The study of the equilibrium and nonequilibrium properties of 
Abelian and non- Abelian plasmas at high temperature as applied to the QGP has as ultimate goal 
to obtain a deeper understanding of the potential experimental signatures of the formation and 
evolution of the QGP in ultrarelativistic heavy ion collisions. Amongst these, photons and dilcptons 
(electron and/or muon pairs) produced during the early stages of the QGP are considered as some 



of the most promising signals |109-115|. Since photons and lepton pairs interact electromagnetically 
their mean free paths are longer than the estimated size of the QGP fireball ^ 10 — 20 fm and unlike 
hadronic signals they do not undergo final state interactions. Therefore photons and dileptons 
produced during the early stages of QGP carry clean information from this phase. 

In this chapter we aim to provide a comprehensive study of several relevant aspects of the 
nonequilibrium dynamics in a QED plasma at high temperature directly in real time. This study is 
of phenomenological importance in ultrarelativistic heavy ion collisions, as many features of a QCD 
plasma at high temperature are similar to those of a QED plasma at high temperature p3|, |62|, |5^ . 
In particular, the leading contributions in the hard thermal loop approximation can be straight- 
forwardly generalized from QED to QCD, thus relaxation of quarks and production of photons in 



the QGP 1 112, 113 1 can be understood to leading order from the study of a QED plasma at high 
temperature. More specifically, we focus on the following nonequilibrium aspects. 

(i) The real-time evolution of gauge mean fields in linear response in the HTL approximation. 
The goal here is to study directly in real time the relaxation of (coherent) gauge field configurations 
in the linearized approximation to leading order in the HTL approximation. Whereas a similar study 
has been carried out in scalar quantum electrodynamics (SQED) js^ and confirmed numerically in 
Ref. ||7^, the most relevant case of fermionic QED has not yet been studied in detail. 

(ii) The quantum kinetic equation that describes the evolution of the distribution function of 
photons in the medium, again to leading order in the HTL approximation. This aspect is relevant 
to study photon production via off-shell effects directly in real time. As explained in detail, this 
quantum kinetic equation, obtained from a microscopic field theoretical approach based on the 
dynamical renormalization group displays novel off-shell effects that cannot be captured via the 
usual kinetic description that assumes completed collisions [ 11C| , 112-114]. 
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(iii) The evolution of fermion mean fields at large times features anomalous relaxation arising 
from the emission and absorption of magnetic photons (giuons) which are only dynamically screened 



by Landau damping |116, 117|. The fermion propagator was studied previously in real time in the 



Bloch-Nordsieck approximation which provides a resummation of the infrared divergences associated 



with soft photon (or gluon) bremsstrahlung in the medium [116, 117 1. In this chapter we implement 
the dynamical renormalization group to study the evolution of fermion mean fields providing an 
alternative to the Bloch-Nordsieck treatment. 

(iv) We obtain the quantum kinetic equation for the fermionic distribution function for hard 
fermions via the implementation of the dynamical renormalization group. There has recently been 
an important effort in trying to obtain the effective kinetic (Boltzmann) equations for hard charged 
(quasi)particles |^,|8^,Q but the collision kernel in this equation features the logarithmic diver- 
gences associated with the emission and absorption of soft magnetic photons (or giuons) . The 
dynamical renormalization group leads to a quantum kinetic equation directly in real time bypassing 
the assumption of completed collisions and leads to a time- dependent collision kernel free of infrared 
divergences. 

A fundamental issue that must be addressed prior to setting up our study is that of gauge 
invariance. In the Abelian case, it is straightforward to reduce the Hilbert space to the gauge 
invariant subspace and to define gauge invariant field operators for the gauge boson and fermion. 
The description in terms of gauge invariant states and operators is best achieved within the canonical 
formulation which begins with the identification of the canonical fields and conjugate momenta as 
well as the primary and secondary (first-class) constraints associated with gauge invariance. The 
physical states are those annihilated by the constraints and physical operators are those commute 
with the constraints. Such a gauge invariant formulation has been implemented explicitly in the 



case of SQED [ 118 | and the fermionic case can be treated in the same manner with a few minor 
technical modifications. 

The final result of this formulation is that the Hamiltonian written in terms of gauge invariant 
fields and acting on the gauge invariant states is exactly equivalent to that obtained in Coulomb 
gauge, which is the statement that Coulomb gauge describes the theory in terms of the physical de- 
grees of freedom. Furthermore, the instantaneous Coulomb interaction can be traded for a Lagrange 
multiplier leading to the following Lagrangian density 

C[At,Ao,^,^ - i [(a,,AT)' + (VAo)'] +§(i^-e7"Ao + e7-AT-m)* 

-I- Jt ■ At + fy* + *?7, (5.1) 

where e is the gauge coupling constant, ^! is charged fermion field, At is the transverse component 
of the gauge (photon) field, and Aq is not to be interpreted as the time component of the gauge filed 
but is a Lagrange multiplier associated with the instantaneous Coulomb interaction. We emphasize 
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that the fields 5", A, and Aq are ah gauge invariant (see Refs. ]86, |118 | for details). In writing 



the above Lagrangian density, we have included the external fermionic (Grassmann) source rj and 
electromagnetic source to provide an initial value problem for studying relaxation of the mean 
fields^ 

Furthermore, we will consider a charge neutral plasma with zero fermion chemical potential at 
a temperature T 3> to, where m is the (physical) fermion mass at zero temperature, hence in what 
follows we neglect the fermion mass unless otherwise stated. 



The chapter is organized as follows. In Sec. we first study relaxation of the photon mean 



field in the hard thermal loop approximation for soft k < eT and semihard eT <^ k <ti T photon 
momentum and obtain the quantum kinetic equation for the (hard and semihard) photon distribution 



function. In Sec. 5.3 we study relaxation of fermionic mean fields for hard momentum and the 



quantum kinetic equation for the fermion distribution function. Our conclusions and some further 



questions are presented in Sec. 5.4 



5.2 Photons Out of Equilibrium 
5.2.1 Relaxation of the gauge mean field 



We begin with the relaxation of soft photons of momentum k ~ eT. As discussed in Sec. 2.2, the 
equation of motion for the transverse photon mean field can be derived from the tadpole method [^5| 
by decomposing the full quantum fields into c-number expectation values and quantum fluctuations 
around the expectation values: 

A|(x,i) = aT(x,f) +.4^(x,t), with (^(x,t))=0. (5.2) 

In momentum space, one obtains the following linearized equation of motion 

{df + aT(k, t)+ f dt' UriK t - t') ariK t') = Jt{K t), (5.3) 

where IlT(k,t — t') is the transverse part of the retarded photon self-energy and 

aT(k,i) = j d^xe-'^ ''aTix,t). 



Using the nonequilibrium Feynman rules and the free real-time fermion propagators [see Eq. (|3.6| )], 
we find 



/ + 00 
dui pTiuJ,k)sm[Lu{t~t')], (5.4) 
-oo 



^In this chapter we will not discuss relaxation of the Lagrangian multiplier Aq associated with the instantaneous 
Coulomb interaction, hence the corresponding external source is neglected. 
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where priiOjli) is the spectral function of the transverse photon self-energy and to one- loop order 
reads 



PT(w,k) 



-{[5{^-p- q)-5{^+p + q)][l + (k • p)(k • q)][l - npip) - npiq)] 



+ [Sicj-p + q)-5{uj+p- q)][nFiq) - riF(p)][l - (k • p)(k • q)]}, 



(5.5) 



with p ~ k + q. 

As before, the source Jt(x, t) is taken to be switched on adiabatically from t = —oo and switched 



off at t = to provide the initial conditions [see Eq. ( 2.22 )] 

aT(k, t = 0) = tp{k, 0), ariK t<0) = 0. 

Introduce an auxiliary quantity 7rT(k, t — t') defined by 

d 



(5.6) 



UriKt-t') = —7TT{k,t-t'), 



or, equivalently, 



/ d 
— PT(t^,k) cos\uj(t - t')], 
-oo ^ 

and impose JT(k, t > 0) = 0, we can rewrite the equation of motion for t > as an initial value 
problem 

{df + k^)aT{k,t)- [ dt'TTT{k,t-t')aT{k,t')+TTT{k,0)aT{k,t):^0, (5.7) 

where the initial conditions specified by Eq. ( |3.6[ ). 

This equation of motion can be solved by Laplace transform as befits an initial value problem. 
In terms the Laplace transformed quantities 

aT(s,k)= / dte~'* aT{k,t), nT{s,k) = / e""* 7rT(k, i), (5.8) 
Jo Jo 

where Re s > 0, the Laplace transformed equation of motion is given by 

s"^ + k'^ + IlT{s,k)] aT(s,k) = [s-^T(s,k)]aT(k,0), (5.9) 
where IIt{s, k) is the Laplace transform of nT(k, t — t') 

UT{s,k) = ttt{KO) - snT{s,k). 



Eq. (5.9) has the solution 

(s , k) = i { 1 - 5t (s , k) [fc2 TTT (k, 0)] } aT (k, 0) , 
where the retarded transverse photon propagator DT{s,k) is given by 

DTis,k) = [s^ +e +IlTis,k)]-\ 



(5.10) 



(5.11) 
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The real-time evolution of aT(k, t) is obtained by performing the inverse Laplace transform along 
the Bromwich contour which is to the right of all singularities of ax(s,k) in the complex s-plane. 
We note that this result is rather different from that obtained in Ref. |86| for SQED in the 



structure of the solution, in particular the prefactor 1/s in Eq. (5.10). This can be traced back to 
the different manner in which we have set up the initial value problem, as compared to the case 
studied in The adiabatic source (2.21) (for to = 0) has a Fourier transform that has a simple 
pole in the frequency plane, this translates into the 1/s factor in Eq. ( 5.10| ). This particular case 
was not contemplated in those studied in Ref. js^ since the source ( ^.21 ) is not a regular function 
of the frequency. This difference will be seen to be at the heart of important aspects of relaxation 
of the mean field condensate in the soft momentum limit as discussed in detail below. 



It is straightforward to see that the residue vanishes at s = in the solution ( 5.10 ), hence the 
singularities of Siris, k) are those arising from the retarded transverse photon propagator Dt{s, k). 



Soft photons k ~ eT: real-time Landau damping 

For soft photons of momenta k ^ eT, the leading (~ e^T^) contribution to the photon self-energy 
arises from loop momenta of order T and is referred to as the hard thermal loop (HTL) in the liter- 
ature ||6^ . In this region of momenta, the contribution of the hard fermion loop is nonperturbative 
(see below). 

In the HTL approximation {s ^ k <ti q), after some algebra we find that the leading contribution 
to 11^(5, k) arises exclusively from the terms associated with 6{u! ^p±q) in Eq. ( |5.5| ) and reads 



nT(s,k) = 



12 



^(l + ^lln. --2t^ 



(5.12) 



The delta functions 6{u! ^ p zL q) have support below the light cone (w^ < k^) and correspond to 
Landau damping processes in which the soft photon scatters a hard fermion in the plasma. 
The analytic continuation of 117^(5, k) in the complex s-plane (physical sheet) is defined by 



nT(t^,k) = Uri-s 



0+,k) = RenT(w,k) +amnT(w,k), 



(5.13) 



where the real and imaginary parts are related by the usual dispersion relations. The analytical 
continuation of Dt(s, k) and ttt(s, k) can be defined analogously, and they are related to nT(w, k) 

by 

AT(t^,k) = -[uj^-k^- nT(w, k)]-\ 
RenT(a;, k) — 7rT(k, Q) — lo Iva-KTi^, k), 

ImnT(t^,k) = a;Re7rT(a;,k). (5-14) 



One finds from Eq. (5.12) that the analytical continuation of n7'(s,k) in the HTL approximation is 
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given by 



RenT(w,k) = 



ImnT(w, k) 









1 In 




















LU — k 





12 



(5.15) 



We note that Imn7'(a;,k) only has support below the light cone and vanishes linearly as lo 



from below. As can be seen from Eq. (5.15) for soft momenta k < eT the HTL photon self-energy 
is comparable in magnitude to or larger than the free inverse propagator, therefore we have to treat 
it nonperturbatively for soft photons. 

Isolated poles of AT(w,k) in the complex w-plane correspond to quasiparticles (or collective 
excitations) |Q. The transverse photon pole LOT{k) is real and determined by|^ 



u)l{k) - fc2 - RenT(wT(fc), k) = 0. 



(5.16) 



For ultrasoft photons k ^ eT, the dispersion relation of the collective excitations is 



uj^ik) =Lol + -k^ + 



fc4 

.2 y2 



(5.17) 



where u)p = eT/3 is the plasma frequency. Consequently, the retarded transverse photon propagator 
Dt{s, k) has two isolated poles at s = ±ia;T(fc) and a branch cut from s = —ik to s = ik. 

Having analyzed the analytic structure of Dt{s,]<.), we can now perform the integral along the 
Bromwich contour to obtain the real-time evolution of slt (k, t) . Closing the contour in the half-plane 
with Res < 0, we find arO^-, t) for t > to be given by 

aT(k,i) =a?,°''=(k,t)+ar(k,i). 



with 



where 



a?,°'^(k,i) 



aT*(k,i) 



Zrik) = 



cos[wT(fc)i] aT(k, 0), 



k^Zrik) 
Lo'^{k) 

fc2 /■^'^;3T(c^,fc)e-^"*aT(k,0), 
-k w 



1 - 



9RenT(w,k) 



/3T{io,k) = - 



ImnT(w, k) 



(5.18) 
(5.19) 



(5.20) 



TT [^2 _ fc2 _ RenT(w, k)]^ + [ImnT(t^, k)]-" 
The solution evaluated at t = must match the initial condition, this requirement leads immediately 
to the sum rule Bsl 



dqo - , , 1 
/ PT{qo,q) = ^, 



^In the following discussion, uj'r{k) is referred to as the positive pole. 
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where pT{qo,q) is the HTL-resummed spectral function for the transverse photon propagator]^ 
PT{qo,q) = -lmAT{qo,q) 

TT 

= Sgn(go) Zriq) S[ql ~ LU^q)] + /3t(<7o, q) 0{q^ - ql). (5.21) 

A noteworthy feature of the cut contribution ( 5.19| ) is the factor lo in the denominator. The presence 
of this factor can be traced back to the preparation of the initial state via the adiabatically switched- 
on external source, which is the switched off at t = 0. The Fourier transform of this source is 



proportional to 1/uo and results in the prefactor of Prii-^, k) in Eq. (5.19). 



For k ~ eT, Eq. (5.19) cannot be evaluated in closed form but its long-time asymptotics can be 
extracted by writing the integral along the cut as a contour integral in the complex a;-plane. The 
integration contour C is chosen to run clockwise along the segment — fc < a; < A; on the real axis, 
the line lu — k — iz with < z < oo, then around an arc at infinity and back to the real axis along 
the line lu = —k — iz. After some algebra we obtain the following expression 

,-»fct f°° dze-^* iPT{k-iz,k) 



'(k,t) 



27ri Res 

poles 
inside C 



k — iz 

f ^T(^,fc) ^^,^f 



c.c. 



LO 



aT(k,0). 



(5.22) 



The contribution from the poles inside the contour C is dominated at long times by the pole closest 
to the real axis and is exponentially suppressed. On the other hand, because of the exponential 



factor 



in the integrals, the dominant contributions to the integral at long times {t ^ 1/fc) 



arise from the end points of the branch cut with z <C fc, thereby leading to a long-time behavior 
characterized by a power law: 



af,"*(k,t) 



fet»i 



12 cosfci 



aT(k,0) 



l + O 



for k ~ eT, 



(5.23) 



g2y2 ^2 

which agrees with the results of Ref. in scalar quantum electrodynamics (SQED) at high tem- 
perature. 

For very soft momenta k <^ eT, the function (3t{(^, k)/u) is strongly peaked at a; = as depicted 



in Fig. 5.1. Furthermore, for lo <^k <^ eT the function I3t{uj, k)/uj features a Breit-Wigner form 

r(fc) 



-PT{i^,k) = -r2 2 , r2fl,\ 
LO Trfc^ LO'' + I ^(k) 



1 12^;"^ 

with T{k) = —. 

TT e^J ^ 



Upon using this narrow- width approximation in evaluating Eq. ( ^.19 ), one obtains 



'(k,i) 



(it) 



aT(k,0) for fc < eT. 



(5.24) 



(5.25) 



The end-point contribution which results in a power law as in Eq.(5.23) is very small compared with 



Eq. (5.25) for times t < l/r(A;). This power law becomes dominant on time scales much longer than 

^In our notation the spectral function for the self-energy is denoted as p , whereas that for the corresponding 
propagator is denoted as p. 
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Figure 5.1: The function /9t(^, k)/u! plotted as a function of u; for ultrasoft photons [k/eT = 0.1). 



l/r(fc), because its amplitude a,t t ^ l/r(fc) is of order (k/eT)^ <C 1 in the soft momentum limit. 
Thus, upon combining the pole and cut contributions, we find the following real-time behavior 

'PZrik) 



aT(k,t) MKO) 



-2^^ cos[ujT{k)t] + e-^C^)* 



for k < eT. (5.26) 



uj^{k) 

Two comments here are in order: (i) The exponential decay is not the same as the usual decay 
of an unstable particle (or even coUisional broadening) for which the amplitude relaxes to zero at 
times much longer than the relaxation time. In this case the asymptotic behavior of the amplitude 
is completely determined by the transverse photon pole Lorik), and to this order in the HTL ap- 
proximation the collective excitation is stable. The exponential relaxation does not arise from a 
resonance at the position of the transverse photon pole but at zero frequency. Clearly, we expect 
a true exponential damping of the collective excitation at higher order as a result of coUisional 
broadening, (ii) The exponential relaxation given by Eq. ( 5.25| ) can only be probed by adiabatically 



switching on an external source [see Eq. (2.21)] whose temporal Fourier transform has a simple pole 
at w = which is the origin of the prefactor 1 /u) in Eq. (|5T9| ). That is the adiabatic preparation of 
the initial state excites this pole in the Landau damping cut. For external sources whose temporal 
Fourier transform is regular at cj = no exponential relaxation arises, in agreement with the results 
of Ref . ||86|] . Thus we find that the exponential relaxation is not a generic feature of the evolution 
of the mean field, but emerges only for particular (albeit physically motivated) initial conditions. 

Our results confirm those found in SQED at high temperature [f76| , wherein a numerical analysis 
of the relaxation of the mean field was carried out. In that reference the authors studied the real 
time evolution in terms of local Hamiltonian equations obtained by introducing a nonlocal auxiliary 
field. The initial conditions chosen there for this nonlocal auxiliary field correspond precisely to the 
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choice of an adiabatically switched-on external source leading to an adiabatically prepared initial 
value problem. 



Semihard photons eT <^ k ^ T: anomalous relsLxation 

Most of the studies of the photon self-energy in the hard thermal loop approximation focused on the 
soft external momentum region k <C eT (with e <C 1) [ p3| , |6^ ,p3[ wherein the contribution of the hard 
thermal fermion loop must be treated nonperturbatively. However, there are important reasons that 
warrant a consideration in the region of semihard photons (eT <^ k <^ T). In particular, from a 
phenomenological standpoint, hard and semihard photons produced in the QGP phase are important 



electromagnetic probes of the deconfined phase |11C, 112 -115|, hence a study of kinetics, relaxation 
and production of photons all over the spectrum to explore potential experimental signatures is 
warranted. 

A more theoretical justification of the relevance of the semihard region is that whereas the 
photon self-energy for this region of momentum is still dominated by the hard thermal fermion loop 



contribution given by Eq. (5.12), its contribution to the full photon propagator is now perturbative as 
compared with the free inverse propagator. The validity of perturbation theory in this regime allows 
us to study the real-time evolution with the tools developed in Chap. |^ that provide a consistent 
implementation of the dynamical renormalization group to study nonequilibrium phenomena directly 
in real time. 



As we have shown in Sec. 4.5, the dynamical renormalization group method is particularly suited 



to study the real-time evolution in the case in which there are (infrared) threshold singularities in the 
spectral function. To understand the potential emergence of anomalous thresholds in the semihard 
and hard limit for the photon self-energy, we note that at leading order in the HTL approximation 



the Laplace transform of the inverse photon propagator can be written as [see Eqs. ( 5.11 ) and ( 5.12 )] 



1 



e^T^ is is 
■ In — 



6 



(5.27) 



12 fc3 is - k 

The transverse photon poles are at s = ±i[k + e^r^/12fc + 0{T^/k'^)], which in the limit k ^ eT 
becomes s zLik. In this limit the poles merge with the thresholds of the logarithmic branch 
cut arising from Landau damping and are no longer isolated from the continuum, leading to the 



enhancement of the spectral function near the threshold as clearly displayed in Fig. 5.2. While the 
perturbative expansion in the effective coupling e^T^/fc^ is warranted in the semihard case under 
consideration, the wave function renormalization constant evaluated as the residue at the pole is a 
nonanalytic function of this coupling and given by 



Zrik) = 1 



e'T' 
12k^ 



3 + In 



e'T' 
24fc2 



O 



e^T\ eT 



(5.28) 



This is obviously a consequence of the logarithmic threshold singularities in the semihard regime. 
This threshold singularity is reminiscent of those studied in detail in Ref. p4] , where it was shown 
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Figure 5.2: The function i3t{(^, k)/uj plotted as a function of uj for semihard photons {k/eT = 10). 

that the Bloch-Nordsieck resummation, which is equivalent to a renormalization group improvement 
of the space-time Fourier transform of the self-energy, leads to a real-time evolution that is obtained 
by the implementation of the dynamical renormalization group in real time . We now implement 
the dynamical renormalization group program to study the relaxation of the photon mean field in 
the semihard limit directly in real time. 

Perturbation theory in terms of the effective coupling e^T^/Zc^ is in principle reliable for semihard 
photons, hence we can try to solve Eq. { p.7\ ) by perturbative expansion in powers of (the true 
dimensionless coupling is e^T^/fc^). Writing 

aT(k, t) = a^°^ (k, t) + ai^^ (k, t) + O(e^), 

7rT(k, t^t')^ (k, t-t')+ O(e^), (5.29) 
and expanding Eq. (^^) consistently in powers of e^, we obtain a hierarchy of equations: 
(92 + fc2)a?)(k,t) = 0, 

{d^ + k^) a^r\'k,t) = / di'4'^(k,t-t')4°^(k,t'), (5.30) 

where we have used the fact that 7r7'(k, i = 0) = in the HTL approximation. Using the solution 
to the zcroth-order equation 

a^\k,t) = ^ [AA(k)e-*^-* + AKk)e*'=*] ^^(k), (5.31) 

A=l,2 

where ^^A(k) is the polarization vector, and the retarded Green's function of the unperturbed problem 

Gg)(fc,t-0- '^"^^^!"''^^ Oit-t'), (5.32) 
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one finds the solution to tlie first-order equation reads 



duj 

— oo 



UJ 



A=l,2 



2ik 




i{uj ^ k) 

where pt(i^, k) is the spectral function of tlic photon self-energy in the HTL approximation 

1 



(5.33) 



PT(w,k) 



-lmIlT{LO, k) 



12 k 



(5.34) 



Potential secular terms arise at long times from the regions in u where the denominators in 
Eq. ( 5.33| ) are resonant. They can be extracted in the long-time limit by using the following formu- 
las Isl: 



y 



y 



dy , X / N 

— (1 - cos yt)p{y) 



t->oo TT 



tp(0)+p'(0) [ln(Mt)+7] 



sin yt 



p{y) 



+ - ^(0) - yp'iO) e{p^ -y)] + (t-i) 

■/o y 

tK0)[ln(A*t)+7-l] 

+ tp [p{y) ~ p{0) 9{pi -y)] + 0{t-^), 

Jo y 



(5.35) 



where 7 = 0.577 ... is the Euler-Mascheroni constant and as can be easily shown the dependence 
on the arbitrary scale /i in the above integrals cancels. For our analysis, the integration variable 
y — Lu ± k and p{y) = pTiuj,'k.)/uj. The terms that grow linearly in time are recognized as those 
emerging in Fermi's golden rule from elementary time-dependent perturbation theory in quantum 
mechanics. We note that p{0) = and p'(0) 7^ 0, thus the first integral above gives a contribution to 
the real part of the mean field (the amplitude) that features a logarithmic secular term, whereas the 
second integral contributes to the imaginary part (the phase) with a linear secular term, which as will 
be recognized below determines a perturbative shift of the oscillation frequency. Upon substituting 



Pt(w, k) into Eq. ( 5.33 ), one obtains 



kt 



12k' 



{\n2kt + j-l) + iSkt + 0{t-^) 



c.c. 



where 



_ RenT(w,k) 



2uj 



12k 



(5.36) 



(5.37) 



While the linear secular terms have a natural interpretation in terms of renormalization of the 
mass (the imaginary part) and a quasiparticle width (the real part) pi], the logarithmic secular 
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term found above is akin to those found in Rcf. |84| that lead to anomalous relaxation. Furthermore 
the origin of these logarithmic secular terms is similar to the threshold infrared divergences and 
threshold enhancement of the spectral function due to the presence of nearby poles studied in 



Sec. 4.5. Thus we implement the dynamical renormalization group by introducing a (complex) 
amplitude renormalization factor in the following manner, 



Ax{k)^ZkiT)Axik,T), Zfc(r) = l + e^ZkiT)+0{e^ 



(5.38) 



where Zk{T) is a multiplicative renormalization constant, .^^(k, r) is the renormalized initial value, 
and T is an arbitrary renormalization scale at which the secular divergences are cancelled p^ . 
Choosing 



e^Zfc(T) 



12fc2 



(ln2fcT + 7-l) + i4 r. 



we obtain to lowest order in e T /k that the solution of the equation of motion is given by 



aT(k,i) = -4A(k,r)£,(k)e-*'=* 
+ nonsecular terms + c.c. 



12A:2 T 



(5.39) 



which remains bounded at large times t provided that t is chosen arbitrarily close to t. The solution 
does not depend on the renormalization scale r and this independence leads to the dynamical 
renormalization group equation [ p4[ , which to this order is given by 

'd_ e^T^ 



+ iSk 



^A(k,T) =0, 



(5.40) 



with the solution ^ ^ 

Ax{k,T)=Ax{k,To)e-''-^^-^'>^(^^^ ^ , (5.41) 
where tq is the time scale such that this intermediate asymptotic solution is valid and physically 



corresponds to a microscopic scale, i.e, tq ^ \/k. Finally, setting r = t in Eq. (5.39) we find that 
a7'(k, f) evolves at intermediate asymptotic times t 3> 1/fc as 



aT(k,i)~ V ^A(k,To)£A(k)e-''(^+^'=)(*-^°) (- 

^ , . V To 



C.C. 



(5.42) 



From Eq. ( 5.37 ) we see that 5k is consistent with the photon thermal mass = e^T^/Q for 
fc2 » m2 II. 

As discussed in detail in Ref. jsi) the dynamical renormalization group solution (|5.42| ) is also 
obtained via the Fourier transform of the renormalization group improved propagator in frequency- 
momentum space, hence the above solution corresponds to a renormalization group improved re- 
summation of the self-energy. 

This novel anomalous power law relaxation of the photon mean field will be confirmed below in 
our study of the kinetics of the photon distribution function in the linearized approximation. We 
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note that this anomalous power law relaxation is obviously very slow in the semihard regime in which 
the HTL approximation and perturbation theory is valid. At higher orders we expect exponential 
relaxation due to coUisional processes, which emerges from linear secular terms ]8^ ] in a perturbative 
solution of the real-time equations of motion. The power law relaxation will then compete with the 
exponential relaxation and we expect a crossover time scale on which relaxation will change from a 
power law to an exponential. Clearly an assessment of this time scale requires a detailed calculation 
of higher order contributions which is beyond the scope of this thesis. 

A related crossover of behavior will be found below for the evolution of the photon distribution 
function and photon production. 



5.2.2 Quantum kinetics of photons 

As mentioned in Chap. ^ the first step towards a kinetic description is to identify the proper degrees 
of freedom (quasiparticles) and the corresponding microscopic time scale. For semihard photons of 
momentum eT ^ fc <C T, it is adequate to choose the free photons as the quasiparticles with the 
corresponding microscopic scale 1/k. A kinetic description of the nonequilibrium evolution of the 
distribution functions assumes a wide separation between the microscopic and the kinetic time scales, 
which is justified in the weak coupling limit under consideration. In the semihard momentum regime 
the effective small coupling is e'^T^/k'^ ^ 1 and both the HTL and the perturbative approximations 
are valid. Furthermore, we assume that the fermions are in thermal equilibrium at a temperature 
T and that there is no initial photon polarization asymmetry. As before, we consider the case in 
which the initial density matrix is diagonal in the basis of the photon occupation numbers but with 
nonequilibrium initial photon distribution functions n^ito). 

We begin by obtaining the photon number operator from the Heisenberg photon field operator 
and its conjugate momentum (in momentum space) 



X=l,2 

UriKt) = zJ2\[l[4i-^^t)exi^k)-ax(k,t)£x(k)\, (5.43) 



where ax(k,t) [a^(k, t)] is the annihilation (creation) operator that destroys (creates) a free photon 
of momentum k and polarization A at time t and Sxik) is polarization vector. The number operator 
Ny (k, t) that counts the (semihard) photons of momentimr k is given by 



A=l,2 



-k, t) ■ nT(k, t) + t=AT(-k, . AT(k, t)] - ^. (5,44) 
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whose expectation value is interpreted as the number of photons per unit phase space volume 

nl{t)^{N,{k,t)) = {2nf^^, (5.45) 

where N is the total number of (semihard) photons in the plasma. Using the Heisenberg equations 
of motion, we obtain to leading order in e the following expression in the CTP formalism 

l-kW-lim^A J _^(^-(_p,t)^.A+(k,t')^-(q,t))+c.c., (5.46) 

where we have separated the time arguments to extract the time derivative from the expectation 
value. 

The nonequilibrium expectation values can be computed perturbatively in powers of e using the 
nonequilibrium Feynman rules and real-time propagators. At 0{e) the right-hand side of Eq. ( ^.46 ) 
vanishes identically. This is a consequence of our choice of initial density matrix diagonal in the 
basis of the photon number operator. To 0{e^), we obtain 

d 

where the r^''''''(t) is the time- dependent photon production (absorption) rate 



^y^{t) = [l + nl{to)]r<{t)-nl{t,)T>{t), (5.47) 



rm ^ P^c.7^f(.,fc) -^"[(;-^)(^-^°^ (5.48) 



and 

^2 r rf3 



2 [1 - (k • p)(k • q)] npiq) [1 - upip)] S{cj+p- q) 



(27r)3 

+ [1 + (k • p)(k • q)] ^^nFip)nFiq) 6(0; ^ p - q) + [1 - npip)] 
X [1 - nF{q)]S{uj +p + q)^ 



(5.49) 



In the above expression, p = k + q and TZ^ {uj, k) is obtained from TZ^ {lu, k) through the replacement 
np <^ ^ — np- A comment here is in order. As explained above we are focusing on the leading HTL 
approximation, consequently, in obtaining TZj-{uj,k) we use the free real-time fermion propagators 
which correspond to the hard part of the fermion loop momentum. In general there are contributions 
from the soft fermion loop momentum region which will require to use the HTL-resummed fermion 
propagators | |6^ for consistency. A detailed study of the contribution from soft loop momentum is 
beyond the scope of this thesis, instead we focus here on the lowest order leading HTL contribution 
in real time. 

Since the fermions are in thermal equilibrium, the Kubo-Martin-Schwinger (KMS) condition 
holds: 

n>{Lo,k)=e^'^n<{Lo,k), (5.50) 

where f3 — 1/T. It is easy to recognize that T?.^ (tj, fc) has a physical interpretation in terms of 
the off-shell (energy-nonconserving) photon production (absorption) processes in the plasma. The 
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first term in TZ^{uj, k) describe the process that a fcrmion (or an antifcrmion) emits a photon, i.e., 
bremsstrahlung from the fcrmions in the medium, the second term describes annihilation of a fcrmion 
pair into a photon, and the third terms describes creation of a photon and a fcrmion pair out of the 
vacuum. The corresponding terms in TZ^{uj, k) describe the inverse processes. 

As argued above, for semihard photons of momenta eT <C fc ^ T, the leading contribution of 
TZj{uj, k) arises from the hard loop momenta q ^ k. An analysis of TZ^ (cj, fc) along the familiar lines 
in the hard thermal loop program |^,|4^ shows that in the HTL approximation {q 3> k) 

7^<(-, fc)|HTL = ^ (l - ^ ) Oik' - (5.51) 

Thus we recognize that in the HTL approximation TZ^'^^^Lu^k) is completely determined by the 
off-shell Landau damping process in which a hard fermion in the plasma emits (absorbs) a semihard 
photon, i.e., bremsstrahlung from the fermions in the medium. 

Dynamical renormalization group and the emergence of detailed balance 



We now turn to the kinetics of semihard photons. To obtain a kinetic equation from Eq. ( 5.47 ), 
we implement the dynamical renormalization group resummation as introduced in Chap. ^. Direct 
integration with the initial condition yields 

nlit) = nl{t,) + [l + nl{t,)]fdt'T<{t')-nl{t,)fdt'T>{t'). (5.52) 

J Iq J to 

The integrals that appear in the above expression 

' dt' T>< it') = /"J d. n>< (., k) (5.53) 

are dominated, in the long-time limit, by the regions of lo for which the denominator is resonant, 
i.e., LO K k. The time dependence in the above integral along with the resonant denominator is 
the familiar form that leads to Fermi's golden rule in time-dependent perturbation theory. In the 
long-time limit t — Iq 3> 1/fc, Fermi's golden rule approximates the above integrals by 



to 



* 1 



dt' r< {t') ^ -{t-to)n^{uj = k,k) = o. (5.54) 



Therefore, the rate of photon production, i.e., the coefficient of the linear time dependence vanishes 
at this order because of the vanishing of the imaginary part of the photon self-energy on the photon 
mass shell. However the use of Fermi's golden rule, which is the usual approach to extract (time- 
independent) rates, misses the important off-shell effects associated with finite lifetime processes. 



These can be understood explicitly by using the first of formulas given in Eq. (5.35). We find that 
for i - io > 

[itTun «'-s'>' ;g!{,„|2M.-wi+7-i) + <:^(^) 

fdt'r>(t'). (5.55) 

Jtn 
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The second equality of this equation displays the condition for detailed balance and holds for time 
scales t — tf) ^ 1/k. We emphasize that this condition is a consequence of the fact that the region 



Lo K k (the resonant denominator) dominates the long-time behavior of the integrals in Eq. (5.53) and 



the KMS condition (5.50) holds under the assumption that the fermions are in thermal equilibrium. 

Thus we find that the time-dependent rates obey detailed balance which emerges in the interme- 
diate asymptotic regime < — to ^ when the secular terms dominate the long-time behavior. This 
is a noteworthy result because it clearly states that detailed balance emerges on microscopic time 
scales t ~ to > 1/k aX which the secular terms dominate the integrals but for which a perturbative 
expansion is still valid. Detailed balance then guarantees the existence of an asymptotic equilibrium 
solution which is reached at times of the order of the kinetic time scale. 



The (logarithmic) secular terms in the time-dependent rates (5.55) can now be resummed us- 
ing the dynamical renormalization group method by introducing a renormalization of the initial 
distribution function 

nl{to)^Zk{to,T)nl{T), Zk{to.T) = l + e^Zk{to,T)+0{e^), (5.56) 



thus rewriting Eq. ( 5.52 ) consistently to order e as 



Kit) = nl{T)+e^z^{to.r)nl{T) + [l + nl{T)] dt'T<{t') 

Jto 

-nlir) f dt'T>{t')+0{e% (5.57) 



The renormalization coefficient Zk{tf),T) is chosen to cancel the secular divergence on a time scale 
t = T. To lowest order in the choice 

e'z,it,r)nliT) = -[l + nlir)] f dt' V<{t') + nUr) f dt'V>{t'), (5.58) 

leads to an improved perturbative solution in terms of the renormalized occupation number ^^(t) 
as 

nl{t) = r^lir) + [l + nliT)]j\t'T<it')-r^lir)j\t'T>{t')+Oie^), (5.59) 

which is valid for large times t ^ Iq provided that r is chosen arbitrarily close to t. A change in 
the time scale r is compensated by a change of the ti^It) in such a manner that n^{t) does not 
depend on the arbitrary scale r. This independence of r leads to the dynamical renormalization 
group equation which, consistently to order e^, is given by 

^nlir) ^ [1 + nlir)] r< (r) - nl{r) r>(r) + 0{e^). (5.60) 



Choosing t to coincide with t in Eq. ( 5.60| ), we obtain the quantum kinetic equation to order to 
be given by 

j^nlit) = [l + nl{t)]r<{t)-nlit)r>{t). (5.61) 
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For intermediate asymptotic times < — to ^ 1/fc at which the logarithmic secular terms dominate 
the integrals for the time-dependent rates and detailed balance emerges, we find 



< fe(t-io)»i e^T^ 1 



l + O ' ^ 



---^'r>(t), (5.62) 



e 



where the detailed balance relation is explicitly displayed and the terms being neglected are os- 
cillatory on time scales ~ 1/fc and fall off faster. This detailed balance relation between the time- 
dependent rates r^(i) guarantees the existence of an asymptotic equilibrium solution of the quantum 
kinetic equation with the distribution function n^{t = oo) = l/(e^'^ — 1). 



The full solution of the quantum kinetic equation (5.61) is given by 

n;^(t) =n;^(to)e" '''o'*'^'=^*'^ +e--f4''*'^''(*') f dt'T<{t')e^^'''"^''^'"\ (5.63) 



Jto 

where Fkit) — ^^(t) — T^{t). However, in order to understand the relaxation to equilibrium of 
the distribution functions, we now focus on the relaxation time approximation, which describes the 
approach to equilibrium of a slightly off-equilibrium initial distribution function n^{to) ~ nsik) + 
dn^{to) while all other modes are in thermal equilibrium. In the relaxation time approximation, one 
obtains 



-Snlit) = -6nl(t) ieP'^ l)T<{t). (5.64) 



d 

dt' 

For semihard photons of momentum eT <^ k <^T, one can simply replace e^'^ — 1 by k/T and upon 
integration one finds 



Snlit) ^6nl{to)[k{t-to)]-''^"/'^''" for fc(t-to)»l- (5.65) 



A noteworthy feature of Eq. (5.65) is that the relaxation of the distribution function for semihard 
photons in the relaxation time approximation is governed by a power law with an anomalous ex- 
ponent rather than by the usual exponential relaxation. This result is similar to that found in 
scalar quantum electrodynamics in the Markovian approximation p6| |. Comparing Eq. ( ^.65 ) and 



Eq. (5.42), we clearly see that the anomalous exponent for the relaxation of the photon distribution 
function in the relaxation time approximation is twice that for the linear relaxation of the photon 



mean field. This relationship is well known in the case of exponential relaxation [see Sec. 3.4| but our 
analysis with the dynamical renormalization group method reveals it to be a more robust feature, 
applying just as well to power law relaxation. 

The dynamical renormalization group approach to derive quantum kinetic equation in field theory 
is different from the one often used in the literature which involves a Wigner transform, assumption 
about the separation of fast and slow variables, and quasiparticle approximation |67|, |68| [zs], 
|85| ]. In particular, this approach reveals clearly the dynamics of off-shell effects associated with 
nonexponential relaxation. This aspect will acquire phenomenological relevance in our study of 
photon production enhanced by off-shell effects from a quark-gluon plasma with a finite lifetime in 
Chap. ||. 
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5.3 Hard Fermions Out of Equilibrium 

To provide a complete real-time picture of nonequilibrium dynamics in a QED plasma at high 
temperature, we now focus on a detailed study of relaxation of fermionic mean fields (induced by 
an adiabatically switched-on Grassmann source) as well as the quantum kinetics of the fcrmion 
distribution function. 

In this section we assume that the photons are in thermal equilibrium, and since we work to 
leading order in the HTL approximation we can translate the results vis a vis to the case of equi- 
librated gluons. In particular we seek to study the possibility of anomalous relaxation as a result 



of the emission and absorption of magnetic photons. In Refs. [116, 117] it was foimd that the re- 
laxation of fermionic excitations is anomalous and not exponential as a result of the emission and 
absorption of magnetic photons that are only dynamically screened by Landau damping. The study 
of the real-time relaxation of the fermionic mean fields in these references was cast in terms of the 
Bloch-Nordsieck approximation which replaces the Dirac gamma matrices by the classical velocity of 
the fermion. In Ref. |Q the relaxation of a charged scalar mean field as well as the quantum kinet- 
ics of the distribution function of charged scalars in scalar electrodynamics were studied using the 
dynamical renormalization group, both the charged scalar mean field and the distribution function 
of charged particles reveal anomalous nonexponential relaxation as a consequence of emission and 
absorption of soft magnetic photons. While electric photons (plasmons) are screened by a Debye 
mass which cuts off their infrared contribution, magnetic photons are only dynamically screened by 
Landau damping and their emission and absorption dominates the infrared behavior of the fermion 
propagator. 

While the dynamical renormalization group has been implemented in scalar theories it has not 
yet been applied to fermionic theories. Thus the purpose of this section is twofold: (i) to implement 
the dynamical renormalization group to study the relaxation and kinetics of fermions with a de- 
tailed discussion of the technical differences with the bosonic case and (ii) to focus on the real-time 
manifestation of the infrared singularities associated with soft magnetic photons. 

5.3.1 Relaxation of the fermion mean field 

The equation of motion for a fermion mean field is obtained by following the strategy described in 
section II. We begin by writing the fermion field as 

*±(x,t) = V(x,t)+x^(x,t), with (x±(x,t))=0. (5.66) 

Then using the tadpole method W5[ with the external Grassmann source 77 (x, t) that is adiabatically 



switched-on from t = —00 and switched-off at i = [see Eq. (2.22)], we find the Dirac equation of 
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the fermion mean field for t > to be given by (in momentum space) 

{il°dt - 7 • k) V^(k, t)~ f dt' S(k, t - t') V-Ck, t') = 0, (5.67) 

J —CO 

where E(k, t — i') is the retarded fermion self-energy. As noted above the relaxation of hard fermions 
is dominated by the soft photon contributions, thus in a consistent perturbative expansion one needs 
to use the HTL-resummed effective photon propagators to account for the screening effects in the 
medium . The HTL-resummed effective photon propagators can be conveniently written in terms 
of the resummed spectral functions: 
(i) Transverse components: 

D++{ci;t,t') = D>{ci;t,t')e{t-t')+D<{(i;t,t')e{t' ~t), 
Di-{q:,t,t') = D>{q:,t,t')e{t' -t) + D<{c^,t,t')e{t-t'), 
D^+{ci;t,t') - D>ici;t,t'), D+- {ci;t,t') ^ D<ici;t,t'), 

/ + 00 
dqoPT{qo,q) [l + nB(9o)]e-*«°(*-*'), 
-OC 

/ + 00 
dqoPT{Qo,q) nB(<7o)e-^'° 
-OO 



(5.68) 



(ii) Longitudinal components: 
D++iq;t,t') = 



1 



6{t - t') + D> (q; t, t')e{t - t') + D< (q; t, t')e{t' - t), 



Dr{q:t,t') = -S{t-t') + D>{q;t,t')e{t' -t)+D<{(i;t,ty{t-t'), 
Dl+{ci;t,t') - D>{q;t,t'), D+- {ci;t,t') = D<{q;t,t'), 

+ 00 



D>{q;t,t') = ^/ dgo Pl(<?o, g) [1 + ns(go)] e""?°^*"*'\ 

J — OO 

/+00 
dqoPLiqo,q)nBiqo)e-"'"'^'~*'l 
-OO 



(5.69) 

Here Pt{qqtQ) is the HTL-resummed spectral function for transverse photon propagator defined in 



Eq. ( 5.21 ) and PLiqo,q) is the HTL-resummed spectral function for longitudinal photon propaga- 
tor El 



PL{qQ,q) = 
Phiqo.q) = 



sgn(go) ZL{q) 5[ql - coUq)] + P^qo, q) 0{q^ - <Zo), 



6 q 



q^ 



2 - *1 In 2+M 

9 9-90 



(5.70) 



where lol(ci) is the plasmon (longitudinal photon) pole and ZL{q) is the corresponding residue. We 
note that pi(2)('^,k) is an even (odd) function of oj, a property that will be useful in the following 
analysis. 

Using these HTL-resummed effective photon propagators, one finds to one-loop order E(k, t — t') 
reads 



I](k,t-<') = 



+ C30 



duj 



—i 7° pi{uj, k) cos[uj{t — t')] + 7 • kp2(w, k) sin[a;(t — t')] 



(5.71) 
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where the spectral functions 



Pi (w,k) 



+ 00 



/02(w,k) 



(2^)^.-00 
X [5{uj - p - qo) + 5{uj + p + qa)] 



(27r)3 



dqo 



k • p 

(k • q)(p • q) priqo, q) PL{qo, q) 



X [1 + ns(f?o) - npip)] [S{u;-p-qo) - 5{uj+p + qo)], 



(5.72) 



(5.73) 



with p = k — q. As before, it proves convenient to introduce the Laplace transform of the retarded 
self-energy I](s,k), whose analytic continuation in the complex s-plane (physical sheet) I](ci;,k) is 
obtained through the replacement s —ito + 0^, or more explicitly, 



S(c^,k) 



-7°Pi(fco,k) + 7 • kp2(fco,k) 



(5.74) 



, ko ~ Lu ~ iO+ 

A comment here is in order. To facilitate the study and maintain notational simplicity, in obtaining 
the S(k, t — t') we have neglect the contribution arising from the instantaneous Coulomb interaction, 
which is irrelevant to the relaxation of the fermion mean field and only results in a perturbative 
frequency shift. 

Following the same strategy in the study of the photon mean field, by introducing an auxiliary 
quantity tT(k, t — t') defined as 

S(k,t-t') = at,cr(k,i-t'), 
one can rewrite Eq. ( ^.67 ) as an initial value problem 



(i7°9t-7-k)i/;(k,t)+ / di'cr(k,<-t')^(k,i') -cr(k,0)7/'(k,t) = 0, (5.75) 

with the initial conditions ^(k, 0) — ^/'o(k) and ^/'(k, f < 0) = 0. 

We are now ready to solve the equation of motion by perturbative expansion in powers of just 
as in the case of the gauge mean field. Let us begin by writing 

V'(k, t) = (k, t) + V'(') (k, t) + 0{e^) , 
(T(k, t-t')^ a^i) (k, t~t')+ O(e^), 

we obtain a hierarchy of equations: 

(z709t-7.k)^(°)(k,i) = 0, 

(^7°9t-7•k)7/.(l)(k,^) = tT(i)(k,0)V^°Hk,t)- / di' ct^i^ (k, i - i') V-^"^ (k, t'), (5-76) 



These equations can be solved iteratively by using the zeroth-order solution 
7^(0) (k,i) = [^(k)iis(k)e-^'=* +S.(k)i;,(-k)e'''^*] , 



(5.77) 



s=l,2 
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and the retarded Green's function of the unperturbed problem 

g^°\-k, t-t') = -i [/i+(k) e-^'^f*-*') + /i_(k) e''^'(*-*')j e{t - t'). 

where /i±(k) — (7° =F 7 • k)/2 and Us(k) and Ws(k) are the free Dirac spinors that satisfy 

;i+(k)u,(k) = 0, /i_(k)w,(k) = 0. 

The solution to the first-order equation is found to be given by 



(5.78) 



(5.79) 



where 



V'^^^^k, t) 



e2V^(i'^)(k,i) = 



7"(k)t [^(k)«.(k)e-''=*-i3,(k)z;,(-k)e^'=*], 

s=l,2 

+ rl,.. ( V I _ g-i(w-fc)t 



dio 
uj — k 



7''(^,k)<^^(k)7.,(k)e- 



-ikt 



i(uj — k) 



•B,(k) z;,(-k) e' 



ifci 



i + 



i{uj — k) 



with 



7"(k) 



7^(0;, k) 



-fOO 7 1 

— P2(w,k) = -Tr ;i+(k)ReE(0,k) 



Trfc 



[pi(w,k) - P2(w,k)] 



2cj 



■Tr 



/i+(k) ImE(w,k) 



(5.80) 



(5.81) 



(5.82) 



We note that in Eq. (5.8C) secular terms are explicitly linear in time and are purely imaginary, 



whereas in Eq. ( 5.81 ) potential secular terms may arise at long times from the resonant denominators. 
From the expression of 7''(a;,k) given in Eq. ( 5.82 ), we recognize that 7''(ijj,k) evaluated at the 



fermion mass shell w = fc is the fermion damping rate computed p3| |. However, it has been shown in 
the literature that due to the emission and absorption of soft quasi-static transverse photons which 
are only dynamically screened by Landau damping, the fermion damping rate exhibits infrared 
divergences near the mass shell in perturbation theory. 

The infrared divergences is easily understood from the following analysis. For soft photons with 
<?o, 9 <C T, we can replace 



l + nB{qa) -npip) ^T/qa, p~fc-gcos( 



(5.83) 



where cos = k • q, thus writing 



7^0;, k) = 



Tk f d^q r+^dgo 



^ J (27r)3 J_q qo 
X S{uJ ~ k + qcos9 — qo), 



(1 - cos26') PT{qa,q) + /3l((Jo, g) 



(5.84) 
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where the subleading pole contributions which corresponding to emission and absorption of on-shell 



photons have been neglected |116]. Recall that for very soft q <C eT the function Pt{<1otQ)/Qo is 



strongly peaked at = [see Eq. (5.24) and Fig. and as q — > it can be approximated by 

(5.85) 



1 S( ') 

— /3t((7o < 9, g) ^ — J- as g 0. 



whose physical origin is the absence of a magnetic photon mass in the plasma. The infrared diver- 
gences near the fermion mass shell become manifest after substituting Eq. ( 5.85D into Eq. (5.84). 



In order to isolate the singular behavior of 7 (a;,k), we follow the steps in Rcf. [116| and write 



— /^T(go,g) = 6{qo) [\ 2 1^ 2 



1 



(5.86) 



where up = eT/3 is the plasma frequency and VT{qQ,q) denotes the regular part of the transverse 



photon spectral function. Substituting Eq. ( 5.86 ) into Eq. ( 5.84 ), we can then separate j^{uj, k) into 
an infrared singular part which is logarithmically divergent near the fermion mass shell 



7sins(^,k) 



(27r)3 



q"^ + Up 



S{uj — k + q cos 6), 



(5.87) 



and a regular part that remains finite near the fermion mass shell 

ne^Tk f So r+1 



7rcg('^>k) 



. fA^ f '^[^^(g„,^)_cos20/3T(<7o,g) 
J {2Try J_ qo 



+ VTiqo,q)] S{iLJ - k + qcos6 - qo). 



Using the delta function S{qcos9 — qo) to perform the angular integrations, one obtains 

LOp 



-r\u,k) = aT In 



-fcl 



I] +0{uj-k), 



where 



I = 



dqq 

J-q qo 



ql 



f^Liqo, q) ^ f3T{qo,q) + i^t(<7o, q) 



(5.88) 



(5.89) 



(5.90) 



The above double integral has been computed analytically in Ref. [116| with the result X = In 3/2. 

We are now in position to find the secular terms in ^p^^^ (k, t) that emerge in the intermediate 
asymptotic regime. The imaginary part of the secular terms in e^V'^^'^H'^i and e^tp''^'''\h,t) can 
be combined into the imaginary part of the secular term 



luiSkit) = TiaSk t 

for the positive (upper sign) and negative (lower sign) energy spinors, where 

1 



a 5fe = - Tr 
2 



/i+(k)ReS(a',k) 



o—k 



(5.91) 



(5.92) 



There are no further secular terms arise from the higher order expansion around the fermion mass 



shell in Eq. (5.89). This (linear) imaginary secular term is thus identified with a perturbative shift 



94 



of the oscillation frequency of the fermion mean field [ p4| . It is a finite quantity determined by a 
dispersive integral of the spectral functions pi.2('^,k). 

The real secular terms are more involved. The finite contribution to 7''(a;, k) as — s- fc leads to 
a linear secular term, whereas for the logarithmically divergent contribution as a; ^ fc the following 
asymptotic result becomes useful ]8^ ]: 

r ^(l-cosyt)ln|y|*=°°7rt(l-7-lnt) + 0(<-i), (5.93) 

J-oo V 

where we have neglected terms that fall off at long times. Thus, one obtains the real part of the 
secular terms to be given by 



In 3 

ReSk{t) = ~aTt ( Incjpi + 7 - 1 + ^ 



Gathering the above results, at large times t ^ 1/cjp the perturbative solution reads 

^(k,i) - ^ + (t)] A, (k)w,(k)e-''=*+ [1 + 5^(^)1 i?,(k)i;,(-k)e''^-* } 



s=l,2 

nonsecular terms, 



(5.94) 



(5.95) 



with 



-at 



Incijpi + 7 — 1 ■ 



In 3 



T + i5, 



(5.96) 



Obviously, this perturbative solution breaks down on a time scale ~ [ann(l/Q;)]~^. To obtain 
a uniformly valid solution for large times, we now implement a resummation of the secular terms in 
the perturbative series via the dynamical renormalization group. As before, introducing (complex) 
renormalization constants for the amplitude 



A,(k) =Zfc(r)A(k,T), B,(k) =Zfc*(r)B,(k,T), ^^(t) = 1 + a Zfc(r) + ^(e^) 



(5.97) 



and choosing azkij) = —Skir) to cancel the secular divergence at time r, one finds to order 0{a) 
the dynamical renormalization group equations to be given by (after setting t = t) 



d dSkjt) 
dt dt 



A(k,i) ^ 0, 



d dSlit) 
dt 



dt 



B4k,t)^0. 



(5.98) 



Solving the above equations, we find that at asymptotic times {ujpt 3> 1) the fermion mean field 
relaxes as 



[As{K to) us{k) e-^^«* + Bs{k, To) vs{-k) e^^W*' 



s=l,2 



-aTt[\nujpt+0.126...] 



(5.99) 



where we have replaced 7 — 1 + ln3/2 = 0.126. . ., tq ^ is the time scale such that this 

intermediate asymptotic solution is valid, and tD(fc) ~ k + Sk is the position of the fermion pole 



(in hard momentum limit) shifted by one-loop corrections. Eq. (5.99) reveals a time scale for 
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the relaxation of the fermion mean field <rci ^ l/ann(l/a), which coincides with the time scale 
on which the perturbative solution given by Eq. ( 5.95| ) breaks down. This highlights clearly the 
nonperturbative nature of relaxation phenomena. 



Our result coincides with that found in Refs. | 116| , 117 1 via the Bloch-Nordsieck approximation 



and in scalar quantum electrodynamics p4| using the dynamical renormalization group method. 
Furthermore, it provides another important and relevant example of the reliability and consistency 
of this novel renormalization group method in study real-time nonequilibrium dynamics. 

5.3.2 Quantum kinetics of fermion quasiparticles 

We now study the quantum kinetic equation for the distribution function of hard fermion quasipar- 
ticles. There has recently been an intense activity to obtain a Boltzmann equation for quasiparticles 
in gauge theories motivated in part by the necessity to obtain a consistent description for 

baryogenesis in non-Abelian theories. Boltzmann equations with a diagrammatic interpretation were 
obtained in | pO| , |4^ ] in which a collision-type kernel describes the scattering of hard quasiparticles. In 
these approaches this collision kernel reveals the infrared divergences associated with the emission 
and absorption of magnetic photons (or gluons) and must be cutoff by introducing a magnetic mass 
Wmag ~ OiT to leading logarithmic accuracy [Q. 

In a derivation of quantum kinetic equations for charged quasiparticles in SQED |]8^ using 
the dynamical renormalization group method, it was understood that the origin of these infrared 
divergences is the implementation of Fermi's golden rule that assumes completed collisions and takes 
the infinite time limit in the collision kernels. The dynamical renormalization group approach leads 
to quantum kinetic equations in real time in terms of time- dependent scattering kernels without any 
infrared ambiguity. 

In this section we implement this program in fermionic QED to derive the quantum kinetic 
equation for hard fermion quasiparticles. There are several important features of our study that 



must be emphasized: (i) as discussed in Sec. 5.1, gauge invariance is automatically taken into 



account by working directly with gauge invariant operators, thus the operator that describes the 
number of fermion quasiparticles is gauge invariant, (ii) a kinetic description relies on a separation 
between the microscopic and the relaxation time scales, this is warranted in a strict perturbative 
regime and applies to hard fermion quasiparticles, (iii) the dynamical renormalization group leads to 
a quantum kinetic equation in real time without infrared divergences since time acts as an infrared 
cutoff. 

As in the case of the photon, we begin by expanding the Heisenberg fermion field in terms of 
creation and annihilation operators as (in momentum space) 



^l^Y. ^(k,i)u,(k) +dt(„k,tK.(-k)] , 
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^t(k,i) 



[bl{-Kt)ul{-k)+d,{k,t)vl{].}] , (5.100) 



where = vk^ + rr? and bs(k,t) [fe|(k, t)] is the annihilation (creation) operator that destroys 
(creates) a free fermion of momentum k and spin s at time t. Here, we have retained the fermion 
mass to avoid the subtleties associated with the normalization of massless spinors and the massless 
limit will be taken below. 

In the hard momentum limit fc ^ T 3> m, the number operator for fermion quasiparticles with 
momentum k is then given by 



Nf{k,t) = ^ 6t(k,i)5^(k,t) 



s=l,2 



= V^(-k,i)/i+(k)7V(k,i). (5.101) 
Taking time derivative of Nf (k, t) and using the Heisenberg equations of motion, we find 



lim^ ze J (-P, i') [7° ^o" (-q, *') - 7 ' A^; (-q, t')] 

X /i+(k)7° V'+(k, t)) + c.c. (5.102) 



We obtain to ©(e^) 



-niit) = e^y J dqn ^ dt" {Ui{to) cos[{k - p - qo){t" - to)][pL{qo,q) 

X /C+(k, q) +2pT{qo, q) /C+(k, q)] + J\f2{to) cos[(fc + p - qo){t" - to)] 
X [pi(go,g)/Cr(k,q) + 2pT(go,g)/C2-(k,q)]}, (5.103) 

where Af and JC denote respectively the following statistical and kinematic factors 

M(i) = [1 - niit)] nl{t) nsiqo) ~ n{it) [1 - n^(t)][l + ns(go)], 
M-2.it) = [1 - 4(<)][1 - n/(0] ns(<Zo) - nl{t) n{,{t) [1 + ns(go)], 

/C±(k,q) = l±k.p, /C2±(k,q)-lTk-p±^— fc^, (5.104) 

l-f(k-q) 

and PTiqojq) Bind PL^qoiq) are the HTL-resummed spectral functions for the transverse and longi- 



tudinal photons given by Eqs.(5.21) and (5.70), respectively. 

Consistent with the use of the HTL-resummed photon spectral functions, we work in the relax- 
ation time approximation in which only the initial distribution function for fermion quasiparticles of 
momentum k is slightly off-equilibrium such that Ti^{to) — np{k) + Sn^{to) while all other (fermion 
and photon) modes are in thermal equilibrium. Then upon integrating over t" , one obtains 

^^Sniit) = -Sniito) fy.nM^) ''''^%'l%''^\ (5.105) 
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where 



7^/(u;,k) 



7re 



(27r)3 



+ 00 



dqa [l + nBiqo) - npip)] |[pL(9o,g) ^^(k,q) 
2pT(go,9)/C^(k,q)] (5(w-p-go) + [pL{qo,q) JC^ (Kq) 
2pT(go,9) ^2"(k, q)] d{uj +P + qo)y 



Eq. (5.105) can be integrated directly to yield 

^+00 



Hit) = Snlito) 



1 



dujTZf{uj, k) 



1 — cos[(a; — k){t — to)] 



(5.106) 



(5.107) 



t:[uj — kY 

The time-dependent contribution above is now familiar from the previous discussions, potential 
secular terms will emerge at long times from the regions in which the resonant denominator vanishes. 
This is the region near the fermion mass shell a; « fc, where 7?./(cl;, k) is dominated by the regions of 
small q and go, which physically corresponds to emission and absorption of soft photons. As before 
for soft photons with q, qo <^ T, we can replace 



l + nsiqa) -npip) ^T/qo, p~k-qcos9, ICf 
thus writing 7?./(w,k) at w « fc as 



1 — cos 



(5.108) 



7^/(c.,k) = 27Te'T J ^ J' ^[^i-cos'e)(3T{qo,q)+/3L{qo,q)] 



X d{u! — k + q cos 6 — qo). 



(5.109) 



Note that the double integral in Eq. (5.10£) is exactly the same as that in Eq. (5.9C). Thus 7?./(w, k) 
features an infrared divergence near the fermion mass shell as shown in the previous subsection. 
Following the analysis carried out in the preceding subsection we obtain 



nf{L^,k) 



-2aT 



In- 



In 3 
~2~ 



+ O[{L0~kf 



(5.110) 



Substituting Eq. ( |5.110D into Eq. ( ^.107D , we find the number of hard fermion quasiparticles at 
intermediate asymptotic times i — io 3> 1 /ojp to be given by 

5nl{t) = (5n((io){l - 2aT(i - to)[\nujp{t - to) + 1 



0.126. 



nonsecular terms. 



(5.111) 



As in the case of the fermion mean field relaxation [see Eq. (5.95)], the perturbative solution contains 
a secular term of the form tint. Obviously, the secular term will invalidate the perturbative solution 
on time scales tj-ci ~ l/[2Q;Tln(l/a)] (which will be identified as the relaxation time scale below). In 
the intermediate asymptotic regime 1/fc <C t — to trci, the perturbative expansion can be improved 
by absorbing the contribution of the secular term on a time scale r into a "renormalization" of the 
distribution function. Hence we apply the dynamical renormalization group method through a 
renormalization of the distribution function much in the same manner as the renormalization of the 
amplitude in the mean field discussed above. 



5nl{to) = Zfc(T, to) 6nl{T), Zk{T, to) ^ 1 + a Zkir, to) + 0{a^). 



(5.112) 
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The independence of the solution on the time scale t leads to the dynamical renormalization group 
equation which to lowest order in a is given by 

4- + 2aT[l + Inwpfr - to) + 0.126 . . .] 1 SnlM = 0, (5.113) 
dr J 

which is the quantum kinetic equation in the relaxation time approximation, but with a time- 
dependent relaxation rate. 

Solving Eq. ( ^.113 ) and choosing r to coincide with t, we obtain the evolution of the fermion 
distribution function at asymptotic times t — to S> 1/fc to be given by 

Sn[{t) ~Snl{to) exp{-2ar(t-to)[lnwp(t- to) +0.126...]}. (5.114) 



Comparing Eq. ( 5.114 ) with Eq. ( 5.99 ) we find that the anomalous exponent that describes the 
relaxation of the fermion distribution function in the linear approximation is twice that for the 
linear relaxation of the mean field. A similar relation is obtained between the damping rate for 
the single quasiparticle relaxation and the relaxation rate of the distribution function in the case of 
time-independent rates and true exponential relaxation . The dynamical renormalization group 
reveals this to be a generic feature even with time-dependent rates. 



5.4 Conclusions 

The goals of this chapter are the study of nonequilibrium dynamics in QED plasmas at high tem- 
perature directly in real time. The focus is a systematic study of relaxation of mean fields as well 
as the distribution functions for photons and fermion quasiparticles. In particular the application 
of the dynamical renormalization group method to study anomalous relaxation as a consequence of 
the exchange of soft photons. 

To begin with, we have cast our study solely in terms of gauge invariant quantities, this can be 
done in an Abelian gauge field theory in a straightforward manner, avoiding potential ambiguities 
associated with gauge invariance. The relaxation of photon mean fields revealed important features: 
For soft momentum k < eT mean fields that are prepared by adiabatically switching-on an external 
source there is exponential relaxation towards the oscillatory behavior dominated by the transverse 
photon pole ^^(fc). The source that induces the mean field in this case has a Fourier transform 
that is singular at zero frequency and excites the resonance in the Landau damping cut near zero 
frequency for the soft photon mean field. Sources that have a regular Fourier transform would not 
lead to the exponential relaxation. For semihard momentum eT <ti k <ti T in principle both the HTL 
and the perturbative approximations are valid, however the spectral function for photons becomes 
sharply peaked at the edge of the Landau damping continuum consistent with the fact that the 
photon pole becomes perturbatively close to the Landau damping cut. This enhancement of the 
spectral function near the bare photon mass shell results in the breakdown of perturbation theory 
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at large times /ct 3> 1. The dynamical renormalization group provides a consistent resummation of 
the photon self-energy in real time and leads to anomalous power law relaxation of the mean field. 
Clearly, higher order terms beyond the HTL approximation will include coUisional contributions 
leading perhaps to exponential relaxation. We then expect a crossover between the anomalous power 
law obtained in lowest order and the exponential relaxation from higher order collision processes, 
the crossover time scale will depend on the details of the different contributions and requires a study 
beyond that presented in this thesis. 

The dynamical renormalization group provides a consistent and systematic framework to obtain 
quantum kinetic equations directly in real time from the underlying microscopic field theory |^^ . 
This method allows to extract information that is not available in the usual kinetic description in 
terms of time-independent collision kernels obtained under the assumption of completed collisions 
which only include on-shell (energy-conserving) processes. The dynamical renormalization group 
approach to quantum kinetics consistently includes off-shell (energy-nonconserving) processes and 
accounts for time- dependent collisional kernels. This is important and potentially phenomenolog- 
ically relevant in the case of the quark-gluon plasma which has a finite lifetime. In the case that 
fermions are in thermal equilibrium, we have implemented this approach to obtain the lowest order 
quantum kinetic equation for the distribution function of semihard photons in the HTL approxi- 
mation. This equation features time-dependent rates and we established that detailed balance, a 
consequence of fermions being in thermal equilibrium, emerges on microscopic time scales. The 
linearization of the kinetic equation describes relaxation towards equilibrium with an anomalous 
exponent, twice as large as that of the photon mean field in the semihard case. 

The relaxation of the fermion mean field for hard momentum is studied with the dynamical 
renormalization group method. This method reveals clearly in real time the emergence of the 
relaxation time scale t^ci ^ l/[Q;Tln(l/a)] at which perturbation theory breaks down. We find an 
anomalous exponential relaxation at large times of the form ^ exp[— Q!rt(lna;pt + 0.126 . . .)] which 



confirms the results of Refs. 1 116, 117], where the Bloch-Nordsicck approximation was used. We then 
obtain the quantum kinetic equation for the distribution function of hard fermion quasiparticles in 
the relaxation time approximation. The collisional kernel is time- dependent and infrared finite as the 
inverse of the time acts like an infrared cutoff. The linearized kinetic equation describes approach 
to equilibrium with an anomalous exponential relaxation, which is twice that of the fermion mean 
field. An important payoff of this approach to quantum kinetics is that it bypasses the assumption 
of completed collisions which leads to collisional kernels obtained by Fermi's golden rule and only 
describes on-shell (energy-conserving) processes as in the usual kinetic approach, which in the case 
under consideration leads to infrared divergent collisional kernels [Q . 

Perhaps the most phenomenologically pressing aspects that requires further and deeper study is 
the photon production by off-shell (energy-nonconserving) processes in a thermalized quark-gluon 
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plasma with finite lifetime. This is important in view of an assessment of experimental electromag- 
netic signatures of the QGP which is expected soon to be produced in ultrarelativistic heavy ion 
experiments at RHIC. The first step in the study of nonequilibrium aspects of photon production 
from a thermahzed QGP must be to include the finite QGP hfetime and to obtain the momentum 
distribution of the photons produced as a function of the temperature and lifetime of the QGP 
directly in real time via a quantum kinetic description. Furthermore, to compare with experimental 
data, an important next step is to include the hydrodynamical expansion of the QGP by coupling 
the quantum kinetic equation for photons to the hydrodynamic equations for the space-time evo- 
lution of the locally thermahzed QGP and to obtain the momentum spectrum of the photon yield 
as a function of the initial temperature of the QGP. These aspects will be studied in detail in the 
succeeding chapter. 



Chapter 6 

Direct Photon Production from 
the Quark-Gluon Plasma 



6.1 Introduction 

The first observation of direct photon production in ultrarelativistic heavy ion collisions has been 
reported recently by the CERN WA98 collaboration in 208p|3_^208p|-, collisions at y/sWN = 158 



GeV at the Super Proton Synchrotron (SPS) |119|. Most interestingly, a clear excess of direct 



photons above the background photons predicted from hadronic decays is observed in the range of 
transverse momentum pT > 1.5 GeV/c in central collisions. As compared to proton-induced results 
at similar incident energy, the transverse momentum distribution of direct photons shows excess 
direct photon production in central collisions beyond that expected from proton-induced reactions. 
These findings indicate not only the experimental feasibility of using direct photons as a signature 
of the long-sought quark-gluon plasma (QGP) but also a deeper conceptual understanding of 
direct photon production in ultrarelativistic heavy ion collisions. 

Unlike many other new phases of matter created in the laboratory, the formation and evolution 
of the QGP in ultrarelativistic heavy ion collisions is inherently a nonequilibrium phenomena. Cur- 
rently, it is theoretically accepted that hard parton scatterings thermalize quarks and gluons on a 
time scale of about 1 fm/c, after which the plasma undergoes hydrodynamic expansion and cools 
adiabatically down to the quark-hadron phase transition. If the transition is first order, quarks, 
gluons, and hadrons coexist in a mixed phase, which after hadronization evolves until freeze-out. 
Estimates based on energy deposited in the central collision region at the BNL Relativistic Heavy Ion 
Collider (RHIC) energies y^s^vw ~ 200 GeV suggest that the lifetime of the deconfined QGP phase 
is of order 10 fm/c with an overall freeze-out time of order 100 fm/c. Different types of signatures 
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are proposed for each different phase As mentioned in Chap. ^, direct photons and dileptons 
emitted by the QGP (electromagnetic probes) are free of final state interactions and can provide 
a clean signature of the early stages of a thermalized plasma of quarks and gluons. Therefore a 
substantial effort has been devoted to a theoretical assessment of the spectra of direct photons and 



dileptons emitted from the QGP [|llO|-yj, |120 1 . 

Conventionally, the theoretical framework for studying direct photon production from a ther- 
malized QGP begins with an equilibrium calculation of the photon emission rate at finite tem- 



perature I lie, 112-114, 12C]. This static rate is then combined with the hydrodynamical descrip- 
tion of QGP to obtain the total yield of direct photons produced during the evolution of the 



QGP |1JX, n^, 114, 121-123|. Whereas this approach is physically intuitive and widely used in 



the literature, it is solely based on semiclassical Boltzmann kinetics, which completely neglects tran- 
sient effects arising from the finite QGP lifetime. The latter is of particular importance in the study 
of experimental signatures of the QGP because the formation and evolution of the QGP at currently 
accessible energy scales is by itself a transient phenomenon. 

In this chapter we focus in particular on experimental signatures associated with processes that 
would be forbidden by energy conservation in a QGP of infinite lifetime. We first argue that the 
finite lifetime of a transient QGP raises a conceptual inconsistency in the calculation of direct 
photon production via an equilibrium rate. We then introduce a real-time kinetic description which 
naturally accounts for the finite- lifetime and nonequilibrium aspects of the QGP and includes energy- 
nonconserving effects. This real-time kinetic approach can be consistently incorporated with the 



widely used Bjorken's hydrodynamics |124, 125| | to obtain the direct photon yield for an expanding 
QGP that is expected to be created in central collisions at RHIC and Large Hadron Collider (LHC) 
energies. It is not our goal to assess photon production from the hadronic phase but to compare the 
real-time kinetic predictions for the QGP phase to those obtained from the equilibrium calculations. 

Recently, Alam et al. |l23| have evaluated the photon yield for Pb+Pb collisions at SPS energies 
from the initial state up to freeze-out by using the equilibrium photon production rate for the QGP 
and hadronic phases within the framework of (3-|-l)-dimensional hydrodynamic expansion. These 
authors also provide an estimate of the photon yield for Au+Au collisions at RHIC energies. Their 
results show that the photon yield from the QGP phase and that from the hadronic phase are of 
comparable order. Hence, our goal here is to compare the nonequilibrium yield from the QGP to 
that obtained from the equilibrium formulation with the QGP phase of the plasma. As we will see 
below, a substantial enhancement arising from nonequilibrium effects associated with the transient 
QGP lifetime indicates that the nonequilibrium yield from the QGP phase will stand out over the 
equilibrium yield from the hadronic phase. 



This chapter is organized as follows. In Sec. 6.2, we first review the usual S'-matrix approach to 
direct photon production from a QGP. We then argue that this approach in obtaining an equilibrium 
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rate has shortcomings and is conceptually and physically incompatible with photon production from 
an expanding QGP with a finite lifetime. In Sec. |6.3| , Bjorken's hydrodynamical model combined 
with the S'-matrix rate calculation of photon production from the QGP is briefly summarized and 



the incompatibility of these two approaches is highlighted. In Sec. 6.4, we introduce the real-time 
kinetic approach to photon production from a longitudinally expanding QGP which is consistently 
combined with Bjorken's hydrodynamics. This section contains our main results and we compare 



these to those obtained from the usual approach. Our conclusions are presented in Sec. p. 5 



6.2 S'-matrix Approach and its Shortcomings 



Production of direct photons from a QGP in thermal equilibrium has been studied extensively [109 



110, 112, 113, 120| | because of its relevance as a clean probe of the QGP. We begin with highlighting 
the main assumptions that are explicit in all previous calculations of photon production from a 
thermalized quark gluon plasma and that are explicitly displayed by the derivation above. First, 
the initial state at ti (which in the usual calculation is taken to — oo) is taken to be a thermal 
equilibrium ensemble of quarks and gluons but the vacuum state for the physical transverse photons. 
Furthermore, the buildup of the photon population is neglected under the assumption that the 
mean free path of the photons is larger than the size of the plasma and the photons escape without 
rescattering. This assumption thus neglects the prompt photons produced during the pre-equilibrium 
stage. Indeed, Srivastava and Geiger [|l26| have studied direct photons from a pre-equilibrium 
stage via a parton cascade model that includes pQCD parton cross sections and electromagnetic 
branching processes. The usual computation of the prompt photon yield during the stage of a 
thermalized QGP assumes that these photons have left the system and the computation is therefore 
carried out to lowest order in a with an initial photon vacuum state. Obviously keeping the pre- 
equilibrium photon population results in higher order corrections in a. In taking the final time 
tf to infinity in the S'-matrix element the assumption is that the thermalized state is stationary, 
while in neglecting the buildup of the population the assumption is that the photons leave the 
system without rescattering and the photon population never builds up. These assumptions lead to 
considering photon production only the lowest order in a, since the buildup of the photon population 
will necessarily imply higher order corrections. Although these main assumptions are seldom spelled 
out in detail, they underlie all previous calculations of the photon production from a thermalized 
quark-gluon plasma. 

We now review some important aspects of the usual S'-matrix calculation so as to highlight 
its shortcomings and establish contact with the real-time kinetic approach to photon production 



introduced in Sec. 6.4. It is convenient to write the total Hamiltonian in the form 



H^Ho + Hint , Ho = Hqcd + H^ , i^int ^ e J d^x J^A^ , (6.1) 
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where Hqqd is the fuU QCD Hamiltonian, is the free photon Hamiltonian, and iJjnt is the 
interaction Hamiltonian between quarks and photons with the quark electromagnetic current, 
the photon field, and e the electromagnetic coupling constant. 

Consider that at some initial time ti the state \i) is an eigenstate of Hq with no photons. The 
transition amplitude at time tf to a final state |/, 7a(p)) = |/) ^ |7a(p)), again an eigenstate of Hq 
but with one photon of momentum p and polarization A, is up to an overall phase given by 

S{tf,U)^{f,j^{p)\U{tf,U)\t), (6.2) 

where U{tf,ti) is the time evolution operator in the interaction representation 



Uitf,U) = Texp 



~ 1-ie J^^ dt J cfxJ^{K,t)Af,j{x,t)+0{e'^), (6.3) 

where the subscript / stands for the interaction representation in terms of Hq. In the above ex- 
pression we have approximated U{tf,ti) to first order in e, since we are interested in obtaining 
the probability of photon production to lowest order in the electromagnetic interaction. The usual 
5-matrix element for the transition is obtained from the transition amplitude S'(i/, ti) above in the 
limits ti — > —oo and — > oo 



Sfi = 5(+oo,— oo) 

dte'P"''^ e^{f\J^'{x)\i), (6.4) 



'2E . 

where E = |p| and = (-E, p) are the energy and four-momentum of the photon, respectively, 
and is its polarization four-vector. Since the states \i) and |/) are eigenstatcs of the full QCD 
Hamiltonian i/qcD, the above 5- matrix element is obtained to lowest order in the electromagnetic 
interaction, but in principle to all orders in the strong interaction. 

The rate of photon production per unit volume from a QGP in thermal equilibrium at tempera- 
ture T is obtained by squaring the S-matrix element, summing over the final states, and averaging 
over the initial states with the thermal weig ht e"''^VZ(/3), where (i = 1/T, E., is the eigenvalue 
of Hq corresponding to the eigenstate and Z{0) ~ J2i^~^^' is the partition function. Using 
the resolution of identity 1 = J2f I /)(/!> the sum of final states leads to the electromagnetic cur- 
rent correlation function. Upon using the translational invariance of this correlation function, the 
two space-time integrals lead to energy-momentum conservation multiplied by the space-time vol- 
ume fl = V{tf — ti) from the product of Dirac delta functions. One can recognize that the limit 
tf — ti oo is the usual interpretation of 27r(5(0) in the square of the energy conservation delta 
function. 



These steps lead to the following result for the photon production rate in thermal equilibrium [110 
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< d^P 



2E{2t:) 



(6.5) 



where W^^ is the Fourier transform of the thermal expectation value of the current correlation 



function defined by 

W<, = j d^xe^P-- {J^{Q)Ux))p. (6.6) 

In the expression above ( • )/3 denotes the thermal expectation value. To lowest order in but to 
all orders in the strong interactions, is related to the retarded photon self-energy 11^^^ by [114| 

2 



Imll^ 



Thus, one obtains the (Lorentz) invariant photon production rate 
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(6.7) 



(6.8) 



'd^pd^x (27r)3 ""^ e-f^/^ - 1 ■ 
Using the hard thermal loop (HTL) resummed effective perturbation theory developed by Braaten 



and Pisarski |62|, Kapusta et al. |112{ and Baier et al. |113| showed that at one-loop order (in effective 
perturbation theory) the processes that contribute to photon production are the giuon-to-photon 
Conipton scattering off (anti)quark q{q)g — s- q{q)"f and quark-antiquark annihilation to photon and 
gluon qq 57. The resultant rate of energetic (E 3> T) photon emission for two light quark flavors 
{u and d quarks) is given by |ll2| 



E 



dN 



d^p d'^x 



5 aUs 
92^ 



T^e-^/^ln 



0.23£;\ 



(6.9) 



one-loop 

where a is the fine-structure constant and — qI/'^t^ with being the strong coupling constant. 
In a recent development Aurenchc et al. |120(| have found that the two-loop contributions to the pho- 
ton production rate arising from (anti)quark bremsstrahlung qq{g) qq{g)l and quark-antiquark 
annihilation with scattering qqq{g) q{g)j are of the same order as those evaluated at one loop. 



The two-loop contributions to the photon production rate read |120 



E- 



dN 



d^p d^x 
1.11 and Jl 



two-loop 
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E 
3f 



(6.10) 



where Jt ~ 1-11 and Jl ~ —1.07 for two light quark flavors [127|. Most importantly, they showed 
that the two-loop contributions completely dominate the photon emission rate at high photon ener- 



gies [120|. We emphasize that the thermal photon production rates (|6.9[) and (6.1C) [or the general 



result (3.S)] has two noteworthy features: (i) the thermal rate is a static, time-independent quantity 
as a result of the equilibrium calculation and (ii) emission of high energy photons is exponentially 
suppressed by the Boltzmann factor e~^^^ . 
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We have reproduced the steps leadmg to Eq. (6.8), which is the expression for the photon pro- 
duction rate used in aU calculations available in the literature, to highlight the important steps in its 
derivation in order to compare and contrast to the real-time analysis discussed below. In particular, 
in the above derivation the transition amplitude is obtained via the time evolution operator U(tf, ti) 
evolved from the initial time ti up to a time tf. In the usual S'-matrix calculation, one takes the 
limits ti — !■ — oo and tf oo. Taking the infinite time limit and squaring the transition amplitude 
lead to energy conservation and a time-independent transition rate, which in turn imply the validity 
of semiclassical Boltzmann kinetics. 

The main reason that we delve on the main assumptions and specific steps of the usual com- 
putation is to emphasize that there is a conceptual limitation of this approach when applied to an 
expanding QGP of finite lifetime. The current theoretical understanding suggests that a thermalized 
QGP results from a pre-equilibrium partonic stage on a time scale of order 1 fm/c after the collision, 
hence for consistency one must choose ~ 1 fm/c. Furthermore, within the framework of hydrody- 
namic expansion, studied in detail below, the QGP expands and cools during a time scale of about 
10 fm/c. Hence for consistency to study photons produced by a quark gluon plasma in local thermal 
equilibrium one must set tf ^ 10 fm/c. Hydrodynamic evolution is an initial value problem indeed 
the state of the system is specified at an initial (proper) time surface (to be local thermodynamic 
equilibrium at a given initial temperature) and the hydrodynamic equations are evolved in time to 
either the hadronization or freeze out surfaces. The calculation based on the S'-matrix approach 
takes the time interval to infinity, extracts a time-independent rate and inputs this rate, assumed to 
be valid for every cell in the comoving fluid, in the hydrodynamic evolution during a finite lifetime. 

As stated in the Introduction, however, the QGP produced in ultrarelativistic heavy ion collisions 
is intrinsically a transient and nonequilibrium state. It is therefore of phenomenological importance 
to study nonequilibrium effects on direct photon production from an expanding QGP with a finite 
lifetime with the goal of establishing potential experimental signatures. 

Our main observation is that the usual computations based on S'-matrix theory extract a time 
independent rate after taking the infinite time interval, which is then used in a calculation of the 
photon yield during a finite time hydrodynamic evolution. While we do not question the general 
validity of the results obtained via the S-matrix approach, we here focus on the signature of processes 
available during the finite lifetime of the QGP and that would be forbidden by energy conservation 
in the infinite time limit. 



6.3 Bjorken's Hydrodynamical Model 

The current understanding of the QGP formation, equilibration, and subsequent evolution through 
the quark-hadron (and chiral) phase transitions is summarized as follows. A pre-equilibrium stage 
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dominated by parton-parton interactions and strong colored fields which gives rise to quark and 
gluon production on time scales < 1 fm/c [^|^- The produced quarks and gluons thermalize via 
elastic collisions on time scales ^ 1 fm/c. Hydrodynamics is probably the most frequently used 
model to describe the evolution of the next stage when quarks and gluons are in local thermal 



equilibrium (although perhaps not in chemical equilibrium) ]124 , 125|. The hydrodynamical picture 
assumes local thermal equilibrium (LTE), a fluid form of the energy- momentum tensor and the 
existence of an equation of state for the QGP. The subsequent evolution of the QGP is uniquely 
determined by the hydrodynamical equations, which are formulated as an initial value problem with 
the initial conditions specified at the moment when the QGP reaches local thermal equilibrium, i.e., 
at an initial time ti ^ 1 fm/c. The (adiabatic) expansion and cooling of the QGP is then followed to 
the transition temperature at which the equation of state is matched to that describing the mixed 



and hadronic phases 1 114, 122, 123 1. 

In order to highlight the conceptual limitation of the S'-matrix calculation for direct photon 
production for a transient QGP, we now review the essential features of the hydrodynamic descrip- 



tion 1 124, 125|. For computational simplicity we work within Bjorken's hydrodynamical model of 



a longitudinally expanding QGP |124]. The main assumption in Bjorken's model is longitudinal 
Lorentz boost invariance in the central rapidity region of the QGP. This is motivated by the obser- 
vation that the particle spectra for the secondaries produced in p+N and N+N collisions exhibit a 
central plateau in the rapidity space near midrapidity. For a longitudinally expanding QGP, it is 
convenient to introduce the proper time r and space-time rapidity rj variables defined by 

1 , t + z 



r = ^¥^\ ^= Ini^, (6.11) 
2 t — z 

where t and z, respectively, are the time and spatial coordinate along the collision axis in the center- 
of-mass (CM) frame. The transverse spatial coordinates will be denoted as xy, hence the space-time 
integration measure expressed in terms of r, 77 and xt is given by d^x — t dr drj cPxt- Invariance 
under (local) longitudinal Lorentz boost implies that thermodynamic quantities are functions of r 
only and do not depend on r]. 



In Bjorken's scenario |124, 125| | the QGP reaches local thermal equilibrium at a temperature 
Ti at a proper time of order ^ Ifm/c after the maximum overlap of the colliding nuclei. The 
initial conditions for hydrodynamical equations are therefore specified on a hypersurface of constant 
proper time r^. The equation of state for the locally thermalized QGP is taken to be that of 
the ultrarelativistic perfect radiation fluid (corresponding to massless quarks and gluons). The 
longitudinal expansion is described by the scaling ansatz — z/t, where is the collective fluid 
velocity of the hydrodynamical flow and describes free streaming of the fluid. Hence the space-time 
rapidity equals to the fluid rapidity. In terms of r and 77 the scaling ansatz implies that the four- 
velocity of a given fluid cell in the CM frame is given by — (cosh?7, 0, 0, sinhjy) with u^u^ — 1. 
The conservation of total entropy leads to adiabatic expansion and cooling of the QGP according to 
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1/3 



(6.12) 



the cooling law |T24[|25| 

T(r)=7^.(^; 

Hence, the QGP phase ends at a proper time t/ — Ti{Ti/Tc)^ , where Tc ^ 160 MeV is the quark- 
hadron transition temperature. At RHIC energies the initial thermalization temperature is estimated 
to be Ti ^ 200 — 300 MeV, which entails that the lifetime of the QGP phase is of order < lOfm/c. 
Within a hydrodynamical model the usual S'-matrix calculation of direct photon production from 



a expanding QGP proceeds as follows [112, 114, 121 - 123 1. 

1. First the rate of direct photon production is calculated within the S'-matrix framework de- 



scribed in the previous section, leading to Eq. (x8) for the invariant rate. This expression for 
the rate describes the photon production rate in the local rest (LR) frame of a fluid cell in 
which the temperature is a function of the proper time of the fluid cell. The rate in the CM 
frame is obtained by a local Lorcntz boost E P^u^ and the replacement T{t) T{t): 



dN 



(PpT dy T dr drj d^XT 
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dN 
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d'^p d^x 



(6.13) 



LR 



where pt and y, respectively, are the transverse momentum and rapidity of the photon. 

2. The direct photon yield is now obtained by integrating the rate over the space-time history of 
the QGP, from the initial hypersurface of constant proper time to the final hypersurface of 
constant proper time Tf at which the phase transition occurs. This leads to the following form 
of the total direct photon yield in the CM frame for central collisions: 



dN 
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(6.14) 



LR 



■-»?co„ d^pd^x 

where Ra is the radius of the nuclei and — r^ccn < t] < ?7ccn denotes the central rapidity region 
in which Bjorken's hydrodynamical description holds. The fact that the S'-matrix calculation 
for the rate results in a time-independent rate determines that the only dependence of the rate 
in the LR frame on the proper time is through the temperature which is completely determined 
by the hydrodynamic expansion. 

It is at this stage that the conceptual incompatibility between the S'-matrix calculation of the 
photon production rate and its use in the evaluation of the total photon yield from an expanding 
QGP of finite lifetime becomes manifest. The hydrodynamic evolution is treated as an initial 
value problem with a distribution of quarks and gluons in local thermal equilibrium on the initial 
hypersurface of constant proper time ~ 1 fm/c. The subsequent evolution determines that the 
QGP is a transient state with a lifetime of order < 10 fm/c. The direct photon yield is obtained 
by integrating the rate over this finite lifetime. The S-matrix calculation of the rate for a QGP in 
thermal equilibrium, on the other hand, implicitly assumes that Ti — > —oo and — > oo as discussed 
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above in detail. Therefore while the rate has been calculated by taking the time interval to infinity, 
it is integrated during a finite time interval to obtain the total yield. 

The question that we now address, which is the focus of this chapter, is the following: Is this 
conceptual incompatibility of physical relevance and if so what are the experimental observables? In 
order to answer this question and to assess the potential experimental signatures from nonequilibrium 
effects, we must depart from the S'-matrix formulation and provide a real-time calculation of the 
direct photon production rate based on nonequilibrium quantum field theory. 

6.4 Real-Time Kinetic Approach 

Having highlighted the conceptual shortcomings of the usual S'-matrix approach to photon produc- 
tion from the QGP, in this section we provide an alternative approach which bypasses the above 
mentioned conceptual shortcomings. This approach, based on the dynamical renormalization group 
method to quantum kinetics discussed in details in previous chapters, treats the production of pho- 
tons as an initial value problem directly in real time without invoking Fermi's golden rule. Compared 
to the usual approach, this real-time kinetic approach has the following noteworthy advantages: (i) 
It is capable of capturing energy-nonconserving effects arising from the finite lifetime of the plasma, 
as completed collisions are not assumed a priori, (ii) Since both the real-time kinetic approach and 
hydrodynamics are formulated as initial value problems, they can be incorporated consistently on 
the same footing. This last point proves very important in the hydrodynamic description of photon 
production. 

6.4.1 Nonexpanding QGP 

We begin our discussion with the calculation of the invariant photon production rate from a nonex- 
panding QGP. This calculation is relevant because the result is interpreted as the invariant photon 
production rate in the local rest frame of a fluid cell. The corresponding rate in the CM frame is 
obtained simply by a local Lorentz boost and the direct photon yield is obtained by integrating the 
rate over the space-time history of the QGP as explained in the previous section. 

Because of the Abelian nature of the electromagnetic interaction, we will work in a gauge invari- 
ant formulation in which physical observables (in the electromagnetic sector) are manifestly gauge 
invariant and the physical photon field is transverse . 

The real-time kinetic approach begins with the initially prepared density matrix and the time 
evolution of the expectation value of the photon number operator (see Chap.^for details). Consistent 
with the hydrodynamical initial value problem, we consider that at the initial time ti quarks and 
gluons are thermalized such that the initial state of the QGP is described by a thermal density 
matrix at a given initial temperature Ti. Furthermore, we consider that photons are not present at 



no 



the initial time. Although this last assumption ean be relaxed allowing an initial photon distribution, 
the usual approach is to assume that the photons produced during the pre- equilibrium stage had left 
the plasma without building up a population. Therefore the initial density matrix p is of the form 

P = PQCD® |0^)(0^|, PQCD = e-^«^°/^'. (6.15) 

We remark that the assumption that the initial density matrix is that of a thermalized system 
(after the strong interactions thermalize quarks and gluons on a time scale ^ 1 fm/c) underlies 
the program that studies the equilibrium properties of the QGP. This is our only assumption, i.e., 
that of a thermalized QGP at an initial time scale ti w 1 fm/c and is consistent with the general 
assumptions behind the equilibrium program. This assumption is elevated to that of LTE, again 
consistent with the hydrodynamical description of an expanding QGP. 

As usual we will focus on hard photon production, since hard photons suffer less from the 



background contamination and has been measured in recent experiments at SPS |119|. By assuming 
photons escape directly from the QGP without further interaction, one can treat them as free 
asymptotic particles. Consequently, the number operator for photon can be obtained from the 
Heisenberg photon field (and the conjugate momentum) in the same manner as we did for the hot 
QED plasma [see Eqs. (5.43) and ( ^.44 )]. The number of photons per unit phase space volume at 



time t is then given by Eq. ( 5.45| ), where the expectation value is now taken with the initial density 



matrix p specified above in Eq. ( |6.15 ). The invariant photon production rate E dN (t) / (Pp d'^x is 



obtained by using the Heisenberg equations of motion and can be written in the CTP formalism as 

3 

I (^(^7(-k,07-Aj(p,O^7(q,t))+c.c., (6.16) 

where k = p + q. In above expression, Nf is the number of quark flavors, e/ is the quark charge 
in units of the electromagnetic coupling constant e. At is the transverse component of the photon 
field, ■0/ is the (Abelian) gauge invariant quark fields (with color index suppressed). 

We assume the weak coupling limit a ^ <C 1. Whereas the first limit is justified and is 
essential for the interpretation of electromagnetic signatures as clean probes of the QGP, the second 
limit can only be justified for very high temperatures, and its validity in the regime of interest can 
only be assumed so as to lead to a controlled perturbative expansion. As mentioned above since 
there are no photons initially and those produced escape from the plasma without building up their 
population, therefore the QGP is effectively treated as the vacuum for photons. Consequently, the 
nonequilibrium expectation values on the right-hand side of Eq. ( |6.16| ) are computed perturbatively 
to order a and in principle to all orders in by using nonequilibrium Feynman rules and real-time 
propagators (but with vacuum photon propagators). The corresponding Feynman diagrams are 



depicted in Fig. 6.1 



We focus here on the lowest order (one-loop) contribution, which is missed by all the previous 
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Figure 6.1: The Feynman diagrams that contribute to the invariant photon production rate from the 
QGP. Figure (a) shows the contribution to order a and to all orders in with the blob denoting the 
full quark electromagnetic current correlation function to all orders in a^. Figure (b) is the lowest 
order contribution, of order a. 



investigations using the S-matrix approach implicitly assuming an infinite QGP lifetime. To lowest 
order in perturbation theory, neglecting the vacuum contribution one obtains the invariant photon 
production rate for a QGP of two light quark flavors {u and d) to be given by 



where 
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(6.17) 



(27r)= 



2[1 - (p • k)(p • q)]n(g)[l - n{k)]S{cj + k - q) 



- [1 + (p • k)(p • q)] n[q)n{k)5{L0 - k - q) 



(6.18) 



with n{q) — \/[e'^/'^' + 1] the quark distribution function at the initial time ti. We note that in 
obtaining the above expressions, we have used the free quark propagators which correspond to hard 
quark loop momentum (g > Ti). Soft quark lines {q <ti Tt) require using the HTL-resummed effective 
quark propagator |l62| , thereby leading to higher order corrections. Indeed, the one- loop diagram 
with soft quark loop momentum is part of the higher order contribution of order aas that has been 
calculated in Refs. [112, IIJ 



We note that apart from color and flavor factors TZ{lu) given by Eq. (6. IS) has the same structure 



as the flrst two terms in TZ'^{uj,k) given by Eq. ( 5.49 ). From the discussion there, one finds that 
the first delta function 6{uj + k — q) has support below the light cone {lo^ < E"^) corresponding to 
the Landau damping cut whereas the second delta function 5{lj ~ k — q) has support above the 
light cone corresponding to the usual two-particle cut. Furthermore, TZ{u]) has a clear physical 
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interpretation in terms of the following energy-nonconserving photon production processes: the first 
term describes (anti)quark bremsstrahhmg q{q) ~^ 9(9)7, and the second term describes quark- 
antiquark annihilation to photon qq 7. 

The dependence of TZ{lu) on the quark distribution at the initial time U is a consequence of 
the fact that in nonequilibrium quantum field theory the important ingredient is the initial density 
matrix, which consistent with hydrodynamics is taken to be thermal for quarks and gluons. Since 
the calculation is carried out consistently in perturbation theory, the quark (and gluon) propagators 
are in terms of the equilibrium free particle distribution functions. In the case of a nonexpanding 



and thermalized QGP, the integral over t" in Eq. ( 6.17 ) can be carried out directly leading to 

At this stage we can make contact with the S'-matrix calculation and highlight the importance 
of the finite-time, nonequilibrium analysis. If, as is implicit in the 5-niatrix calculation, the QGP is 
assumed in thermal equilibrium and with an infinite lifetime (entailing that ti —00 and tj — > 00), 
then we can take the infinite time limit — s- — 00 in the argument of the sine function in Eq. ( |6.19 ) 
and use the approximation 

Siu. E). (6.20) 

7r(aj — E) 

This is the assumption of completed collisions that is invoked in time-dependent perturbation the- 
ory leading to Fermi's golden rule and energy conservation. The delta function 5{lo — E) in a. 
manifestation of energy conservation for each completed collision. Under this assumption one finds a 
time-independent photon production rate proportional to TZ{E), provided that the latter is finite. In 
the present situation, however, the delta functions in TZ(lu) cannot be satisfied on the photon mass 
shell. Therefore, under the assumption of completed collisions the lowest order energy-nonconserving 
contribution to the photon production rate simply vanishes due to kinematics. Therefore this lowest 
order contribution is absent (by energy conservation) in the S'-matrix calculation, but is present at 
any finite time. 

The relevant question to ask is that how this finite-time contribution of order a compares to the 
higher order S-matrix contribution to the photon yield. 



For finite QGP lifetime the time-dependent rate given in Eq. (3.19) is finite and nonvanishing, 
thus leading to a nontrivial contribution to direct photon production. The photon yield (per unit 
volume) is obtained by integrating the rate over the lifetime of the QGP. Using Eq. ( 6.19| ), one 
obtains 

where tf ~ ti ^ lOfm/c is the lifetime of the QGP. 

Before proceeding further, we give an analytic estimate of the behavior of the photon yield in 
the HTL approximation. In this approximation the leading contribution of TZ(uj) is dominated by 
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Figure 6.2: The ratio of the nonequihbrium and equihbrium hard {E ^ T) photon yield plotted as 
a function of time for a = 1/137, as = 0.3 and T = 200 MeV (see text for details). 



the Landau damping cut, which corresponds to ofF-shcU (anti)quark bremsstrahlung, and is found 
to be given by 



20 TT^aT^ uj 



1 



3 12 E 

where use has been made of the KMS condition, valid for quarks in thermal equilibrium. 



(6.22) 



For i? <C T, 7lini^{uj) can be further simplified as [see Eq. ( |5.5l| )] 



E^T 20 TT^aT^ 
~ y 12£: 



1 - 



E^ 



e{E^-uj^). 



(6.23) 



The dominant contribution of the uo integral in Eq. (6.21) for E <^T arises from the region where 



the resonant denominator vanishes {to ~ E). Using Eq. (3.23), we obtain 



(6.24) 



where 7 — 0.577 ... is the Euler-Mascheroni constant. Figure |6.2| displays the ratio of the hard 



{E ~ T) photon yield extrapolated from Eq. (6.24) and that obtained from the equilibrium rate 
of Kapusta et al. given by Eq. ( |6.9[ ). We see clearly that the lowest order 0{a) nonequihbrium 
contribution to the photon yield is comparable at early times to the higher order 0{aas) equilibrium 
contribution. 



For E '3> T, as depicted in Fig. 3.3, 7?.htl('j^) is exponentially suppressed in the region T < uj < E. 
From this observation we emphasize the following important features of the nonequihbrium yield 



given by Eq. ( 6.21 ). (i) Because for ^ T the threshold contribution near the photon mass shell 
UJ K E is exponentially suppressed, the uj integral in Eq. (6.21) is now dominated by the interval 
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Figure 6.3: The function 7?.htl(w) plotted in the limit E ^ T. Here we take E = 2 GcV and 
T = 200 MeV. 



-E < u! < T, which corresponds to highly off-shell (anti)quark bremsstrahlung. (ii) As the integrand 



in Eq. ( 3.21 ) is positive-definite and the function 1 — cos[([x' — E)t] averages to 1 for large t in the 
region —E<lu<T, for fixed E ^ T the yield approaches a constant at large times, (iii) In contrast 
to those obtained from equilibrium rates, the yield for E ^ T is not suppressed by the Boltzmann 
factor e~^/"^. 



Figure 5.4 shows a comparison of the nonequilibrium and equilibrium contributions to the hard 
photon yield in the range of energy T < E < 5 GeV from a QGP of temperature T = 200 MeV 
and lifetime t = 10 fm/c. Whereas the two- loop contribution dominates the direct photon yield for 
smaller values of i?, a significant enhancement of the direct photon yield due to the nonequilibrium 
contribution is seen at i? > 2 GeV as a consequence of its power law falloff for E ^ T. Numerical 
evidence shows that the nonequilibrium yield EdN{t)/(Pp<fix falls off with a power law E^'^ with 
V w 2.14 for T ^ E < 5 GeV. In particular, the nonequilibrium contribution is larger than the 
equilibrium contributions by several orders of magnitude for E > 3 GeV. As the nonequilibrium 
contribution for fixed E ^ T approaches a constant at large times, the equilibrium contributions, 
which grow linearly in time, will eventually dominate the yield if the QGP has a very long lifetime. 
However, the linear growth in time of the equilibrium contributions has to compensate the Boltzmann 
suppression for E ^ T. Therefore we emphasize that for E ^ T the equilibrium contributions could 
dominate the yield only if the QGP lifetime is of order 10^ — 10^ fm/c or larger, which nevertheless 
is very unrealistic at RHIC energies. 

While these results in the nonexpanding case revealed the importance of the nonequilibrium and 



115 




1 2 3 4 5 

E (GeV) 



Figure 6.4: Comparison of various contributions to the direct photon yield from a QGP of temper- 
ature T = 200 MeV and hfetime i = 10 fm/c. The inset shows the figure on a log-log plot. 



finite-lifetime aspects, the most experimentally relevant case to study is that of an expanding QGP. 
This experimentally relevant case can only be investigated by a detailed numerical study. 



6.4.2 Longitudinally expanding QGP 



As a prelude to photon production from an expanding QGP, we first focus on the invariant photon 
production rate from each individual fluid cell of the QGP. Since the proper time equals to the local 
time in the local rest frame of any fluid cell, we can follow the same real-time nonequilibrium analysis 
presented in the proceeding subsection to calculate the nonequilibrium invariant photon production 
rate in the local rest frame of each fluid cell. We obtain 



E 



dN{T) 



lim 



LR 



Nj 

7E 



ee/ 



2(27r)3 J (27r) 



(V;7(-k,r)7-A+(p,T')^7(q,T)) 



(6.25) 



As before, the nonequilibrium expectation values on the right-hand side of Eq. ( 6.25 ) is computed 
perturbatively to order a and in principle to all orders in by using real-time Feynman rules and 
propagators. To lowest order in perturbation theory, the result for two light quark flavors {u and d) 
reads 

dN{T) 



E 



d^p d'^x 



LR 



(27r)3 



dr" 



+ 00 



dio 



n{uj) cos[{uj~E){t-t")], 



(6.26) 



where TZ{uj) is the same as that given in Eq. ( |6.18| ), but now with n{q) — XjY^I^^ + 1] being the quark 
distribution function at initial proper time Ti , at which the QGP reaches local thermal equilibrium. 
In principle in the case of an expanding QGP under consideration, the photon production rate 



given by Eq. (3.26) has to be supplemented by kinetic equations that describe the evolution of the 
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quark and gluon distribution functions so as to setup a closed set of coupled equations. However 
the assumption of the validity of (ideal) hydrodynamics entails that the quarks and gluons form 
a perfectly coupled fluid. This in turn implies that the mean free paths of the quarks and gluons 
are much shorter than the typical wavelengths and the relaxation time scales are much shorter 
than the typical time scales, i.e., the quark and gluon distribution functions adjust to local thermal 
equilibrium instantaneously. Thus the assumption of the validity of (ideal) hydrodynamics bypasses 
the necessity of the coupled kinetic equations: quarks and gluons are in local thermal equilibrium 
at all times. Therefore, within the framework of hydrodynamics we obtain the invariant photon 



production rate by directly replacing the initial quark distribution function n{q) in Eq. ( 3.18 ) by the 
"updated" distribution function n[q,T{T")] — l/[e^/-^'^'^ ^ + 1] at proper time t" > r^, where T{t") 
is determined by the cooling law Eq. (|6.12|) . The nonequilibrium invariant photon production rate 
that is consistent with the underlying hydrodynamics is then given by 

dN 



E- 



with 



[T,E,T{r)] = — ^ / dr" / 7^[c.,T(r")] cos[(c. - i?)(r - r")], (6.27) 



20 71^0- r d-^n r 

n^,T{r)] = __y _|_[2[l-(p.k)(p.q)]n[g,r(r)Mfc,r(r)]5(^ + 

+ [1 + (p • k) (p • q)] n[q, T{T)]n[k, T{t)]S{l^ -k-q)], (6.28) 



where n[q,T{T)] ~ 1/ [e'^/-'"*^'^^ + 1] and n — I — n. The momentum q integrals in Eq. (6.28) are 
calculated in the LR frame and hence are equivalent to those of the nonexpanding case above. 

Before proceeding further, we emphasize two noteworthy features of the nonequilibrium invariant 
photon production rate: (i) The photon production processes do not conserve energy, (ii) The rate 
depends on (proper) time not only implicitly through the local temperature but also explicitly. Fur- 



thermore, this explicit (proper) time dependence is non-Markovian as clearly displayed in Eq. (6.27). 
Whereas these two features seem rather uncommon in transport phenomena in heavy ion collisions 
(see, however, Refs. | |7^ , |12^ ), they are not unusual in nonrelativistic many-body quantum kinetics 
under extreme conditions [ p9| . In particular, energy- nonconscrving transitions and memory effects 
which cannot be explained by usual (semiclassical) Boltzmann kinetics have been observed recently 
in ultrafast spectroscopy of semiconductors with femtosecond laser pulses . 

The real-time kinetic approach when incorporated consistently with the hydrodynamical evo- 
lution of the QGP reveals clearly that photon production is inherently a nonequilibrium quantum 
effect associated with the expansion and finite lifetime of the QGP. 

At RHIC and LHC energies the quark distribution function T(t)] depends on the proper 
time T very weakly through the temperature T(r) within the lifetime of the QGP phase, hence a 
Markovian approximation (MA) in which the temperature in TZ[lu, T{t")] is taken at the upper limit 
of the integral is reasonable and hence the memory kernel may be simplified. In this Markovian 
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approximation, TZ[uj,T{t")] in Eq. (3.27) is replaced by TZ[uj,T{t)] and taken outside of the r 
integral. Thus Eq. ( |6.27 ) becomes 

dN 



E 



LR (2")' J-^ ' ' "(^ - E) 

A computational advantage of this Markovian nonequilibrium production rate is that it provides 
the "updated" quark distribution functions locally in time. Physically, the motivation for this 
approximation is that the most important aspect of the nonequilibrium effect is the nonconservation 
of energy originated in the finite lifetime of the QGP, a feature that is missed by the usual 5-matrix 
calculation, while the proper time variation of the temperature is a secondary effect and accounted 
for in the S'-matrix approach. 

It is worth noting that a connection with the Boltzmann approximation can be obtained by 
assuming completed collisions, i.e., taking the limit — oo in the argument of the sine function 



in Eq. ( p.2g| ) and using the approximation given by Eq. (6.2C). Consequently, the lowest order 
nonequilibrium photon production rate vanishes in the Boltzmann approximation due to kinematics. 



This highlights, once again, that the usual approach to photon production outlined in Sec. 6.2 and 
used in the literature corresponds to the Boltzmann approximation and therefore fails to capture 
the energy-nonconserving and memory effects that occur during the transient stage of evolution of 
the QGP. 

We are now in a position to calculate the direct photon yield from a longitudinally expanding 
QGP. In the CM frame the invariant production rate for photons of four-momentum from a fluid 



cell with four-velocity u'^ can be obtained from Eq. (6.29) through the replacement E P'^u^. In 
terms of the photon transverse momentum px and rapidity y, one finds P^u^^ = px cosh(?/ — rj). The 
photon yield is obtained by integrating the nonequilibrium rate over the space-time evolution of the 
expanding QGP. Assuming a central collision of identical nuclei, in the Markovian approximation 
we find the invariant nonequilibrium photon yield to lowest order in perturbation theory to be given 

by 



dN 



(PpT dy 



MA ^.Tf ^>,,„„ 

= TrR\ I dTT I drjE- 

CM Jri J -ri^on 



MA 

[T,pt^u^,TiT)], (6.30) 

LR 



d^pd'^x 

where Ra is the radius of the nuclei and — ?7con < rj < ?7ccn denotes the central rapidity region within 
which Bjorken's hydrodynamical model is valid. 

As remarked above, at RHIC and LHC energies the quark distribution function n[q, T{t)] depends 
on the proper time r very weakly through the temperature T(r) within the lifetime of the QGP phase, 
therefore the resultant photon yield is expected to qualitatively resemble the nonequilibrium photon 
yield from a nonexpanding QGP studied in above. This will be numerically verified below. We 
note that the expanding and nonexpanding cases differ mainly by the Jacobian r in the t integral 



in Eq. (6.3C) that accounts for the longitudinal expansion of the QGP, and by the replacement 
E P^Ufi in the argument of the invariant photon production rate that accounts for the shift of 
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Figure 6.5: Comparison of nonequilibrium (solid) and equilibrium (dotted) photon yields at midra- 
pidity (y ~ 0) from a longitudinally expanding QGP at RHIC energies with initial conditions given 
by = 1 fm/c and Ti = 200 (top), 300 (bottom) MeV. 



the photon energy under local Lorentz boost . 

6.4.3 Numerical analysis: central collisions at RHIC and LHC energies 

We now perform a numerical analysis of the nonequilibrium photon yield and compare the results 
to the equilibrium one obtained from higher order equilibrium rate calculations given by Eqs. ( |6.9[ ) 
and (6.10). The nonequilibrium photon yield in the Markovian approximation given by Eq. ( |6.30 ) 



contains a four-dimensional integral that is performed numerically for the values of parameters of 
relevance at RHIC and LHC energies. 

For central ^^'^Au+^^'^Au colhsions at RHIC energies ^/s^n ^ 200 GeV, we take Ra ^ 1.2 A^/^ 



fm « 7 fm |43| and ?7con = 2. The initial thermalization time is taken to be = 1 fm/c [123, 124| 



the final proper time Tf is determined when the critical temperature for the quark-hadron transition 



is reached at T{Tf ) ~ 160 MeV and is obtained from the cooling law given by Eq. (3.12) for a given 
initial temperature at proper time r^. 

The nonequilibrium photon yield at midrapidity {y = 0) in the range of transverse momentum 
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Figure 6.6: Rapidity distribution of the nonequilibrium photon yield aX pT — 1 (sohd), 2 (dotted), 
and 3 (dashed) GeV for a longitudinal expanding QGP at RHIC energies with initial conditions 
given by = Ifm/c and Tj = 200 (top), 300 (bottom) MeV. The distribution is symmetric at 
2/ = 0. 



0.3 < < 5 GeV/c is shown on a log-log plot in Fig. 6.5 for initial temperatures T,; = 200 and 
300 MeV. For comparison we also plot the corresponding equilibrium yield obtained by integrating 
Eqs. (6.9) and ( |6.1C ) with the transformation to the CM frame as specified by Eq. (6.14), and using 



the value of [122, 123, 129| 



«.[T(r)] = 



67r 



(33-27V/)ln[8r(T)/Td' 



(6.31) 



with N j = 2. Furthermore, Fig. 6.6 shows the rapidity distribution of the nonequilibrium photon 



yield at different values of px for Ti = 200 and 300 MeV. Several noteworthy features are gleaned 
from these figures: 

1. Whereas the equilibrium yield dominates the total yield for small pT, the nonequilibrium yield 
becomes significantly dominant in the range px > 1.0 — 1.5 GcV/c. Perhaps coincidentally, 
this is the range in which the CERN WA98 data for central Pb+Pb collisions at SPS energies 
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shows a distinct excess [119 



2. While the equihbrium yield leads to a transverse momentum distribution that falls off ap- 
proximately with the Boltzmann factor e~P'^/'^\ the noncquilibrium yield is not Boltzmann 
suppressed and falls off algebraically. It is found numerically that for px 2> Ti and within the 
region 1 ^ Pt ^ 5 GeV/c the noncquilibrium yield at midrapidity falls off with a power law 
with V ~ 2.47 and 2.77 for Ti = 200 and 300 MeV, respectively. This is a remarkable 
consequence of the fact that bremsstrahlung is the dominant process. As discussed in detail 
above at large energies the dominant process is bremsstrahlung corresponding to the contri- 
bution from the term n(q)[\ — n{k)] with fc = |p + q| in Eq. ( |6.1£ ). For large photon energy 



p the important contribution, which is not exponentially suppressed arises from the small q 
region. 

3. A central plateau in the range of rapidity y < 2 is seen clearly. The rapidity distribution begins 
to bend down when y ~ ?7ccnj i-e., when the photon rapidity probes the fragmentation region. 

4. The numerical analysis of the expanding case reveals features very similar to those found in the 
nonexpanding case, where we estimate that for pT > 1.5 GeV/c the higher order equilibrium 
contribution to the direct photon yield (which grows linearly in time) becomes of the same 
order as the lowest order noncquilibrium contribution (which grows at most logarithmically 
in time) only if the lifetime of the QGP is of the order longer than 1000 fm/c. Therefore, 
we conclude that these noncquilibrium effects dominate during the lifetime of the QGP in a 
realistic ultrarelativistic heavy ion collision experiment. 

Similar numerical analysis can be performed for central ^°*Pb+^"*Pb collisions at higher LHC 



energies ^/snn ^ 5500 GeV, for which we take Ra ~ 7 fm, r/ccn — 5 and Ti ~ \ fm/c [123, 124[. The 
comparison of noncquilibrium and equilibrium photon yields at midrapidity is displayed in Fig. 5.7 
for Ti — 450 and 500 MeV. The dominance of the noncquilibrium yield at high pT remains but 
now at higher transverse momentum pt ^ '2 GeV. This can be understood as a consequence of the 
longer QGP lifetime resulting from higher initial temperature at LHC energies. Furthermore, it is 
found numerically that for 1 ^ Pt ^ 5 GeV/c the noncquilibrium yield at midrapidity falls off with 
a power law p^'^ with i' ~ 2.52 and 2.56 for Ti — 450 and 500 MeV, respectively. 

These results indicate a clear manifestation of the noncquilibrium aspects of direct photon pro- 
duction associated with a transient QGP of finite lifetime. The most experimentally accessible signal 
of the noncquilibrium effects revealed by this analysis is the power law falloff of the transverse mo- 
mentum distribution for direct photons in the range 1 ^ Pt J; 5 GeV/c with an exponent 2.5 < < 3 
for temperatures expected at RHIC and LHC energies. Our numerical studies reveal that this ex- 
ponent increases with initial temperature and therefore with the initial energy density of the QGP 



and the total multiplicity rapidity distribution dNT^/dy [124[. This could be a clean experimental 
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Figure 6.7: Comparison of nonequilibrium (solid) and equilibrium (dotted) photon yields at midra- 
pidity (y ~ 0) from a longitudinally expanding QGP at LHC energies with initial conditions given 
by = 1 fm/c and Ti = 450 (top), 500 (bottom) MeV. 



nonequilibrium signature of a transient QGP since the photon distribution from the hadronic gas is 



expected to feature a Boltzmann type exponential suppression for px ^ Ti 1 112, 121 -123 1 



6.5 Conclusions 

In this chapter we study an important phenomenological application of the real-time quantum kinetic 
techniques developed in previous chapters. Our goal is to search for clear experimental signatures 
of direct photon production associated with nonequilibrium aspects of the transient QGP created in 
ultrarelativistic heavy ion collisions at RHIC and LHC energies. 

We first argue that the usual S'-matrix approach to direct photon production from an expanding 
nonequilibrium QGP has conceptual limitations. Instead we introduce a real-time kinetic approach 
that allows a consistent treatment of photon production from a transient nonequilibrium state of 
finite lifetime. 

We focus on obtaining the direct photon yield from a thermalized QGP undergoing Bjorken's 
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hydrodynamical expansion in central Au+Au collisions at RHIC {^/Tnn ~ 200 GeV) and LHC 
{^/snn ~ 5500 GeV) energies. The lifetime of a QGP for these collisions is of order 10 — 30 
fm/c. We find that energy-nonconserving (anti)quark bremsstrahlung q{q) — s- q{q)j and quark- 
antiquark annihilation qq — > 7, both of lowest order in a, dominate during such short time scales 
with bremsstrahlung dominating for pT S> T^. The contribution from these processes is a conse- 
quence of the transient nature and the finite lifetime of the QGP. As compared to the equilibrium rate 
calculations, the energy-nonconserving processes lead to a substantial enhancement in direct photon 
production for 1 — 2 < ^ 5 GeV/c near midrapidity. In striking contrast with the equilibrium 
calculation that predicts an exponential suppression of the transverse momentum distribution for 
PT S> Ti {Ti is the initial temperature), the nonequilibrium processes lead to a power law behavior 
instead. We find that at RHIC and LHC energies the direct photon transverse momentum distribu- 
tion near midrapidity is of the form with 2.5 ^ v ^ S for 1 ^ ^ 5 GeV/c and that photon 
rapidity distribution (for fixed px) is almost flat in the interval \y\ < ?7ccni where \ri\ < y^ccn denotes 
the central rapidity region of the QGP. The exponent v is numerically found to increase with the 
initial temperature, hence increases with the total multiplicity rapidity distribution dN-j^/dy, which 
is an experimental observable. 

Thus, as the main conclusion, we propose that direct photons are distinct experimental signatures 
of the transient nonequilibrium QGP created at RHIC and LHC energies, both in the form of a large 
enhancement at 1 — 2 < pT ^ 5 GeV/c as well as a power law transverse momentum distribution 
p^'^ with an exponent v that is within the range 2.5 ^ ^ 3 and increases with total multiplicity 
rapidity distribution dN^^/dy. 

Our study of direct photon production focuses solely on the QGP phase and neglects contribu- 
tions from pre-equilibrium stage as well as the mixed and hadronic phases, because we focus on a 
comparison between the usual equilibrium approach and the kinetic real-time approach that allows 
for following the time evolution of the density matrix as an initial value problem consistent with 
hydrodynamics. Furthermore, as discussed elsewhere in the literature, since most of the high-p^ pho- 
tons originate from the very early, hot stage in the QGP phase 110, 112, 114| and photons produced 
in the mixed phase as well as the subsequent hadronic phase are mainly in the lower-py region |123[ , 
the nonequilibrium yield from the QGP phase contributes dominantly to the total high-p^ photons. 

We have provided a systematic real-time description compatible with the initial value problem 
associated with hydrodynamic evolution, however more needs to be understood for a complete 
description of all the different stages. The parton cascade approach to describe the early pre- 
equilibrium stage after the collision is an important first step in a full microscopic description, but 
perhaps a more consistent description of the initial condition for QGP formation must be based on 



the recent notions of a color glass condensate [131, 132 1. Hence, a complete treatment of the direct 
photon yield must, in principle, begin from the initial stage, possibly a color glass condensate, and 
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obtain the real-time evolution of photon production. Clearly there must be many more advances in 
this field before such a program becomes feasible. The finite-temperature equilibrium calculations 
assume that the equilibrium thermal state always prevailed, thus ignores completely not only the 
initial stages but also the time evolution. Our approach while incorporating the time evolution 
consistently, also neglects the initial stage. We have focused on direct photons from a transient 
QGP for a direct comparison with equilibrium calculations, but obviously photons will continue to 
be produced during the mixed and hadronic phases. A detailed understanding of the transition as 
well as the hadronic photon production matrix elements is necessary for a more reliable estimate of 
the potential noncquilibrium effects after hadronization. 

Although we do not claim to have provided a completely detailed understanding of potential 
noncquilibrium effects due to the lack of knowledge of initial conditions in heavy ion collisions, 
we do claim that our approach provides a more systematic description of the finite-lifetime effects 
associated with a transient QGP. These effects lead to a distinct experimental prediction: a power law 
falloff of the distribution dN/(PpTdy near the central rapidity region which is distinguishable from a 
thermal tail for 1 — 2 < 5 GeV/c. We note that the parton-cascade results of photon production 

from the pre-equilibrium stage |126| also seem to reveal a power law falloff of the photon distribution 
in this range of transverse momentum with comparable order of magnitude and exponent, as can 
be gleaned from Figs. 2 and 3 of Ref. |126[ . Perhaps these two effects, namely the prompt photons 
from the pre-equilibrium stage and the direct photons from the hydrodynamically expanding QGP 
with a finite lifetime cannot be distinguished experimentally. Nevertheless, we emphasize that a 
power law departure from a thermal tail in the direct photon spectrum may very well be explained 
by noncquilibrium effects, either by a finite QGP lifetime as advocated in this chapter or by prompt 
photons from the pre-equilibrium stage. 

An important and very relevant question is that why not treat high energy particle collisions in 
the same manner, i.e., with the real-time evolution rather than with the S-matrix approach. The 
answer to this question hinges on the issue of time scales. In a typical high energy particle collider 
experiment, the colliding "beams" are actually bunches or packets with a typical spatial extent of 
order 1 — 10 cm and hence the typical time interval for the collision is of order 10^^ sec. This time 
scale is many orders of magnitude longer than the typical hadronic interaction time scale ^ 10^^'^ 
sec, thus taking the infinite time limit is amply justified. This of course is the basis for using the 
S'-matrix calculation in terms of asymptotic in and out states (at t — ^c>o): the total interaction 
time is much, much longer than the typical hadronic time scale. With regard to the heavy ion 
collision, the consensus is that after an initial pre-equilibrium stage, a locally thermalized QGP 
is formed. It then evolves hydrodynamically, hadronizes and eventually the freeze-out of hadronic 
gas ensues. Current theoretical estimates indicate a total time between formation and freeze-out 
of order 50 — 100 fm/c, with a QGP phase lasting for about 10 fm/c. These time scales are not 
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several orders of magnitude larger than the hadronic interaction time scale, thus the infinite time 
limit taken in the 5-matrix calculation is at best questionable. Therefore while in typical particle 
collider experiments the 5-matrix approach is valid (as has been demonstrated by half a century 
of experiments), the transient QGP with a finite lifetime of the order of the hadronic time scales 
merits a different treatment. 



As many recent investigations [130| have suggested that the QGP produced at RHIC and LHC 
energies is not expected to be in local chemical equilibrium, i.e., the distribution functions of quarks 
and gluons will probably be undersaturated, an important extension of our work will be the study of 
nonequilibrium effects on direct photon production from a chemically nonequilibrated QGP. Another 
important aspect that requires further investigation is the finite size of the QGP. Much in the same 
way as the finite lifetime introduces nonequilibrium effects associated with energy nonconserving 
transitions, we expect that the finite size ^ 7 fm of the QGP will introduce uncertainty in the 
momentum of the emitted particles. This is clearly an important topic that deserves further and 
deeper study but is beyond the scope of this thesis. However, at this stage we speculate that if such 
effects are present they could bear an imprint in transverse flow pq] . 



Chapter 7 



Summary 



The central theme of this thesis is the study of real-time nonequilibrium quantum dynamics with 
specific focus on the real-time relaxation of the mean fields and quantum kinetics in a variety of 
physically relevant quantum field theories with applications to cosmology and ultrarelativistic heavy 
ion collisions. This study could be generalized to nonequilibrium problems in condensed matter 
physics, quantum optics and atomic physics. To conclude this thesis, we summarize our work and 
highlight the main results. 

Fermion dcimping in a fermion-scalar plasma. We study the real-time relaxation of the 
fermion mean field in a fermion-scalar plasma which is relevant to models of baryogenesis at the 
electroweak scale. We begin by obtaining directly in real time the effective in-medium Dirac equation 
for the fermion mean field induced by an adiabatically switched-on external source. The in-medium 
Dirac equation is fully renormalized, retarded and causal, thereby leading to an initial value problem 
for the fermion relaxation. 

For the light scalar which cannot decay into fermion pairs, we find that to lowest order in the 
Yukawa coupling the only medium effect on fermions is dispersive and in-medium propagation of 
fermion excitations is undamped. On the other hand, for the heavy scalar such that its decay into 
fermion pairs is kinematically allowed, we find a novel dissipative medium effect arises from scalar 
decay in the medium and leads to damping of fermion excitations. That is, the fermions acquire a 
width due to the decay of the heavy scalar in the medium and become quasiparticles with a finite 
lifetime. Solving the effective Dirac equation by Laplace transform, we obtained the time evolution 
of the fermion mean field which allows a clear identification of the damping rate. We calculate 
the damping rate in the narrow width approximation to one loop order for arbitrary values of the 
fermion and scalar masses (provided the scalar is heavy enough to decay), temperature and fermion 
momentum. An all-order expression for the damping rate in terms of the exact quasiparticle wave 
functions is established by analyzing the structure of the fermion self-energy. 
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Finally, a Boltzmann-type kinetic equation for the fermion distribution function in the relaxation 
time approximation confirms the damping of fermion excitations as a consequence of the decay of 
heavy scalars in the medium, and displays the relationship between the fermion damping rate and 
interaction rate to lowest order in the Yukawa coupling directly in real time. 

Dynamical renormalization group approach to quantum kinetics. We provide an al- 
ternative first-principle derivation of quantum kinetic equations in quantum field theory by im- 
plementing a diagrammatic perturbative expansion improved by a resummation via the dynamical 
renormalization group directly in real time. 

This method begins by obtaining the equation of motion for the quasiparticle distribution func- 
tion in perturbation theory. The solution of this equation of motion at large times reveals secular 
terms that grow with time, as a consequence of the fact that the perturbative analysis breaks down 
at large times. We implement the dynamical renormalization group to resum these secular terms 
directly in real time, the corresponding dynamical renormalization group equation, which describes 
the time evolution of the quasiparticle distribution function insensitive to the physics on microscopic 
time scales, is interpreted as the quantum kinetic equation. A remarkable feature of this method 
is that it allows to include consistently medium effects via resummation akin to the hard thermal 
loops but away from equilibrium. 

We establish a close relationship between the dynamical renormalization group approach and the 
renormalization group in Euclidean field theory. In particular, coarse graining, stationary solutions, 
relaxation time approximation and relaxation rates have a natural parallel as irrelevant operators, 
fixed points, linearization near the fixed point and stability exponents, respectively, in the usual 
Euclidean renormalization group. This identification brings a new and rather different perspective 
to kinetics and relaxation that will hopefully lead to a new insights. 

We then apply this novel method to study quantum kinetics of pions and sigma mesons in the 
0(4) chiral linear sigma model. We derive the quantum kinetic equations that describes the time 
evolution of the pion and sigma meson distribution functions. The pion and sigma meson relaxation 
rates calculated in relaxation time approximation and agree with the corresponding damping rates 
found in the literature. We also find that in large momentum limit the emission and absorption of 
massless pions result in threshold infrared divergence in sigma meson relaxation rate and lead to a 
crossover behavior in relaxation that is unambiguously captured by the dynamical renormalization 
group analysis. 

Furthermore by establishing a direct correspondence between pinch singularities and secular 
terms, we show that the dynamical renormalization group approach to quantum kinetics provides a 
natural framework to interpret and resolve the issue of pinch singularities in nonequilibrium quantum 
field theory. 

Nonequilibrium dynamics in a QED plasma at high temperature. We study real-time 
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nonequilibrium dynamics in a QED plasma at high temperature within the hard thermal loop (HTL) 
approximation by implementing the dynamical renormalization group. The goal is to understand 
the relaxation of photon and fermion mean fields directly in real time and quantum kinetics of the 
photon and fermion distribution functions beyond the usual Boltzmann description of kinetics. As 
many features of a hot QCD plasma are similar to those of a hot QED plasma, in addition to its 
direct impacts on astrophysics and cosmology, this study is also of phenomenological importance in 
ultrarelativistic heavy ion collisions. In particular, the leading contributions in the hard thermal 
loop approximation can be generalized easily from hot QED to hot QCD, thus relaxation of quarks 
and production of photons in the QGP can be understood to leading order from the study of a hot 
QED plasma. 

For semihard photons of momentum eT <C fc <C T, we find to leading order in the HTL approxi- 
mation that the photon mean field relaxes with a power law ~ f-e^T^/i2fc^ large times t 1/k, as 
a result of infrared enhancement of the photon spectral density near the Landau damping threshold. 
The quantum kinetic equation for the distribution function of semihard photons, which includes 
off-shell effects missed by the usual Boltzmann description, is derived directly in real time from 
first principles in the relaxation time approximation by using the dynamical renormalization group 
method. We find that the distribution function also relaxes with a power law, with an exponent 
which is twice that for the photon mean field. The dynamical renormalization group analysis reveals 
the emergence of detailed balance on microscopic time scales larger than 1/k while the rates are still 
varying with time. 

We show directly in real time that as a consequence of the emission and absorption of soft, quasi- 
static magnetic photons, hard fermion mean field with momentum fc ~ T exhibits a nonexponential 
relaxation at large times as ~ e-o:Tt{\nii>pt+o.i26...) ^ where cjp = eT/3 is the plasma frequency and 
a = /A-K. A quantum kinetic equation for the distribution function of hard fermions in the 
relaxation time approximation is obtained directly in real time by implementing the dynamical 
renormalization group resummation. Unlike that of the usual Boltzmann kinetic equation derived 
in the quasiparticle approximation, the collision term of this quantum kinetic equation is time- 
dependent and infrared finite. The distribution function is found to relax nonexponentially with an 
anomalous exponent twice as large as that for the hard fermion mean field. 

Our real-time analysis reveals clearly the limitation of the usual Boltzmann description of kinetics 
and the quasiparticle approximation in the study of nonequilibrium quantum dynamics in gauge field 
theory and has important impacts on direct photon production from the quark-gluon plasma. 

Direct photon production from the quark-gluon plasma. As an important phenomeno- 
logical application of real-time nonequilibrium dynamics in quantum field theory, we study direct 
photon production from a longitudinally expanding QGP at RHIC and LHC energies. Instead of 
focusing on obtaining an equilibrium photon production rate using the usual 5-matrix formalism, we 



128 



study photon production as an initial value problem with the real-time quantum kinetic description 
developed in this work. 

This real-time approach allows us to incorporate consistently the kinetic description of photons 
with the hydrodynamical description of the QGP. Within Bjorken's hydro dynamical model, we show 
that the lowest order energy-nonconserving (anti)quark bremsstrahlung q{q) — > q{q)'j and quark- 
antiquark annihilation 9? — > 7 contribute dominantly to photon production during the transient 
hfetime of the QGP. For central collisions at RHIC (y^sTTiv ~ 200 GeV) and LHC {^/sWN ~ 5500 
GeV) energies, we find a significant excess of direct photons in the range of transverse momentum 
1 — 2 < ^ 5 GeV/c as compared to higher order equiUbrium estimates. The rapidity distribution 
of photons is fairly flat near the midrapidity region. The remarkable result is that the transverse 
momentum distribution of photons at midrapidity falls off with a power law p^" with 2.5 < < 3 
as a consequence of the energy-nonconserving processes taking place during the transient lifetime of 
the QGP, providing a distinct nonequilibrium experimental signature of the QGP formation. 

As summarized above, in this thesis we have advanced several aspects of real-time nonequilibrium 
dynamics in quantum field theory. In particular, we have estabUshed a theoretical framework for 
studying nonequilibrium dynamics of quantum multiparticle systems directly in real time from first 
principles. This framework is used to study a variety of physically relevant quantum field theories. 
It not only reproduces the usual results in the Hterature when the same approximations and assump- 
tions are invoked, but also leads us to a new territory which has not been explored before. Most 
importantly, it is capable of making concrete predictions that can be tested experimentally at BNL 
RHIC and CERN LHC. We envisage our work will also have direct impacts on the understanding 
of nonequilibrium phenomena in condensed matter physics, quantum optics and cosmology. 



Appendix A 



Renormalization Group and 
Asymptotic Analysis 



In this appendix we present a pedagogical introduction to the renormahzation group (RG) method 
in studying asymptotic analysis of ordinary differential equations. In order to introduce the basic 
idea, we consider a simple but illuminating example: a weakly damped harmonic oscillator. The 
reader is referred to Refs. for generalization to the much more complicated problems. 

The differential equation governing a weakly damped harmonic oscillator is given by 



y + y = -ey, e < 1 



(A.l) 



with initial conditions y{to) and {/(to) specified at time to, where the overdot denotes the derivative 
with respect to t. Note that Eq. ( [A.l|) can be solved exactly with the exact solution given by 



y{t) = Roe 



-e(t-to)/2 



COS 



i--(t-eo) 



(A.2) 



where Rq and Qo are constants determined by the initial conditions. However, in order to demon- 



strate the renormalization group method, we attempt to solve Eq. (A.l) in a perturbative expansion 
in powers of e 



y = yo + m + <^ 2/2 + 



(A.3) 



Upon substituting the perturbative expansion into Eq. (AJ), one obtains a hierarchy of equations: 



yo + ya = 0, 
yi + yi = -yo, 
2/2 + 2/2 = -yi, 



(A.4) 
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Starting from the solution to the zeroth-order equation 



yo{t) = Rocos{t - Go), 



(A.5) 



and the zeroth-order retarded Green's function 



G^{t-t') = sin{t~t')0{t-t'), 



(A.6) 



one obtains 

2/1 W = - / dt'sm{t-t')yoit') 

J to 

Ro 



2/2 (t) 



[{t - to) cos(< - 9o) - cos(to - 6o) sin(i - to)], 
dt' sm{t-t')yi{t') 



to 



Rq 



{t - to? + n cos{t - Qo) - (i - to) sm{t - 2to + Qo) - cos(t - 2to + Qo) 



(A.7) 



Hence, to order 0{e^) the solution is given by 



yit) = Ro< 



l-|(i-io) + ^(i-io)' 



cos(t - eo) - -(t - to) sin(t - 2to + Oq) 



+ I cos(to - Oo) sin(t ^ ^o) + ^ 



- cos(t - Oo) - cos(t - 2to + Qo) 



Oie'). 



(A.; 



It is noted that this solution contains secular terms that grow in t. This perturbative expansion 
breaks down for e(t — tg) > 1 because of these secular terms. The presence of secular terms in 
a perturbative expansion indicates that perturbation method is only suitable to study short time 
behavior of a nonequilibrium system. To extract information of the long time behavior, one needs to 



go beyond the perturbation method. A close inspection of the secular terms in Eq. (A. 8) suggests the 



prefactor multiplying cos(i — Oq) resembles a time-dependent amplitude R{t) and that multiplying 
sin(< — 2to + Oo) resembles a time-dependent phase 6(t). 

The way to solve this difficulty of the perturbation method is to take account of its limitation 
and hence apply it only to describe short time behavior. The perturbative solution is evolved in time 
from to only up to a time, say ti, such that e(ti — to) ^ 1 and hence the perturbative solution is still 
reliable. Then the amplitude i?(ti) and phase Q{ti) is used as the updated initial conditions at ti 
to iterate the perturbative solution forward in time to t2 say ti, such that e(t2 ~ ti) <^ 1. Repeating 
this procedure, one obtains an improved perturbative solution up to an arbitrary time r. The idea 
of renormalization group comes in when one wants to find out how the amplitude R{t) and phase 
6(r) change as r changes. 
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Let us assume that R{t) and 0(t) are related to Rq and Oq, respectively, by a multiplicative 
amplitude renormalization Zji(t) and an additive phase renormalization Zq(t) such that 



i?o = ^i?(r)i?(T), 60 - e(r) + Z0(t), 
with Zfilr) and Zq{t) being expanded in power of e as 

^fl,(T) = 1 + eai(T) + 6^02 (r) + • • • , Zeir) = e5i(r) + e^feaCr) + 



(A.9) 



(A.IO) 



Inserting the above expressions into Eq. (A. 8), one obtains to order 0{e) 



yit) 



R \ cos(i - e) - e 



^{t - to) - ai{T) 



cos(i - e) 



- [cos(e - to) sin(i - to) + 26i(r) sin(i - 6)] J> + O(e^). 



(A.U) 



Now we want to remove the term that grows with to- Because for a fixed (eventually large) t, this 
is a secular term. Choosing 

(A.12) 



aiir)^^iT-to), &i(t) = 0, 



we can rewrite Eq. (A.U) as 

y{t, t) = R{t) { [1 - |(i - ^)] cos[t - e(r)] + | cos[e(T) - to] sm{t - to)] + 0{e^). (A.13) 



A remarkable feature of Eq. ( A.13| ) is that it is valid for large times as long as (i) t is chosen 
such that e(i — r) ^ 1, and (ii) R{t) and Q(t) are uniquely determined by the initial conditions. 



The first condition is trivial as the variable t introduced through Eq. (A.9) is complete arbitrary. 
To determine R{t) and 0(t), we notice that since t done not appear in the original problem, the 
solution y{t, r) cannot depend on the arbitrary variable r. This leads to 



dr 



dR 



-R 



cosh - e(r)] + R^ sink - QM] = 0. 
or 



(A.14) 



dr ' 2' 

Using the fact that cos[t — 6(t)] and sin[t — Q{t)] are linearly independent functions, to this order 



in e, Eq. ( [A.14| ) reduces to 



^ = -|i?(r) + 0(a, ^^0 + O(a. 



(A.15) 



This set of equations are analogous to the renormalization group equations in Euclidean field theory 
in the sense that they describe the flow of the amplitude R{t) and phase 8(t) under the change of 
T. The time scale r can be thought of as an arbitrary renormalization scale, and the RG equation 
determine how the "renormalized" amplitude and phase depend on the renormalization scale in such 
a way that the change in the renormalization scale is compensated by the change in the renormalized 
amplitude and phase. It is easy to see that the configuration R{t) = and 0(r) = const corresponds 
to the fixed point of the RG equations. 
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Solving the above RG equations and setting t — t, one obtains 

R{t) = R{0)e-<'-'°^/^ + 0{eH), e{t) = 6(0) + Oieh), 



(A.16) 



where R{0) and 9(0) are constants determined by initial conditions at Iq. Upon setting t = t in 
Eq. ( A.13| ), we obtain an "improved" perturbative solution y{t) which to this order does not contain 
any secular terms and hence is valid for large times. With the initial conditions y{to) = 1 and 
2/ (to) — 0, one finds i?(0) = 1, 0(0) = to and the RG- improved solution reads 



e(t-to)/2 



cos(t - ^o) + 2 sin(< - to)\ + 0{e^). 



(A.17) 



This renormalization group method can be extended to higher order in e. Specifically, to order 
0{e^) the choice of renormalization constants 



a2(T)-i(r-to)', hir) 



1 



{r-to), 



(A.18) 



can be shown to completely remove all of the secular terms to this order. This is a sign of the 
renormalizability to all orders in perturbation theory. After amplitude and phase renormalization 
the solution is found to be given by 



Rir) 

^2 



1 



cos[t — 8(t)] + - cos[0(t) — Iq] sin(t — Iq) 



1 



{t - r) shi[t - 2ta + e(T)] + - cos[< - 2to + e(r)] 



+ 0(e3). 



The fact that y{t, r) is independent of r leads to 

cos[t - e(r)] + 



dR e 



de e dR 
dr 4 dr 



3in[t-e(T)] +0(e^) =0. 



Hence, the RG equations to order 0{e^) read 



dR 
d7 



The solution of the above RG equations (A. 21) is given by 



R{t) = i?(0) e-'(*-*«'/2 ^ oieH), 0(1) = 6(0) + + OieH), 



(A.19) 



(A.20) 



(A.21) 



(A.22) 



where R{0) and 6{0) are constants determined by initial conditions. Using the initial conditions 
specified above, to order e^, one finds the RG-improved solution is given by 



^ g-«(t-to)/2 



16 



1 H cos 

16 ' 



1 + ^ ) - io) 



I- J jit -to) 



e 

- cos 
2 



^{t-to 



sin(t — io) 



(A.23) 



The interpretation of the renormalization group method is very clear in this simple example: 
the perturbative expansion is carried out to a time scale ^ 1/e within which perturbation theory 
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is valid. The correction is recognized as a change in the ampUtude and phase, so at this time scale 
the correction is absorbed into the renormalization of the amplitude and phase. The perturbative 
expansion is carried out to a later time but in terms of the amplitude and phase at the renormalization 
time scale. The dynamical renormalization group equation is the differential form of this procedure 
of evolving in time, absorbing the corrections into the amplitude and phase, and continuing the 
evolution in terms of the renormalized amplitudes and phases. 
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